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PREFACE 


The objective of this contract (NAS5-24067) is the development of 
computer, software for the preflight mission analysis of missions to earth- 
sun libration points. This software, designated STEAP-L, extends the capa- 
bility of the Space Trajectories Error Analysis Programs (STEAP) developed 
under contracts NASl-9745, NAS5-11795, and NAS5-11S73 and begins the inte- 
gration of STEAP with the Goddard Trajectory Determination System (GTDS). 

The software produced consists of two related programs, both of which 
use the GTDS Cov7ell propagator for the computation of the trajectory and state 
transition matrices. The first program, NOMNAL, is responsible for the gen- 
eration of the nominal trajectory from launch at earth to insertion into halo 
orbit about the desired libration point. NOMNAL uses a Newton-Raphson 
iteration (moving backward in time from the insertion maneuver) to perform 
the targeting of both impulsive and finite burn insertions into halo orbit. 

A user-controlled launch profile allows the transfer to be tied to a realis- 
tic launch and injection. NOMNAL stores the targeted trajectory and state 
transition matrices on a file for later analysis by the second program 
ERRAN. 

The program ERRAN performs generalized linear error analyses along 
specific targeted trajectories. Knowledge and control covariances are pro- 
pagated along the trajectory through a series of measurements and guidance 
events in a totally integrated fashion. The knowledge covariance is pro- 
cessed through measurements using a Kalman-Schmidt recursive filter with 
arbitrary solve- for/consider/ignore state augmentation. Probabilistic 
midcourse corrections are computed using an exact analytic formulation. 

ERRAN obtains the trajectory and state transition matrices from a file 
generated by NOMNAL for program efficiency. 

A major conclusion of this effort is that the complementary features 
of the GTDS and STEAP systems may be effectively combined to yield a signifi- 
cantly improved system. Thus the Cowell file generator/reader capability of 
the GTDS has been combined with the generalized covariance analysis of STEAP 
to yield a more efficient, extended error analysis capability than either 
system had previously. Other conclusions reflect the efficacy of the back- 
ward targeting algorithm developed for the libration mission targeting and 
the analytic formulation implemented for the midcourse correction sizing. 

The general recommendations for future effort identified during this 
study are two-fold.. Because of the success of this preliminary integration 
of the GTDS and STEAP systems it is recommended that this effort be continued 
and enlarged. In the specific area of libration point mission analysis, it 
is recommended that more detailed models (e.g., pulsing thrust insertion into 
halo orbit) be developed and continued studies be made of critical problems - 
(e.g. station-keeping error analysis) for these peculiar missions which are 
neither interplanetary, lunar, nor earth-orbiting. 
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1 . INTRODUCTION 


This Programmer's Manual summarizes the structure and coding of the STEAP-L 
programs NOMNAL and ERRAN and the subroutines they include. The discussions 
are intended to provide the reader with sufficient information to effectively 
modify or extend those programs. An accompanying volume, the Analytic and 
User's Manual, summarizes the mathematical assumptions and analysis of the 
programs and details the actual usage (input and output requirements) of those 
programs. 

This volume is divided into three major parts. This introductory chapter 
discusses the general development of the STEAP library of programs, describes 
the libration point missions toward which the current effort is directed, and 
summarizes the capability of the two programs developed for this application, 
NOMNAL and ERRAN. The second and third chapters summarize the overall 
programming and storage of the NOMNAL and ERRAN programs respectively. The 
fourth chapter forming the bulk of this volume provides documentation of each 
subroutine of STEAP-L in alphabetical order. 

1.1 Development of STEAP 

/ 

STEAP is an acronym for Space Trajectory Error Analysis Programs. Rather 
than a single computer program, STEAP is a library of related programs for the 
analysis of the navigation and guidance charac ter ist ics of space missions. 

These programs have been developed, modified, and extended over a number of 
years by the Martin Marietta Corporation (MMC) under the direction of NASA in 
a variety of contracts. ' 

There are two primary unifying elements in the development of the STEAP 
system. The first is in the underlying philosophy of STEAP. STEAP has always 
been directed toward the performance of a totally-integrated analysis of the 
navigation and guidance processes of space missions. Thus interaction is 
continually forced between the tracking uncertainties and the maneuver 
execution errors to determine the evolving uncertainties in the knowledge and 
control of the spacecraft trajectory. The second element is in general program 
structure. The STEAP software has continually been divided into three distinct 
operational modes responsible for nominal trajectory targeting and generation 
(NOMNAL), linear error analyses (EREIAN), and single-case or Monte Carlo 
simulations (SIMUL) . . The current effort does not address the third of these 
types of programs. 

The mathematical foundation for the STEAP system was initially developed 
under Contract NAS8-21120 for Marshall Space Flight Center. The first version 
of STEAP (Contract NASl-87,45) was constructed for general interplanetary 
ballistic missions for Langley Research Center to support the Viking mission 
analysis and design. Later development of STEAP was performed for Goddard 
Space Flight Center (Contracts NAS5-11795 and NAS5-11873) where specific 
extensions required for Planetary Explorer (later known as Pioneer Venus) and 
general lunar missions were added in a version called STEAP-II. More recently, 
programs for the navigation and guidance analysis of low thrust inter-planetary 
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and near-Earth missions have been developed for Langley Research Center (NASl- 
11686) and Marshall Space Flight Center (Contract NAS8-29666) . Throughout 
this time, improvements in the analytical techniques and program structure 
have been continually identified and incorporated into the STEAP series of 
programs. (References 1-5). 

Under the current contractual effort, versions of NOMNAL and ERRAN 
appropriate for missions to Earth-Sun libration points have been developed 
(termed STEAP-L) . A very significant feature of this effort is that the 
Goddard Trajectory Determination System (GTDS) Cowell propagator is being 
integrated into the STEAP-L programs. The Cowell propagator permits the 
generation of a file containing trajectory and state transition matrix 
(computed by integration of the variational equations) data during the NOMNAL 
run. This data may then be efficiently retrieved in subsequent ERRAN runs, 
thereby eliminating the costly integration cycle from ERRAN. 

A number of new analytical features have been added to STEAP under this 
contract. An unusual approach has been used in the targeting of the libration 
point missions. Backward integration is used in computing the successive 
trajectory iterates and targeting matrices required by the Newton Raphson 
targeting algorithm. This backward targeting scheme efficiently produces a 
targeted transfer trajectory that is consistent with realistic launch and 
injection constraints. The approach is well-suited to cometary or lunar 
missions as well. 

An exact computation of the probabilistic midcourse correction require- 
ments using the recently published technique of Lee-Boain (Reference 6) has been 
added to ERRAN. This replaces the previous model which employed the Hoffman- 
Young approximation (Reference 7) and which could lead to significant errors 
at the higher probability levels. This technique is applicable to lunar or 
interplanetary trajectories as well as the libration point missions.' 

A third significant item developed during this effort has been the re- 
formulation of the variable time-of-arr ival (VTA) guidance policy for the 
libration point mission application. The guidance policies available in 
previous versions of STEAP always assumed that the target state was referenced 
to a gravitational body such as the moon or a planet. This restriction has 
now been removed. 

The characteristics of the libration point missions necessitating these 
extensions are described in the summary of the libration point missions given 
in the next section. The capabilities of the resulting programs NOMNAL and 
ERRAN are then detailed in the next two sections. 

1.2 Libration Point Mission Application 

The STEAP-L programs developed under this contract are designed for use 
primarily for the analysis of missions to the two Earth-Sun libration points 
near the Earth. These are designated Lj^ and L 2 in Figure 1.1, which shows 
schematically the location of all five classical Lagrangian or libration points. 
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Figure 1.2 Details of and L 2 
Libration. Points 


Figure 1.2 shows in more detail the location of points and L 2 with 'respect 
to the Earth, with the orbit of the Moon, the classical or Laplacian sphere 
of influence of the earth, and an enlarged version occasionally used in tar- 
geting of swingby missions. The two spheres of influence are defined by 

2/5 
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where Rgg is the Earth-Sun distance and Mg and Mg are the masses of the Earth 
and Sun respectively. 


Figure 1.1 Earth-Sun Libration 
Points 



Efficient transfers from circular Earth parking orbit to the L-j^ and L 2 
points have been sho^'m (Reference 8) to fall into at least two major families; 
those with short (~25 to 50 day) transfer times and those with long (-100 to 
135 day) transfer times. The fast transfers require from 341 to about 400 
meters/second AV to insert into orbit near the libration point, with the 
minimum AV at about 36.4 days. The slow transfers require insertion AV of 
from 272 to about 400 meters/second, with the minimum AV at about 116.8 days. 
These optimum insertion values are based upon the Earth in a circular orbit 
around the Sun and will vary slightly due to the ellipticity of the orbit of 
the Earth,. The influence of the moon will affect them also. Both of the 
families discussed above assume a posigrade transfer orbit upon leaving the 
Earth; corresponding families exist for retrograde departures, but these re- 
quire higher insertion AV at the libration point. For long flight times at 
least two other families of trajectories exist but have higher AV require- 
ments. Even more families exist with longer flight times (-175 days) that 
have lower AV requirements ('200 meters/ second) (Reference 8) . 
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The primary feature of the libration points is that they are equilibrium 
points of the system; i.e., if a spacecraft is placed exactly at a libration 
point with no motion relative to the system, it will remain at that point 
relative to the two-body configuration. The collinear points L 2 , L 3 ) 

are unstable while the equilateral triangle points (L^, L 5 ) are only quasi- 
stable. Thus, some form of station-keeping is necessary to maintain the 
spacecraft in that location. However, the fuel required is still much less 
than it would be. at arbitrary points of the system. Thus, the and L 2 
points offer attractive stations for spacecraft for monitoring solar or solar/ 
earth phenomena- (Reference 9). To facilitate communications, the spacecraft 
would generally be placed in a "halo -orb it” about the libration point so that 
the sun would not obstruct the view of the spacecraft from earth. A typical 
halo-orbit in the plane normal to the rotating earth-sun line is illustrated 
in Figure 1.3. 



X-axis from Sun to Earth 

Z-axis Normal to Ecliptic 

Y-axis Completes Right Hand 
System 


Figure 1,3 Typical Halo-Orbit as Viewed from Earth 


The current effort is directed toward the study of the transfer and in- 
sertion phases of the libration point mission; the station-keeping while in 
the halo-orbit was not addressed in this effort. The two programs developed 
for the analysis of libration point transfers include the nominal trajectory 
and maneuver targeting program I\'0>u'JAL and the navigation and guidance error 
analysis program ERRAN sutranarixed in the next two sections. 

1.3 Summary of NOMNAL 

The computer program NOIiNAL is responsible for the generation of a nomi- 
nal trajectory from injection at earth to insertion into a halo orbit about 
a libration point in the earth-sun system. 

NOMa'IAL uses a specialized version of the GTDS Cowell propagator for the 
integration of the trajectory equations. The dynamic model used in the re- 
duced Cowell propagator includes the accelerations, on the spacecraft pro- 
duced by a central body, up to two non-central bodies, and finite thrust 
engines. The Cowell propagator generates state and control transition ma- 
trices by integration of variational equations simultaneously with the equa- 
tions of motion. These matrices are then used in the targeting of the li- 
bration point missions within NO'-T’AL and in the propagation of covariance 
matrices and the error analysis of the finite burn insertion maneuver in ERRAN, 
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NOMNAL has the capability to target transfer trajectories to libratlon 
points using both impulsive and finite thrust insertion maneuvers. In either 
case a backward targeting scheme is employed where conditions at the libra- 
tion point are iteratively improved to yield trajectories which when propa- 
gated backwards in time from the desired arrival point and time to the earth 
satisfy desired target conditions. The three target conditions at the earth 
are radius of closest approach, equatorial inclination at closest approach, 
and time at closest approach. These three conditions are normally selected 
to be consistent with the desired parking orbit radius, launch site latitude, 
and desired trip time. 

In impulsive targeting the three controls at the libration point are the 
three components of velocity on the transfer trajectory. In finite thrust 
targeting the controls are the right ascension and declination of the thrust 
direction and the duration of the burn; the thrusf magnitude, engine specific 
impulse, and initial spacecraft mass are held constant at the user-supplied 
values, A Newton-Raphson algorithm is used to iteratively improve the con- 
trol parameters to determine their required values. 

The program includes three options for the determination of the zero 
iterate values to begin the targeting process; table interrogation, conic 
approximation, and user-specification. Tables defining targeted velocities 
have been constructed for transfers to the and L 2 points with trip times 
in the vicinity of either optimal transfer (tabulated AVs for trip times of 
from 25 to 50 days and from 102 to 130 days at 1 day intervals). Initial 
values of velocity may then be interpolated from the data stored In these 
tables. The second option computes the initial libration point velocity by 
solving Lambert's theorem for the geocentric conic connecting the libration 
point radius and the injection radius in the desired time. The third option 
accepts a user-supplied zero iterate vector computed by the user outside the 
program. 

NOMNAL can adjust the injection time of the transfer to correspond to a 
realistic launch profile specified by the user. It then adjusts the arrival 
time by the same amount to hold the trip time at the user-desired value, 

NOMNAL computes and records such information as the required launch azimuth, 
coast time, and whether or not a coplanar injection maneuver is required. 

1.4 Summary of ERRAN 

The error analysis/generalized covariance analysis program ERRAN is a 
preflight mission analysis tool that is used to determine how selected error 
sources influence the orbit determination process for libration point missions. 

In the error analysis mode, ERRAN provides three primary quantitative 
results: (1) knowledge covariance matrices, which provide a measure of how 
well the actual trajectory is known, (2) control covariance matrices, which 
when propagated forward to the target provide a measure of how well the 
nominal target conditions will be satisfied by the actual trajectory, and 
(3) statistical midcourse AVs, which provide a measure of the amount of fuel 
required for a successful mission. 



In the generalized covariance analysis mode, ERRM provides all of the 
above information plus corresponding ’’actual" statistical information. The 
three results discussed in the previous paragraph are all computed on the 
basis of statistical distributions assumed by the navigation filter to des- 
cribe the significant error sources. In the generalized covariance analysis 
mode, "actual" knowledge covariances, control covariances, and statistical 
midcourse AVs are computed on the basis of statistical distributions that 
actually describe both error sources acknowledged by the navigation filter 
and the error sources ignored. The primary use of the generalized covariance 
analysis program is to study the sensitivity of filter performance to off- 
design conditions. 

Up to 15 measurement parameters may be solved-for or considered by the 
navigation filter employing a Kalman-Schmidt sequential formulation. Param- 
eters not acknowledged in design of the filter may be treated as ignore pa- 
rameters vjhen ERRAN is run in the generalized covariance analysis mode. 
Measurement biases include biases in the locations of the three earth-based 
tracking stations, and biases in all measurements. Available measurement 
types are range, Doppler, and a simple optical model. Measurement noise 
for each measurement type is assumed to be constant. 

The computational procedure in ERRAN is divided into basic cycle com- 
putations and event computations. Basic cycle computations are concerned with 
the propagation of covariances fon-rard to a measurement time and processing 
the measurement. Events refer to a set of specialized computations, not 
directly concerned with measurement processing, that can be scheduled to occur 
at arbitrary times along the trajectory. State transition matrices Interpo- 
lated from the file created by NOMNAL are used for all covariance matrix 
propagation. 

The four events available in ERRAN are eigenvector, prediction, guid- 
ance, and final insertion into halo orbit. At an eigenvector event the 
position and velocity partitions of the knowledge co\^arlance matrix are 
diagonalized to reveal geometric information about the size and orientation 
of the position and velocity navigation uncertainties. At a prediction event 
the most recent covariance matrix is propagated forward to some critical tra- 
jectory time to determine predicted navigation uncertainties in the absence 
of further measurements. 

The guidance event is the most complex event and yields much useful 
information for preflight mission analysis. Several types of guidance events 
are available in ERRAN. At a midcourse guidance event the user can choose 
from cither fixed or variable time of arrival guidance policies (FTA or VTA), 
Execution error statistics are generated using an impulsive error model de- 
fined by a proportionality error, a resolution error, and two pointing angle 
errors. The execution errors of the insertion maneuver may be modeled as 
either an impulsive maneuver (defined above) or a finite thrust maneuver 
(component errors modeled as two pointing errors and a thrust magnitude 
uncertainty). The target condition covariance matrix both before and 
after the maneuver is printed out for midcourse and insertion maneuvers. 
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2c nomnal program structure 


2c 1 nomnal program description 

The steap-l/nomnal program is a digital computer program 
WRITTEN IN fortran IV COMPATIBLE WITH THE IBM 360/370 SYSTEM. 

THE program CONTAINS ABOUT EIGHTY SUBROUTINES OF WHICH ONE-HALF 
are associated with THE COWELL PROPAGATOR* ONE THIRD ARE GENERAL 
utility routines* and THE REMAINDER ARE SPECIFIC ROUTINES FOR THE 
LIBRATIon POINT APPLICATION. THE PROGRAM ACCESSES FOUR DATA SETS 
DURING operations THE STANDARD FORTRAN INPUT (SYSIN) AND THE 
presentation output (SYS0UT=A) data SETS* THE DIRECT ACCESS 
solar/lunar/planetary ephemeris file used by the gtos 

{DSN=6To5.SLP1950o JAN71 ) AND THE SEQUENTIAL ORBIT 'FILE 
ON WHICH THE TRAJECTORY AND STATE TRANSITION MATRIX DATA 
are stored for later USE® the program <N0N-0VERLAID) 

REQUIRES 255K OR 2A5K BYTES STORAGE DEPENDING UPON 
WHETHER THE ORBIT FILE IS GENERATED OR NOT, 


2.2 NOMNAL SUBROUTINE HIERARCHY 

figure 2®1 illustrates the GENERAL PROGRAM STRUCTURE. THE 
MAIM PROGRAM CALLS THREE EXECUTIVE ROUTINES FOR THE THREE MAJOR 
ACTIVITIES OF THE PROGRAM, HPRELM IS RESPONSIBLE FOR PRELIMINARY 
DATA MANIPULATION AND THE GENERATION OF THE ZERO ITERATE SOLUTION, 
MGIDNS is RESPONSIBLE FOR THE TARGETING OF EITHER THE IMPULSIVE 
OR FINITE BURN INSERTION TRAJECTORY, HTRJTY PROPAGATES AND 
STORES THE REQUIRED TRAJECTORIES AND STATE TRANSITION MATRICES 
USING The GTDS COWELL PROPAGATOR SPECIALIZED TO THIS APPLICATION. 
CHAPTER 4 PROVIDES AN INDEX OF ALL ROUTINES USED IN NOMNAL AND 

GIVES Individual routine documentation including descriptions* 
analyses* and Flowcharts. 


2.3 NOMNAL COMMON VARIABLES 

the following two subsections provide definitions OF THE 
variables appearing in common Blocks in the nomnal programs, 
subsection 2.3.1 LISTS THE COMMON BLOCKS IN ALPHABETICAL ORDER AND 
GIVES The SIZE OF THE BLOCKS AND THE DEFINITIONS OF THE VARIABLES 
appearing within each BLOCK. SUBSECTION 2.3.2 LISTS ALL THE 
COMMON variables IN ALPHABETICAL ORDER AND DEFINES THE BLOCK TO 

WHICH They belong and their definition. 
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Figure 2. 1 NOMNAL SUBROUTINE HIERARCHY 












8.3.1 common blocks IN ALPHABETICAL ORDER 


NAMEIDIM) 


DISP 

DEFINITION 

BDATA {SIZE 
N0OO 

10) 

0 

number of gravitating BODIES 

N8<3) 

R 


4 

VECTOR OF BODIES (STEAP CONVENTION) 

FLAGSl (SIZE 
IBODY 

10) 

0 

CENTRAL BODY NUMBER (GTOS CONVENTION) 

INDXA 


4 

NOT USED 

INDYA 


8 

NOT USED 

NEQ 


C 

number of variational EQUATIONS 

FLA6S2 (SIZE 

E8) 



ICENTB 


0 

CENTRAL BODY NUMBER (6TDS CONVENTION) 

NT5F(5S 


4 

total number of SECTIONS 

NSECTN 


8 

current SECTION NUMBER 

INDSECnOfS) 


C 

flags for section 

INDIES) 


84 

flags for current section 

FLAGS3 (SIZE 

30) 



NCNM 


0 

NUMBER OF C(N.M) TO BE ESTIMATED 

NSNM 


4 

NUMBER OF S(N.M) TO BE ESTIMATED 

N6 


8 

TOTAL NUMBRE OF POTENTIAL COEFFICIENTS TO BE 

NSTftTE 


C 

ESTIMATED 

number of STATE UNKNOWNS 

KSTATE(6) 


10 

VECTOR OF LABEL NUMBERS FOR STATE UNKNOWNS 

NCONDT 


88 

condition NUMBER 

HANOON 


2C 

MANEUVER OCCURENCE FLAG 

FLA6S4 (SIZE 

») 


! 

NCURS 


0 

gtds parameter 

NENTRV 

V • V S A «> CB 0 0 0 O 0 ^ 0 ■ 


4 

gtds parameter ^ 


) 

2-3 



NAMEJOIMI 

FLAGS 5 JSJZE A! 
WPS 


OISP DEFINITION 


0 gtds parameter 


FLAGS6 (SIZE 85 

} 

ICENT 0 NUMBER OF CENTRAL BODY (GTDS CONVENTION) 

ISUN 4 SUN-CENTRAL BODY FLAG 


FLAGS? (SIZE 10) 

NB0PT(3) 0 ARRAY OF GRAVITATING BODIES (GTDS CONVENTION) 

NB C NUMBER OF BODIES 


FLAGS9 (SIZE 4) 
NSTR 


MAX NUMBER OF ACCELERATION POINTS 


FSUNLT (SIZE 8) 

fsunlt 0 occultation flag 
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NAME (DIM) 

DISP 

DEFINITION 




GOaTAI (SIZE A«) 






XL (6! 

0 

STATE VECTOR OF 

S/C 

AT 

LIBRATION POINT 

DL 

30 

JULIAN DATE OF 

S/C 

AT 

LIBRATION POINT 

Xfi(6t 

38 

STATE vector of 

S/C 

AT 

BURN INITIATION 

TB 

68 

BURN time 




XF (ft)' ■ 

70 

state vector of 

S/C 

AT 

TCA 

TDUR 

AO 

flight time 





GDOiT^a (SIZE 2DO) 

BST'^(6®6) 0 
ASTM(6®3) 120 
ST'^<&«6) 180 


st&te partials 

BURN PARTIALS 

state Partials 


OVER burn phase 

OVER COAST PHASE 


GO®TA3 (SIZE 18) 


ACCTHOl 



BURN PARTIALS AT BURN EVENT INITIATION 


GDATAA (SIZE 18) 

V0LP(3) 0 INERTIAL VELOCITY OF LIBRATION POINT 


GFL'AG (SIZE 18) 


IBURN 

0 

FLAG TO 

INDICATE 

FINITE BURN 

INTEGRATION 

INCFLG 

4 

NOT USED 



ITOL 

e 

FLAG TO 

INDICATE 

TARGETING CONVERGENCE 

ITMAX 

c 

MAXIMUM*! 

NUMBER OF 

TARGETING ITERATION 

iter 

10 

CURRENT 

TARGETING 

ITERATION 


ibtype 

1 A 

FLAG TO 

INDICATE 

TYPE OF BURN 

STRATEGY 


GTAR (SIZE ^0) 


OTaR(35 0 TARGET VECTOR OF DESIRED VALUES 

DT0L(3) 18 TARGET VALUE TOLERANCE VECTOR 

PERT(2) 30 PERTURBATION VECTOR 
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NAME (DIMI 


DISP 


definition 


HLDATA (SIZE 

44) 



ZBIAS(6! 


0 

VECTOR added to LIBRATION STATE FOR TARGETING 

RLIBR(2> 


30 

EARTH STATE SCALING CONSTANTS REQUIRED TO 
generate state vector of LIBRATION POINT 

LIBR 


40 

LIBRATION POINT OF INTEREST 

lOFLAG (SIZE 

338) 



I DON 


0 

FLAG TO INDICATE ORBIT FILE IS BEING WRITTEN 
ON THIS PASS 

lOISK 


4 

FLAG TO INDICATE IF AN ORBIT FILE IS TO BE 
WRITTEN FOR THIS CASE 

NCPR 


0 

NOT USED 

NPOINT 


C 

number of SPECIAL PRINT POINTS 

TP 


10 

VECTOR OF SPECIAL PRINT POINTS 

TMPR 


330 

NOMINAL PRINT INTERVAL 

KNSTN2 (SIZE 

28) 



PI 


0 

PI 

PI2 


8 

TWO PI 

RPO 


1 0 

RADIANS PER DEGREE 

SPD 


18 

SECONDS PER DAY 

XKMPAU 


20 

kilometers PE ASTRONOMICAL UNIT 

KNSTN3 (SIZE 

108) 



SMU(ll) 


0 

VECTOR OF GRAVITATIONAL CONSTANTS (STEAP 
CONVENTION) 

RSOKIU 


58 

RADIUS OF SPHERE OF INFLUENCE (STEAP CONVENTION) 

RP ( 1 1) 


80 

RADIUS OF PLANETS (STEAP CONVENTION) 

LPROF (SIZE 8) 



LAUNCH 


0 

FLAG TO INDICATE TARGET-LAUNCH COMPATIBILITY 

LNON 


4 

FLAG TO INDICATE THAT LAUNCH PROFILE IS TO 
BE GENERATED ON THIS PASS 
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name (DIM) 


DISP 


DEFINITION 


LA6S10 JSIZE 8) 

NOCOWL 0 NUMBER OF TIMES CSTEP WAS CALLED 

ITERS A total NUMBER OF CORRECTOR ITERATIONS 


LAGSll (SIZE ICl 


N1 

0 

number 

OF 

TERMS 

IN 

INTEGRATION 

formulas 

FOR 



POSITION 

AND velocity 1 



N2 

A 

NUMBER 

OF 

terms 

IN 

INTEGRATION 

FORMULAS 

FOR 


PART I AS 

MAXIT(S) 8 MAXIMUM NUMBER OF ITERATIONS USED BY INTEGRATOR 


LAGS12 (SIZE A) 

K 0 current acceleration POINT 


LAGS13 (SIZE fl) 

N0FC(2) 0 number OF TIMES FROCES WERE REQUESTED 


FLAGSIA (SIZE fl) 

N0STEP(2) 0 GTDS COUNTER 

f 

LAGS15 (SIZE A) 

lELEVN 0 ACCUMULATED ACCELERATION POINTS TO BE WRITTEN 

TO THE ORBIT FILE 
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DEFINITION 


NAME(OIM) 

DE 

LOATA (SIZE 50 


FI 

0 

psn 

8 

PSI2 

10 

TIMl 

16 

TIM2 

20 

THELS 

28 

PHILS 

30 

THEDOT 

30 

RPRAT 

AO 

SIGNAL 

^8 

RETRO 

50 

KOAST 

58 


NOT USED 

FIRST INJECTION BURN ARC 

SECOND INJECTION BURN ARC 

duration of 1ST INJECTION BURN 

duration Of 2ND INJECTION BURN 

longitude OF LAUNCH SITE 

LATITUDE OF LAUNCH SITE 

POTATION RATE OF LAUNCH PLANET 

INVERSE OF PARKING ORBIT RATE 

OESIRED LAUNCH AZIMUTH 

FLAG TO INDICATE DIRECTION OF MOTION 

NOT USED 


LINKll (SIZE lOJ 

H 0 signed step SIZE (SEC) 

T a TIME UP TO WHICH INTEGRATION HAS PROGRESSED 


LINK12 (SIZE IB) 

SECTIM 0 VECTOR OF TIME WITH SECTION 


LINKU (SIZE 8) 

RM 0 RADIUS MAGNITUDE 


LINK19 (SIZE 50) 


TOL(IO) 0 VECTOR OF TOLERANCES USED BY INTEGRATOR 


LINKE3 (SIZE 160) 


ALPHA( 1 1 ) 

0 

PREDICTOR 

COEFFICIENT 

FOR 

POSITION 

TERMS 


ALPHAS ( 1 1 ) 

56 

CORRECTOR 

COEFFICIENT 

FOR 

POSITION 



ALPHAB( 1 1 ) 

BO 

PREDICTOR 

COEFFICIENT 

FOR 

POSITION 

partial 

TERMS 

BETAB( 11 > 

100 

PREDICTOR 

COEFFICIENT 

FOR 

VELOCITY 

partial 

TERMS 
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NAME (DIM) 


OISP 


DEFINITION 


LENKPR 

CETOL 0 TRUNCATION ERROR TOLERANCE 


LINK29 

X0(6«2) 0 array of STATE VECTORS OF NON-CENTRAL BODIES 

WRT central body 


LINK33 (SI^E 6E05 


Bcsniflio) 

ACS ( 1 1 « 10) 


0 starter corrector coefficients for POSITION 

370 STARTER CORRECTOR COEFFICIENTS FOR VELOCITY 


LINK35 (SIZE BO) 


SACS (U ) 
SBcsnu 


0 STARTER COEFFICIENTS FOR FIRST SUM 

58 starter COEFFICIENTS FOR SECOND SUM 


LINK36 (SIZE 210) 


X1D(1 1«3) 


0 array OF starter POSITION VECTORS 

8 ARRAY OF STARTEP VELOCITY VECTORS 


LINK37 (SIZE 60) 


ACCB (3) 

0 

ACCELERATION 

DUE 

TO 

CENTRAL BODY 

ACME (3) 

18 

acceleration 

DUE 

TO 

NON-CENTRaL BODIES 

ACTH (3) 

30 

acceleration 

DUE 

TO 

THRUST 

ANCF (3) 

■^8 

ACCELERATION 

DUE 

TO 

thrust and non-central BODIES 


LINK30 (SIZE 3F0) 


SRI (3) 

SX2 (3) 

SVl (3»20 ) 
SV2(35201 


0 1ST SUM FOR EQUATIONS OF MOTION 

18 2ND SUM FOR EQUATIONS OF MOTION 

30 1ST SUM FOR VARIATIONAL EQUATIONS 

210 2ND SUM rOR variational EQUATIONS 
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NAME (DIM) 


OISP 

DEFINITION 

LINK39 (SIZE 

10) 



YMDIC 

HMSIC 


O GO 

YEAR» MONTHt DAY OF INITIAL CONDITIONS 
HOUR* minute* SECONDS OF INITIAL CONDITIONS 

LINK40 (SIZE 

1F8) 



X0L0(3> 

YOLD(3.20) 


0 

18 

PREVIOUS POSITION AND VELOCITY 
PREVIOUS POSITION AND VELOCITY PARTIALS 

LINK41 (SIZE 

30) 



SPV(3) 

SPC<3) 


0 

IB 

INITIAL PREDICTED POSITION 
FINAL CORRECTED POSITION 

LINK42 (SIZE 

BO) 



BETAdl ) 
BETAS(H) 


O CO 

PREDICTOR COEFFICIENT FOR VELOCITY TERMS 
CORRECTOR COEFFICIENT FOR VELOCITY TERMS 

LINK43 (SIZE 

90) 



ACC8C<3»3) 

ANCBV(3i3) 



central BODY ACCELERATION FOR VARIATIONAL TERMS 
NON-CENTRAL BODY ACCELERATION FOR VARIATIONAL 
TERMS 



LINK44 (SIZE 108) 

RM2 0 RADIUS MAGNITUDE SQUARED 

RMS 6 RADIUS MAGNITUDE CUBED 

6MBM3(3) 10 GRAVITATIONAL CONSTANT 

BP(6«31) 28 ARRAY OF POSITIONS OF NON-CENTRAL BODIES WRT S/C 

BPM2(3) 08 square OF DISTANCE BETWEEN NON-CENTRAL BODIES 

AND S/C 

BPM(3) 00 DISTANCE BETWEEN NON-CENTRAL BODIES AND S/C 

6MQPM3(3) Ee GRAVITATIONAL CONSTANT 

GMRM3 100 gravitational CONSTANT 



2-10 



NAME (DIM) 


DlSP 


DEFINITION 


MATRIX (SIZE ICB) 



A (3®3) 

0 

matrix to convert SELENOCENTRIC TO SELENOGRAPHIC 

ADOT(3»3) 

46 

DERIVATIVE OF MATRIX * 

B(3o3) 

90 

MATRIX TO CONVERT EARTH INERTIAL TO EARTH BODY 

fixed 

C(3»3) 

08 

MATRIX TO CONVERT MEAN 1950,0 TO TRUE OF DATE 
COORDINATES 

6HA 

120 

GREENWICH HOUR ANGLE 

0X0(18) 

128 

NOT USED 

XP 

188 

X POLAR MOTION ANGLE 

YP 

ICO 

Y POLAR MOTION ANGLE 

MECEQ (SIZE 48) 



ECEO (393) 

0 

ECLIPTIC TO EARTH E(3UATORIAL ROTATION MATRIX 

MEQEC (SIZE 48) 



EQEC(393) 

0 

EARTH EQUATORIAL TO ECLIPTIC ROTATION MATRIX 

MLINK (SIZE 8) 



IPRE 

0 

FLAG TO INDICATE INITIAL RUN 

KWIT 

RLINK9 (SIZE 80) 

4 

FLAG TO INDICATE TERMINATION OF CURRENT CASE 

TWOPI 

0 

2.0 * PI 

GM(U ) 

8 

GRAVITATIONAL CONSTANTS OF CENTRAL BODIES 



(GTOS NUMBERING CONVENTION)^ 

MSGS (SIZE A) 


V 

HSGLVL 

0 

FLAG TO INDICATE PRINTING OF DEBUG DATA 

NEWLXi (SIZE 20) 



0PP(3) 

0 

VECTOR FROM SUN TO SPACECRAFT 

VBP 

18 

MAGNITUDE OF VECTOR 0PP 


BRAD 20 SOLAR RADIATION CONSTANT 


NAME (OIM> 


DISP 


DEFINITION 


NEWLK2 {SIZE 10) 

SCMREA 0 ■ AREA OF SPACECRAFT 

CSUBR 8 SPACECRAFT REFLECTIVITY CONSTANT 


PRT (SIZE 80) 

M0NTH(12) 
planet (11) 


0 HOLLERITH NAME OF MONTHS 

30 HOLLERITH NAME OF PLANETS 


RLINKA (SIZE 130) 


PVINT(6) 0 
AEINT(6! 30 
SPINT(fj) 60 
OBLINT(20) 90 


S/C INITIAL state VECTOR 
KEPLERIAN ELEMENTS 
spherical ELEMENTS 
AUXILLARY ELEMENTS 


RLINK5 (SIZE 5200) 


X(3J 

0 

S/C POSITION 

VECTOR 


XD(3) 

18 

S/C VELOCITY 

VECTOR 


XDD(40®3) 

30 

array of S/C 

ACCELERATION 

VECTORS 

XV(3»?0) 

3F0 

S/C POSITION 

PARTIALS 


XVD(3»20) 

500 

S/C velocity 

PARTIALS 


XVDO (40»3f 20) 

7B0 

array of S/C 

ACCELERATION 

PARTIALS 


RLINK6 (SIZE 270) 


PX(3»3) 
PXD(3t3) 
ACCPAR (3?20) 


0 acceleration PARTIALS WRT POSITION 

A8 ACCELERATION PARTIALS WRT VELOCITY 

90 ACCELERATION PARTIALS WRT PARAMETERS 


SLPOPT (SIZE 38) 
OJ. 

0 

EVAL internal 

IDAYl 

B 

DAY OF FIRST RECORD ON EPHEMERIS FILE 

lYEAR 

C 

EVAL internal 

I SPAN 

10 

eval internal 

N8EPM(3) 

14 

BODIES FOR polynomial COEFFICIENTS 

NDEGRE O) 

20 

DEGREE OF POLYNOMIALS 

NCFDAY 

2C 

number of days PER CURVE FIT 

ISLP50 

30 

EVAL internal 

NBSLP 

34 

EVAL internal 



DISP 


DEFINITION 


SLPRFC (SIZE D7C» 


TSEC ; 

PPOLY«392092) 

VPOLY (392092} 

APOL Y ( 3 9 3 9 1 0 } 
CPOLY (393o10} 
PDELH ( 10 > 

IDAY 


0 TIME IN SECONDS FROM START OF THIS YEAR TO 
MIDPOINT OF THIS RECORD TIME INTERVAL 
8 POLYNOMIAL COEFFICIENTS FOR THE POSITION 
coordinates of the two BODIES 
3CB polynomial COEFFICIENTS FOR THE VELOCITY 
COORDINATES OF THE Tao BODIES 
700 POLYNOMIAL COEFFICIENTS FOR THE MATRIX 

A58 POLYNOMIAL COEFFICIENTS FOR THE 'C« MATRIX 

D28 polynomial COEFFICIENTS FOR DELTA H 

070 beginning OAY of THIS RECORD 


THRUST (SIZE 60> 


THRMA6 

0 

THRUST MAGNITUDE 

XISP 

8 

SPECIFIC IMPULSE 

SCMASS 

10 

S/C MASS 

ALPHA 

la 

right ascension of burn vector 

BETA 

20 

declination of burn vector 

TBURN 

28 

BURN DURATION 

OMASS 

30 

MASS RATE 

CURMAS 

30 

S/C MASS AT CURRENT INTEGRATION TIME 

,COSA 

AO 

COSINE OF 5alPHA» 

COSR 

A0 

COSINE OF »0ETA8 

SINA 

50 

SINE OF ^ALPHAS 

SINB 

50 

SENF OF obeTA® 

TIMCOF (SIZE 

F10J 


COEF (72«6) 

0 

TIME DIFFERENCE COEFFICIENTS 

POCOF (A« 90 } 

6C0 

POLAR motion COEFFICIENTS 

JARG(72) 

CAO 

JULIAN DATES FOR TIME DIFFERENCE COEFFICIENTS 

lARGf-^^) 

D60 

JULIAN date INTERVALS FOR POLAR MOTION 



COEFFICIENTS 

NDAYS 

ElO 

SIZE OF VECTOR JARG 

NNOAYS 

E14 

SIZE OF VECTOR lARG 


USUN (SIZE 185 

USUN 0 UNIT VECTOR FROM CENTRAL BODY TO- SUN 
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NAME (DIM) 


DISP 


DEFINITION 


ZOATA 

ATRY 

NFR. 


(SIZE CP 

0 INITIAL GUESS FOR SOLUTION TO LAMBERTS PROBLEM 
8 number of FULL REVOLUTION ABOUT THE EARTH 
PRIOR TO encounter 




2-lA 


2.3„2 COMMON variables IN ALPHABETICAL ORDER 


variable (DIM) 

BLOCK 

Al3s>3> 

MATRIX 

ACCB131 

LINK37 

ACNB «3) 

LIWK37 

ACCBC<3«3» 

LINKA3 

ACCPAR(3«20) 

RLINK6 

ACCTH(3) 

GD A T A 3 

ACSai »10> 

LINK33 

ACTH(3> 

LINK37 

AOOTOoS) 

matrix 

AEINT(6J 

rlink« 

ALPHA 

THRUST 

ALPHA (U) 

LINK23 

ALPHAS U 1) 

LINK23 

ALPHABUl) 

LINK23 

ANCBV ( 3p3) 

LINKA3 

ANCF(3J 

link 37' 

APOLY I3e39 10) 

slprec 

ASTM<6«3) 

GDATA2 

ATRY 

ZDATA 

B(3»3) 

matrix 

BCS(1J»10) 

LINK33 

BETA 

Thrust 

BETA n 1 ) 

LINK«2 

BETABIU) 

LINK23 

BETASai) 

LINKA2 

0P(6*3) 

LINKA<> 

8PM(3) 

LINK^<v 

BPM2 J3J 

link^-^ 

BPP(35 

NEWLKl 

BRAD 

NEWLKl 

BSTM(!696) 

GDATA2 

COsS) 

matrix 

CETOL 

LINK28 

C0EF<7296) 

TIMCOF 

COSA 

THRUST 

COSB 

THRUST 

CPOLY!3t39lO} 

SLPREC 

CSUBR 

NEWLK2 

curmas 

THRUST 

DJ 

■SLPOPT 

OL 

GD A T A 1 

DM ASS 

THRUST 

DTaR(3> 

GTAR 

OTOL<3) 

gtar 

OXQ < 18) 

matrix 

ECEQ (3«3J 

MECEO 

EQEC(3c3j 

MEQEC 


DEFINITION 

matrix to convert selenocentric to selenographic 

ACCELERATION DUE TO CENTRAL BODY 

ACCELERATION DUE TO NON-CEnTRAL BODIES 

central body ACCELERATION FOR VARIATIONAL TERMS 

ACCELERATION PARTIALS t«pT PARAMETERS 

BURN PARTIALS AT BURN EVENT INITIATION 

starter corrector COEFFICIENTS FOR VELOCITY 

ACCELERATION DUE TO THRUST 

derivative OF MATRIX “A® 

keplerian elements 

right ascension of burn vector 

PREDICTOR COEFFICIENT FOR POSITION TERMS 

CORRECTOR COEFFICIENT FOR POSITION 

predictor COEFFICIENT FOR POSITION PARTIAL TERMS 

non-central body acceleration for variational 

TERMS 

ACCELERATION DUE TO THRUST AND NON»CENTRAL BODIES 
polynomial COEFFICIENTS FOR THE ® A MATRIX 
BURN PARTIALS 

INITIAL GUESS FOR SOLUTION TO LAMBERTS PROBLEM 
MATRIX TO CONVERT EARTH INERTIAL TO EARTH BODY 
FIXED 

starter CORRECTOR COEFFICIENTS FOR POSITION 
DECLINATION OF BURN VECTOR 
PREDICTOR COEFFICIENT FOR VELOCITY TERMS 

PREDICTOR COEFFICIENT FOR VELOCITY PARTIAL TERMS 

CORRECTOR COEFFICIENT FOR VELOCITY TERMS 

ARRAY OF POSITIONS OF NON-CENTRAL BODIES WRT S/C 

DISTANCE between NON-CENTR AL BODIES AND S/C 
square of DISTANCE BETWEEN NON-CENTRAL BODIES 
AND S/C 

VECTOR FROM SUN TO SPACECRAFT 

SOLAR RADIATION CONSTANT 

state PARTIALS OVER BURN PHASE 

MATRIX TO CONVERT MEAN 1950o0 TO TRUE OF' DATE 

COORDINATES 

truncation error tolerance 

time DIFFERENCE COEFFICIENTS 
COSINE OF «ALPHA« 

COSINE OF «BETA» 

POLYNOMIAL COEFFICIENTS FOR THE “C® MATRIX 
SPACFCPAFT REFLECTIVITY CONSTANT 
S/C MASS AT CURRENT INTEGRATION TIME 
FVAL internal 

JULIAN DATE OF S/C AT LIBRATION POINT 
MASS RATE 

target VECTOR OF DESIRED VALUES 
target VALUE TOLERANCE VECTOR 
NOT USED 

ECLIPTIC TO EARTH EQUATORIAL ROTATION MATRIX 
EARTH EQUATORIAL TO ECLIPTIC ROTATION MATRIX 


2-15 



VARIABLE JOIMJ 

BLOCK 

DEFINITION 

FI 

LDATA 

NOT USED 

FSUWLT 

fsunlt 

occultation flag 

GHA 

matrix 

GREENWICH HOUR ANGLE 

GM< U ) 

RLINK9 

gravitational CONSTANTS OF CENTRAL BODIES 
(Gins NUMBERING CONVENTION) 

GMBM313S 

LINX44 

gravitational CONSTANT 

GMBPM3<3) 

LIWK44 

gravitational constant 

GMRM3 

LINK44 

gravitational constant 

H 

LlNKll 

SIGNED step size 5SECI 

HHSIC 

L!NS<39 

HOUR» minute* seconds of initial conditions 

IARG(44; 

TIMCOF 

JULIAN DATE INTERVALS FOR POLAR MOTION 
COEFFICIENTS 

IBODY 

FLA6S1 

CENTRAL BODY NUMBER <GTDS CONVENTION) 

IRTVPE 

GFL AG 

FLAG TO INDICATE TYPE OF BURN STRATEGY 

IBURN 

GFL AG 

FLAG TO INDICATE FINITE BURN INTEGRATION 

ICENT 

FLAGS6 

number of CENTRAL BODY JGTDS CONVENTION) 

ICENTB 

FLAGS2 

CENTRAL BODY NUMBER (GTOS CONVENTION) 

IDAY 

SLPREC 

BEGINNING DAY OF THIS RECORD 

IDAYl 

SLPOPT 

DAY OF FIRST RECORD ON EPHEMERIS FILE 

IDISK 

lOFLAG 

FLAG TO INDICATE IF AN ORBIT FILE IS TO BE 
WRITTEN FOR THIS CASE 

IDON 

I OFL AG 

FLAG TO INDICATE ORBIT FILE IS BEING WRITTEN 
ON THIS PASS 

lELFVN 

LAGS15 

ACCUMULATED ACCELERATION POINTS TO BE WRITTEN 
TO THE ORBIT FILE 

incflg 

-GFLAG 

NOT USED 

IND(25> 

FLA6S2 

FLAGS FOR CURRENT SECTION 

IN0SEC< 10«3) 

FLAGS2 

FLAGS FOR SECTION 

IN0X4 

FLA6S1 

not used 

INDYfc 

FLAGSl 

not used 

3SLP50 

SLPOPT 

EVAL internal 

ISPAW 

SLPOPT 

EVAL INTERNAL 

ISUN 

FLAGS6 

SUN-=CENTRAL BODY FLAG 

ITER 

gflag 

CURRENT TARGETING ITERATION 

ITERS 

LAGSlO 

total number of CORRECTOR ITERATIONS 

UMAX 

GFL AG 

MAXIMUM number OF TARGETING ITERATION 

ITOL 

GFLAG 

FLAG TO INDICATE TARGETING CONVERGENCE 

lYEAR 

SLPOPT 

EVAL INTERNAL 

JARGU2J 

TIMCOF 

JULIAN DATES FOR TIME DIFFERENCE COEFFICIENTS 

K 

LA6S12 

CURRENT ACCELERATION POINT 

KOAST 

LDATA 

NOT USED 

KSTATE ibl 

FLAGS3 

VECTOR OF LABEL NUMBERS FOR STATE UNKNOWNS 

KW3T 

ML INK 

FLAG TO INDICATE TERMINATION OF CURRENT CASE 

LAUNCH 

LPROF 

FLAG TO INDICATE TARGET-LAUNCH COMPATIBILITY 

LE0R 

MLDATA 

LI8RATI0N POINT OF INTEREST 

LNON 

LPROF 

FLAG TO INDICATE THAT LAUNCH PROFILE IS TO 
BE GENERATED ON THIS PASS 

MANDON 

FLAGS3 

maneuver OCCURENCE FLAG 

MAXIT(S) 

LA6SH 

MAXIMUM number OF ITERATIONS USED BY INTEGRATOR 

month U?J 

PRY 

HOLLERITH NAME OF MONTHS 

MSGLVL 

MSGS 

FLAG TO INDICATE PRINTING OF DEBUG DATA 
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VARIABLE SDIMJ) 
NB 

NB(3) 

NBEPMO) 

NBSLP 

NROD 

NBOPTO) 

NCFOAY 

NCNM 

ncondt 

ncpr 

NCURS 

MOAYS 

NDE6RE (3> 

NENTRY 

NEQ 

NFR 

N6 

NNDAYS 

NOCOWL 

NOFC {21 

N0STEP(2> 

NPOINT 

NPS 

NSECTW 

NSNM 

nstate 

NSTR 

NTSFQS 

Nr 

N2 

OBLINT{20> 

POELHUOii 

PERT{2) 

PHILS 

PI 

PI2 

planet (1 1 J 
POCOF 

PPOLY(3o20o2) 

PSIl 
PSI2 
PVINT (6) 

PK 53o3! 

PXD (3®3li 


BLOCK^ DEFINITION 

FLAGS? number OF BODIES 

BDATA VECTOR OF BODIES {STEAP CONVENTION) 

SLPOPT bodies for POLYNOMIAL COEFFICIENTS 

slpopt eval internal 

0DA7A number of GRAVITATING BODIES 

FLAGS? 'ARRAY OF GRAVITATING BODIES «GTDS CONVENTION) 

SLPOPT number of DAYS PER CURVE FIT 

FLAGS3 number OF C{N«M) TO BE ESTIMATED 

FLAGS3 condition NUMBER 

lOFLAG NOT USED 

FLAGSA 6TDS PARAMETER 

TIMCOF SIZE OF VECTOR JARG 

SLPOPT DEGREE OF POLYNOMIALS 

FLAGSA 6TDS PARAMETER 

FLAGSl number OF VARIATIONAL EQUATIONS 
ZDATA NUMBER OF FULL REVOLUTION ABOUT THE EARTH 
PRIOR TO ENCOUNTER 

FLAGS3 total NUM9RE OF POTENTIAL COEFFICIENTS TO BE 
ESTIMATED 

IIMCOF SIZE OF VECTOR EARG 
LAGSIO NUMBER OF TIMES CSTEP WAS CALLED 
LAGS13 NUMBER OF TIMES FROCES WERE REQUESTED 
LAGS14 GTDS COUNTER 

IOFLAG NUMBER OF SPECIAL PRINT POINTS 

FLAGSS GTDS PARAMETER 

FLAGS2 CURRENT SECTION NUMBER 

FLAGS3 number OF S(N«M) TO BE ESTIMATED 

FLAGS3 number OF STATE UNKNOWNS 

FLAGS9 MAX NUMBER OF ACCELERATION POINTS 

FLAGS2 total NUMBER OF SECTIONS 

LAGSU NUMBER OF TERMS IN INTEGRATION FORMULAS FOR 
POSITION AND velocity 

LAGSU NUMBER OF TFRMS IN INTEGRATION FORMULAS FOR 
PARTIAS 

RLINKA AUXILLARY ELEMENTS 

SLPREC POLYNOMIAL COEFFICIENTS FOR DELTA H 
GTAR PERTURBATION VECTOR 
LDATA LATITUDE OF LAUNCH SITE 
KNSTN2 PI 
KNSTW2 TWO PI 

PRT HOLLERITH NAME OF PLANETS 

TIMCOF polar MOTION COEFFICIENTS 
SLPREC polynomial COEFFICIENTS FOR THE POSITION 
COORDINATES OF THE TWO BODIES 
LDATA FIRST INJECTION BURN ARC 
LDATA SECOND INJECTION BURN ARC 
RLINK4 S/C INITIAL STATE VECTOR 
RLINK6 ACCELERATION PAR7IALS WRT POSITION 
RLINK6 ACCELERATION PARTIALS WRT VELOCITY 


Z) 
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VARIABLE(DIM) 

BLOCK 

RETRO 

LOATA 

RLI9R(2) 

HLDATA 

RM 

LINK14 

RM2 

LINKA4 

RM3 


RPdi) 

KNSTN2 

RPD 

KNSTN2 

RPRAT 

LOATA 

RSOI (11) 

KNSTN2 

SACS( 11 > 

LINK35 

SBCS( U ) 

LINK35 

SCAREA 

NEWLK2 

SCMASS 

THRUST 

SECT I M 

LINK12 

SIGNAL 

LDATA 

SINA 

thrust 

SINB 

thrust 

SMU(ll) 

KNSTN3 

SPC<3) 

LINKAl 

SPD 

KNSTN2 

SPINT(6) 

RLINKA 

SPV(3) 

LINK41 

STM(6s6) 

GDATA2 

5V1 (3f20> 

LINK38 

SV2(3»20) 

LINK38 

SXU3) 

LINK38 

SX2(3) 

LINK38 

T 

LINKll 

TB 

GDAT A 1 

TBURN 

THRUST 

TOUR 

GDAT A 1 

thedot 

LOATA 

thels 

LDATA 

thrmag 

thrust 

TIMl 

LDATA 

TIM2 

LOATA 

TMPR 

lOFLAA 

tsec 

SLPREC 

TOL (10 ) 

LINK19 

TP 

lOFLAG 

TWOPI 

RLINK9 


DEFINITION 

FLAG TO INDICATE DIRECTION OF MOTION 

EARTH state SCALING CONSTANTS REQUIRED TO 

GENERATE STATE VECTOR OF L10RATION POINT 

RADIUS MAGNITUDE 

RADIUS MAGNITUDE SQUARED 

RADIUS magnitude CUBED 

RADIUS OF planets 5STEAP CONVENTION) 

RADIANS PER DEGREE 

INVERSE OF PARKING ORBIT RATE 

RADIUS OF SPHERE OF INFLUENCE (STEAP CONVENTION) 

starter COEFFICIENTS FOR FIRST SUM 

starter COEFFICIENTS FOR SECOND SUM 

AREA OF SPACECRAFT 

S/C MASS 

VECTOR OF TIME WITH SECTION 
DESIRED LAUNCH AZIMUTH 
SINE OF »ALPHA« 

SINE OF *BETa« 

VECTOR OF GRAVITATIONAL CONSTANTS (STEAP 
CONVENTION) 

final CORRECTED POSITION 

SECONDS PER DAY 

SPHERICAL ELEMENTS 

INITIAL PREDICTED POSITION 

STATE PARTIALS OVER COAST PHASE 

1ST SUM for variational EQUATIONS 

2ND SUM FOR variational EQUATIONS 

1ST SUM FOR EQUATIONS OF MOTION 

2ND SUM FOR EQUATIONS OF MOTION 

TIME UP TO WHICH INTEGRATION HAS PROGRESSED 

BURN time 

BURN DURATION 

FLIGHT TIME 

potation RATE OF LAUNCH PLANET 
longitude of LAUNCH SITE 
THRUST magnitude 
duration of 1ST INJECTION BURN 
DURATION OF 2ND INJECTION BURN 
NOMINAL PRINT INTERVAL 

TIME IN SECONDS FROM START OF THIS YEAR TO 
MIDPOINT OF THIS RECORD TIME INTERVAL 
VECTOR OF TOLERANCES USED BY INTEGRATOR 
VECTOR OF SPECIAL PRINT POINTS 
2.0 * PI 
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VARIABLE (DIM) 


BLOCK 


DEFINITION 


USUN 

VRP 

VOLPO) 

VPOL Y ( 3 5 20 9 2 ) 

X(3) 

XBC6) 

XB(6«2) 


XDO) 

XDOCA093) 

SF (6) 

XISP 
XKMPAU • 

XL (6) 

XP 

X0L0(3) 

XV(3920! 

XVDOoZO) 

XVDD (AO 9392 O ) 
XI (1193) 
XlD(n93S 
YMDIC 
YOLD(3s20) 



USUN UNIT VECTOR FROM CENTRAL BODY TO SuN 
NEWLKl MAGNITUDE OF VECTOR 3PP 
60ATAA inertial VELOCITY OF LIBRATIOn POINT 
SLPREC polynomial COEFFICIENTS FOR ThE VELOCITY 
COORDINATES OF ThE TWO BODIES 
RLINK5 S/C POSITION VECTOR 

6DATA1 STATE VECTOR OF S/C AT BURN INITIATION 
LINKER ARRAY OF STATE VECTORS OF NON-CENTrAL BODIES 
t^RT CENTRAL BODY 
RLENKS S/C velocity VECTOR 
RLINK5 array OF S/C ACCELERATION VECTORS 
GDATAl .state VECTOR OF S/C AT TCA 
THRUST SPECIFIC IMPULSE 
KNSTN2 KILOMETERS PE ASTRONOMICAL UNIT 
GDATAl state VECTOR OF S/C AT LIBRATION POINT 
MATRIX X POLAR MOTION ANGLE 
LINK40 PREVIOUS POSITION AND VELOCITY 
RLINK5 S/C POSITION PARTIALS 
RLINK5 S/C velocity PARTIALS 
RLINKS array OF S/C ACCELERATION PARTIALS 
LINK36 array of STARTER POSITION VECTORS 
LINK36 ARRAY OF STARTER VELOCITY VECTORS 
LINK39 YEAR? MONTH? DAY OF INITIAL CONDITIONS 
LINKAO PREVIOUS POSITION AND VELOCITY PARTIALS 
MATRIX Y POLAR MOTION ANGLE 

HLDaTA VECTOR ADDED TO LIBRATION STATE FOR TARGETING 
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ESRAN RH0GR6M STRUCTURE! 


3ol EHRAN program DESCRIPTION 

The STEAP-L/ERRAN program is a digital computer program (WRITTEN 

IN fortran IV COMPATI0LE WITH THE IRM 360/370 SYSTEM® THE PROGRAM 
CONTAINS ABOUT SlATV SUBROUTINES OF WHICH ABOUT A TENTH ARE ASSOCIATED 
WITH THE COWELL FILE HEADER* ABOUT A SIXTH ARE GENERAL UTILITY ROUTINES* 
ABOUT A THIRD ARE COMMON TO MEASUREMENT AND E^ENT PROCESSING* AND THE 
REST ARE SPECIFIC TO MEASUREMENT PROCESSING OR TO PARTICULAR EVENT 
'PROCESSING* The program ACCESSES FIVE DATA SETS DURING OPERATIONS 
THE standard fortran INPUT tSYSiN)* THE PRESENTATION OUTPUT (SYSOUT=A) 
AND PUNCH OUTPUT (5YS0UT=8) DATA SETS* THE DIRECT ACCESS SOLAR/LUNAR/ 
planetary EPhEMERIS file used by the GTDS JDSn=GTDS«SLP19S0® JAN71 ) » 

AND THE sequential ORBIT FILE ON WHICH THE TRAJECTORY AND STATE TRANSI- 
TION MATRIX data are STORED. THE {OVERLAID) PROGRAM REQUIRES ABOUT 
310K BYTES STORAGE. 


3.2 ERRAN SUBROUTINE HIERARCHY 

FIGURE 3.1 illustrates THE GENERAL PROGRAM STRUCTURE® THE MAIN 
PROGRAM CALLS TWO EXECUTIVE ROUTINES: DATA* FOR DATA INPUT AND INI- 

TIALIZATION* and ERRaN? for MEASUREMENT AND EVENT PROCESSING®. ERRAN 
CONTROLS measurement PROCESSING DIRECTLY* BUT CALLS SPECIFIC (0\/EHLAID) 
subroutines TO PROCESS PARTICULAR EVENTS® 'ERRAN CALLS SETEVN FOR THE 
EIGENVECTOR EVENT* WHOSE PROCESSING IS COMMON TO ALL EVENTS® ERRAN 
THEN WILL CALL PRED IF PROCESSING A PREDICTION EVENT* OR WELL CALL GUIDM 
IF processing a guidance event OH A FINAL INSERTION EVENT® ERRAN ThEN 
MAY CALL GENGID IF PROCESSING A GUIDANCE EVENT WITH GENERALIZED 
covariance^® chapter 4 provides an index OF ALL SUBROUTINES USED IN 
ERRAN AND GIVES INDIVIDUAL ROUTINE DOCUMENTATION INCLUDING DESCRIPTIONS* 
ANALYSES* and FLOWCHARTS® 


3,3 ERRaN COMMON VARIABLES 


The following two subsections provide definitions of 

APPEARING IN COMMON BLOCKS IN THE ERRAN PROGRAMS® SUBSECT 
LISTS the common BLOCKS IN ALPHABETICAL ORDER* GIVES THE S 


BLOCKS* AND DEFINES 
3®3®2 LISTS all THE 
BLOCK TO WHICH EACH 


THE VARIABLES APPEARING IN EACH BLOCK® 
common variables in ALPHABETICAL ORDER 
BELONGS* AND DEFINES THE VARIABLES® 


THE VARIABLES 
ION 3®3.l 
IZE OF THE 
SUBSECTION 
* GIVES THE 
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Figure 3. 1 ERRAN SUBROUTINE HIERARCHY 
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Figure 3. 1 ERRAN SUBROUTINE HIERARCHY (Cont) 
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Figure 3,1 ERRAN SUBROUTINE HIERARCHY (Cent) 
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Figure 3, 1 ERKAN SUBROUTINE HIERARCHY (Cont) 







3c3«,l common slocks IN ALPHABETICAL ORDER 

PREVIOUS versions OF STEaP REQUIRED THE FOLLO»<ING COMMON BLOCKS 
WHICH ARE NO LONGER NECESSARY: C0NST3(SIZE 40), ENCKEISiZE 4), 
PRELCM(S17E 3CA), STMG(SIZE 4), TRAJCD(SIZE 301, TRJ(SIZE IBO), 
and update (SIZE 54) e 


NAME (DIM) 

DISP 

DEFINITION 

BLK (SIZE ‘J54) 



T 

0 

trajectory TIME IN DAYS 

PMASS( 1 1 ) 

8 

GRAVITATIONAL CONSTANTS ( A . U , **3/DAY**2 ) 

COMEG (4,9) 

60 

* NO LONGER USED 

C INC (4,9) 

180 

* NO LONGER USED 

C0MEGT(4,9) 

2A0 

*** NO LONGER USED 

SMJR (2,9) 

3C0 

* NO LONGER USED 

CECC(4»9) 

450 

* NO LONGER USEO 

CMEAN(4,9) 

570 

* NO LONGER USED 

MUPLAN(ll), 

690 

GRAVITATIONAL CONSTANTS ( KM ( ( 3/SEC»*2 ) 

CSAX(2,9) 

6Efl 

NO LONGER USED 

EMN( 15) 

770 

NO l,.ONGER USED 

RADIUS ( n ) 

7F0 

RADII OF THE PLANETS (A.U.) 

RMASS (11) 

848 

GRAVITATIONAL CONSTANTS (RELATIVE TO SUN) 

SPHERE (11) 

SAD 

SPHERES OF INFLUENCE (A.U.) 

XP(S) 

8F8 

STATE vector OF PLANET 

NO( n ) 

928 

NO LONGER USED 


COM (SIZE 4,0) 



PI 

0 

MATHEMATICAL CONSTANT 

RAD 

8 

number OF degrees per 

ITRAT 

10 

* NO LONGER USED 

KOUNT 

14 

* NO LONGER USEO 

INCMNT 

18 

* NO LONGER USEO 

INCPR 

IC 

* NO LONGER USEO 

INC 

20 

*«» NO LONGER USED 

IPR 

24 

* NO LONGER USED 

N80DYI 

20 - 

* NO LONGER USED 

NBODV 

2C 

,« NO LONGER USED 

IPRT (4) 

30 

* NO LONGER USEO 
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NAMEIDIM) niSP DEFINITION 


CONST (SIZE D^) 




OMFSA 

0 

POTATION PATE OF EARTH 


EPS 

8 

EARTH OBLIQUITY 


SAL(3) 

10 

STATION altitudes (ABOVE 

RADIUS OF EARTH) 

SLAT O) 

28 

STATION latitudes 


SLONO) 

^0 

station longitudes ' 


ONCN (3) 

58 

dynamic NOISE constants 


HNCN(12) 

TO 

MEASUREMENT NOISE CONSTANTS 

NST 

DO 

NUMBER OF stations TO BE 

USED (MAXIMUM 3) 


CONST2 (SIZE 58) 


UST(3) 

0 

X-DIRFCTIOM 

COSINE FOP 

STARS 


VST (3) 

18 

Y-DIRECTION 

COSINE FOP 

STARS 


i«iST(3) 

30 

z-otrection 

COSINE for 

stars 


FOP 

48 

off-diagonal 

annihilat I 

ON VALUE 

(POSITION) 

FOV 

50 

off-diagonal 

ANNIHILATION VALUE 

(VELOCITY) 



DPNU^) (SIZE 88} 


ZERO 

0 

dourle-precision 

VALUE 

OF 

ZERO (0»0) 

ONE 

8 

DOUBLE-PRECISION 

VALUE 

OF 

one (1»0) 

TsiO 

10 

DOUBLE-PRECISION 

VALUE 

OF 

TWO (2.0) 

HALF 

18 

DOUBLE-PRECISION 

VALUE 

OF 

HALF (0.5) 

THREE 

20 

double-precision 

VALUE 

OF 

Three (3.0) 

EMI 

28 

DOUBLE-PRECISION 

VALUE 

OF 

1 .E-l 

EM2 

30 

DOUBLE-PRECISION 

VALUE. 

OF 

1 .E-2 

EM3 

38 

DOUBLE-PRECISION 

VALUE 

OF 

1 ,E-3 

EM4 

40 

DOUBLE-PRECISION 

VALUE 

OF 

1 ,E-4 

EM5 

48 

double-precision 

VALUE 

OF 

1 =E-5 

EM6 

50 

DOUBLE-PRECISION 

VALUE 

OF 

1 .£-6 

EM7 

58 

DOUBLE-PRECISION 

VALUE 

OF 

l.E-7 

EM8 

&0 

DOUBLE-PRECISION 

VALUE 

OF 

1 .E-8 

EM9 

68 

DOUBLE-PRECISION 

VALUE 

OF 

1 ,E-9 , 

EP50 

70 

DOUBLE-PRECISION 

VALUE 

OF 

1 ,E*50 

TisfOPI 

T8 , 

double-precision 

VALUE 

OF 

2.<»PT 

EM 13 

80 

double-precision 

VALUE 

OF 

1 .E-13 - 
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1 


NAME (DIM) 

DISP 

DEFINITION 


EVENT (SIZE 4R4) 




TEV (50) 

0 

scheduled TIMES OF EVENTS 


TPT2 (20) 

190 

TIMES PREDICTED TO IN PREDICTION 

EVENTS 

SI6RES 

230 

VARIANCE OF RESOLUTION ERROR 


SIGPRO 

238 

VARIANCE OF PROPORTIONALITY ERROR 

SIGALP 

240 

VARIANCE OF ERROR IN POINTING ANGLE 1 

SIGBET 

248 

variance of error in POINTING ANGLE 2 

HP7 

250 

« NO LONGER USED 


P7 

258 

* NO LONGER USED 


TAU7 

260 

* NO LONGER USED 


AINC7 

268 

NO LONGER USED 


ANODE7 

270 

* NO LONGER USED 


PERP7 

270 

* NO LONGER USED 


ECC7 

260 

* WO LONGER USED 


DV8 (35 

288 

* NO LONGER USED 


BRNTIM 

2A0 

duration (DAYS) OF FINAL INSERTION BURN 

NEV 

2A8 

number of EVENTS SCHEDULED 


lEVNT(SO) 

2AC 

COOED event TYPES CORRESPONDING 

TO REV 

IHYPl 

374 

NO LONGER USED , 


lEIG 

378 

NO LONGER USED 


ICDT3(20> 

37C 

CODES FOR guidance POLICIES 


NPE 

3CC 

COUNT OF prediction EVENTS 


N6E 

3D0 

count of guidance events 


IPOL 

3D4 

NO longer USFO 


I IPOL 

3D0 

NO LONGER USED 


IC003 (20) 

30C 

NO longer used 


NEVl 

42C 

NUMBER OF SCHEDULED EIGENVECTOR 

EVENTS 

NEV2 

430 

NUMBER OF SCHEDULED PREDICTION 

EVENTS 

NEV3 

434 

NUMBER OF SCHEDULED GUIDANCE 

EVENTS 

NEVA 

438 

NUMBER OF SCHEDULED INSERTION 

EVENTS 

NQE 

436 

NO LONGER USED 


NEV5 

440 

NO LONGER USED 


NEV6 


NO LONGER USED 


NAE 

44B 

NO LONGER USED 


NAF6 (20) 

44C 

NO LONGER USED 


NEV7 

49C 

NO LONGER USED 


I0P7 

4A0 

NO LONGER USED 


NEV8 

4 A4 

NO LONGER USED 


NEV9 

4A8 

NO longer USED 


NEVIO 

4 AC 

NO LONGER USED 


NEVH 

400 

NO LONGER USED 



TI^'ES 


FLfi6S2 (SIZE E8) 


ICENTB 
NTSEQS 
NSECTW 
INOSEC (10*3) 
IND (25) 


0 CENTRAL BODY NUMBER 

total number of sections 

8 CURRENT SECTION NUMBER 

C SECTION FLAGS 

84 CURRENT SECTION FLAGS 
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DI5P 


definition 


NAME (DIM) 


GAINC (SI7E 1628) 

PMIN(6f6) 0 POS/VEL covariance BEFORE MEASUREMENT (WLS) 

PSMINdSslS) 120 SOLVE-FOR COVARIANCE BEFORE MEASUREMENT (WLS) 

CMIN(6d5) 028 CORRELATION MATRIX, POS/VEL AND SOLVE-FORS 

PPLU(6,6) AF8 POS/VEL COVARIANCE AFTER MEASUREMENT (WLS) 

PSPLU(15*15) Cle SOLVE-FOR COVARIANCE AFTER MEASUREMENT (WLS) 

CPIU(6,15) 1320 CORRELSTTON MATRIX, POS/VEL AND SOLVE-FORS 

RSAVE(6) 15F0 state VECTOR AT TLAST 

TLAST 1620 TIME WHEN MEASUREMENT LAST PKOCESSEO 


GCA (SIZE EO) 

XIGdS) 0 IGNORE PARAMETER LABLES 

IAUGW(24) 70 IGNORE PARAMETER AUGMENTATION VECTOR 

NDIMA or dimension of IGNORE PARAMETER STATE 

I6EN DC =0, PERFORM NO GENERALIZED COVARIANCE ANALYSIS 

=1, PERFORM GENERALIZED COVARIANCE ANALYSIS 



GENGD (SIZE 40) 


EE (4) 
EEE (4) 

0 

20 

ACTUAL MEANS OF EXECUTION ERROR PARAMETERS 
VARIANCES OF EXECUTION ERROR PARAMETERS 

GENGDl (SIZE 

23E8) 




ACTUAL CONTROL SECOND MOMENT MATRICES 

GPG(6»6) 

0 

state 

GCXXSG(6,15) 

120 

STATE/SOLVE-FOR VECTOR 

6CXUG (650) 

3F0 

NO LONGER USED 

GCXVG(65l5)' 

570 

STATE/MEASUREMENT CONSIDERS 

GCXWG (6dS) 

840 

STATE/IGNORE PARAMETERS 

GPSGdbdS) 

BIO 

SOLVE-FOR VECTOR 

6CXSUG ( 15«8) 

1210 

NO LONGER USED 

GCXSV6 ( 15o 15) 

1508 

SOLVE-EOR/MEASUREMENT CONSIDERS 

GCXSWG (15,15) 

ICEO 

solve-for/ignore parameters 
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GENRL <SIZE 57081 


GP 1 6 0 6 1 

0 

GCXXS«6*15> 

'120 

GCXU(6*8> 

3F0 

GCXV (6* 15) 

570 

GPSn5»15) 

840 

GCXSU(15«8) 

F4B 

GCXSVU5* IS) 

1308 

GCXS«M15*1S) 

lAlO 

JPR (A*4) 

2118 

TXW (6* 15) 

2198 

AN (A* 15) 

2468 

GCUV(8*15) 

2648 

GCUW(8* 15) 

2A08 

6CVW(15,15) 

2DC8 

GU(8«8) 

34D0 

6V(15,15) 

36D0 

GW ( 15*15) 

3DD8 

GDNCN (3) 

44E0 

GMNCN (15) 

44F8 

EXI (6) 

4570 

EXSI (15) 

45A0 

EU(fl) 

4618 

EV (15) 

4658 

EW(15) 

46D0 

UPR (6*6) 

4748 

RPR(4*A) 

4866 

GCXW(6*15) 

48E8 

EXT(6) 

4BB8 

EXST(15) 

4BE8 

EXTP (6) 

4C60 

EXSTPdS) 

4C90 

GCXWP (6* 15) 

4D08 

GCXSWP ( 15*15) 

4F08 

EMRES(A) 

56E0 

IGDNF 

5700 

IGMNF 

57 0 4 




ACTUAL STATE 2ND MOMENT MATRIX 

ACTUAL 2ND MOMENT MATRIX* STATE/SOL VE-FORS 

NO longer used 

ACTUAL 2ND MOMENT MATRIX* STATE/MEaS. CONS. 
ACTUAL 2ND MOMENT MATRIX* 50LVE-F0R VECTOR 
NO longer USED 

ACTUAL 2ND MOMENT MATRIX* SOLVE-FOR/MS , CN. 
ACTUAL 2ND MOMENT MATRIX* SOLVE-FOR/IGNORES 
ACTUAL 2nd MOMENT MATRIX* MEASUREMENT RESIDUAL 
STM, PARTITION ASSOCIATED WITH IGNORE PARAMETERS 
OBSERVATION MATRIX PARTITION ASSOCIATED WITH 

IGNORE Parameters 
NO longer used 
NO longer used 

ACTUAL 2ND MOMENT MATRIX* MEAS. CONS ./ IGNORES 

NO longer used 

ACTUAL 2ND MOMENT MATRIX* MEAS, CONS. VECTOR 
ACTUAL 2ND MOMENT MATRIX* IGNORE PARAMETERS 
ACTUAL DYNAMIC NOISE CONSTANTS 
ACTUAL MEASUREMENT NOISE CONSTANTS 
ACTUAL MEANS OF INITIAL STATE DEVIATIONS 
ACTUAL MEANS OF INITIAL SOLVE-FOR DEVIATIONS 
NO longer USED 

ACTUAL MEANS OF INITIAL MEAS.CONO, DEVIATIONS 
ACTUAL MEANS OF INITIAL IGNORE DEVIATIONS 
ACTUAL 2nd MOMENT MATRIX* DYNAMIC NOISE 
ACTUAL 2ND MOMENT MATRIX* MEASUREMENT NOISE 
ACTUAL 2nd MOMENT MATRIX* STATE/ IGNORES 
ACTUAL MEANS* UPDATED EST, ERRORS* STATE 
ACTUAL MEANS* UPDATED EST. ERRORS* SOLVE-EORS 
ACTUAL MEANS* PROPAGATED EST, ERRORS* STATE 
ACTUAL MEANS* PROPAGATED EST, ERRORS* SOLVE 
ACTUAL 2ND MOMENT MATRIX* STATE/IGNORE 
BEFORE PROCESSING A MEASUREMENT 
ACTUAL 2ND MOMENT MATRIX* SOLVE-FOR/IGNORE 
BEFORE PROCESSING A MEASUREMENT 
ACTUAL means* MEASUREMENT RESIDUALS 
ACTUAL DYNAMIC NOISE FLAG 
NO longer USED 
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NAME (DIM) 

DISP 

definition 


GUI (SIZE lacs) 

PG(6s6) 

0 

CONTROL 

covariance? 

1 

CXXS6(6? 15) 

120 

CONTROL 

COVARIANCE? 

state/solve-for 

CXUG (6»8) 

3E0 

CONTROL 

covariance? 

NO LONGER USED 

CXVG(6«-15) 

570 

CONTROL 

covariance? 

STATE/MEAS. CONS, 

PSG(15«15) 

■ 840 

CONTROL 

COVARIANCE? 

SOLVE-FOR VECTOR 

CXSUG( 1590) 

F4fl 

CONTROL 

COVARIANCE? 

NO longer used 

CXSVGdSelS) 

1308 

CONTROL 

covariance? 

solve-foh/meas, cons. 

XG(fe) 

TG 

1 AlO 
1 A40 

state vector at TG ■ 

TIME OF LAST GUIDANCE EVENT 

EM (2,6) 
PHIG(696) 

1 A4R 
1 AA8 

NO LONGER USED 

STM FROM INITIAL TIME ON FILE TO TG 


LAGS12 (SUE 4) 

K 0 CURRENT ACCELERATION POINT 


LAGS15 


(SUE 


4) 



lELVN 


ACCUMULATED ACCELERATION POINTS TO BE WRITTEN 


LINKll (SUE 10) 


H 

T 


0 signed STEP51ZE (SEC) 

R TIME UP TO WHICH INTEGRATION HAS PROGRESSED 


LINK.1R (SUE TFO ) 


SXl 

(3) 

0 , 

1ST 

SUM 

FOR 

FQNS OF MOT 

ION 

SX2 

(3) 

1 8 

2ND 

SUM 

FOR 

EONS OF MOT 

ION 

svi 

(3?20) 

30 

1ST 

SUM 

FOR 

VARIATIONAL 

EONS 

SV2 

(3?20) 

210 

2ND 

SUM 

FOR 

VARIATIONAL 

EONS 



3-12 


DISP 


DEFINITION 



NAME (DIM) 


LINK39 JSIZE 10) 


YMDK 

HMSIC 

o oc 

YEAR? MONTH? DAY IN CODE (INITIAL FILE TIME) 
HOUR? MINUTE? SECOND IN CODE (SAME TIME) 

MATRIX (SIZE 
A(3?3) 

1C8) 

0 

CONVERSION* SELENOCENTRIC TO SELENOGRAPH I C 

AOOT (393) 


DERIVATIVE OF A 

B ( 3 « 3 ) 

90 

CONVERSION* EARTH INERTIAL TO EARTH BODY 

C (393) 

D8 

FIXFO 

CONVERSION* MEAN 1950 TO TRUE OF DATE 

GHA 

120 

GREENWICH HOUR ANGLE 

0X0(18) 

128 

NOT USED 

XP 

188 

X POLAR MOTION ANGLE 

YP 

ICO 

Y POLAR MOTION ANGLE 


MEAS 


{SIZE 2EE«?) 


TMN( 1000) 
MCODE ( 1000 » 
NMN 
MCNTR 


0 scheduled times of measurements 

IFAO CORRESPONDING TYPES OF MEASUREMENTS SCHEDULED 

2EE0 NUMBER OF MEASUREMENTS SCHEDULED 

2EE4 NUMBER OF MEASUREMENT TO BE PROCESSED NEXT 


MEQEC (SIZE <>H) 

EQEC(3?3) 0 equatorial TO ECLIPTIC TRANSFORMATION MATRIX 


MISC (SIZE 80 


ACC 

0 

FaCP 

0 

FACV 

10 

aiA (12) 

18 

IDNF 

78 

ICOOR 

7C 

ITR 

80 

IMNF 

04 

ISP? 

88 



NO LONGER USED 
NO LONGER USED 
NO LONGER USED 
NO LONGER USED 

FLAG FOR ASSUMED DYNAMIC NOISE 
NO LONGER USED 
NO LONGER USED 
NO LONGER USED 
NO LONGER USED 


3 - 13 , 


NAME(DIM) DI5P DEFINITION 


NAME (SIZE ?FR) 

FVNM(ll) 

WNNAME ( 12o3) 
CMPNMOO) 

0 

58 

178 

EVENT NAMES 
MEABUPEMENT NAMES 
NO LONGER USED 

NOVENT (SIZE B) 



NEVENT 

0 

NUMBER OF NEXT EVENT 

II 

4 

NO LONGER USED 


OVERPR (SIZE 10) 



MMCODE 

0 

NEXT MEASUREMENT TYPE 

NR 

4 

number of rows in OBSERVATION MATRIX 

TRIM? 

8 

TIMF OF NEXT MEASUREMENT OR EVENT 


OVERZ (SIZE 198) 





RI( 6 ) 

0 

TEMPORARY STORAGE FOR 

STATE 

vector. 

RF ( 6 ) 

30 

TEMPORARY STORAGE FOR 

STATE 

vector 

ADA O 96 ) 

60 

variation MATRIX 



TINJ 

Fn 

NO LONGER USED 



TEVN 

F 8 

TIMf OF NEXT EVENT 



G A ( 3 « 6 ) 

100 

GUIDANCE MATRIX 



IGP 

190 

GUIDANCE POLICY CODE 



NOGFN 

194 

NO longer used 




PHISAV 

(SIZE BAB) 







TOLD 


0 

T1 FOR PHIOLD 





PHIOLD 

(6®20) 

8 

STM from initial 

FILE 

TIME 

TO 

T1 

PHINEW 

(6b20) 

3CP 

STM FROM INITIAL 

FILE 

TIME 

TO 

T2 (NEW TIME) 

FISAVE 

( 61 . 2 O) 

788 

TEMPORARY STORAGE 

FOR 

STM 




\ 
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OISP 


DEFINITION 


NAME (DIM) 


PHT (SIZE 58) 

PLANET(ll) 0 NAMES OF PLANETS 


PSAVE (SIZE UlO) 

PSAV (834) 0 STORAGE FOR COVARIANCES DURING PREDICTION 

AND GUIDANCE EVENTS 


PUNK (SIZE 10) 


IPUN 

0 

PUNCH 

FLAG 

TO 

CORREL 


IPUNE 

4 

punch 

FLAG 

FOR 

EIGENGECTOR 

EVENTS 

IPUNP 

8 

PUNCH 

FLAG 

FOR 

PREDICTION 

EVENTS 

I RUNG 

C 

punch 

flag 

FOR 

guidance 

events 


RLINK4 (SIZE 130) 


PVINT (6) 

0 

PLANET state 

AT GIVEN 

AEINT(6) 

30 

•NOT 

USED 

THIS 

PROGRAM 

SPIN! (6) 

60 

NOT 

USED 

THIS 

PROGRAM 

0BLINT(14> 

90 

NOT 

USED 

THIS 

PROGRAM 

OUMY (6) 

100 

NOT 

USED 

THIS 

program 
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NAME (DIM) 


DISP 


DEFIf'JiTION 


WLINK5 (SIZE 52RD) 



Q(3) 

0 

S/C POSITION 

OD (3) 

18 

S/C VELOCITY 

ZOO(<»053) 

3n 

ARRAY OF S/C ACCELERATION VECTORS 

XV(3«20) 

3F0 

S/C POSITION PARTIALS 

XVD (3»20) 

5D0 

S/C velocity PARTIALS 

WLINK9 (SIZE 60) 



Tt^OPI 

0 

DOUBLE PRECISION VALUE OF 2*PI 

GM ( 1 1 ) 

8 


SCALE (SIZE 8) 



SKALE 

0 

FACTOR USED TO SUBTRACT FRACTION OF KNOWLEDGE 



COVARIANCE FROM CONTROL COVARIANCE IN GUIDM 

SLPOPT (SIZE 38) 



DJ 

0 

EVAL INTERNAL 

IDAY 

8 

DAY OF 1ST RECORD ON EPHEMERIS FILE 

lYEAW 

C 

EVAL internal 

ISPAN 

10 

EVAL INTERNAL 

NBEPM (3) 

1 4 

BODIES FOR POLYNOMIAL COEFFICIENTS 

NOEGRE (3) 

20 

DEGREE OF POLYNOMIALS 

NCFOAY 

2C 

NO, OF DAYS PER CURVE FIT 

ISLP50 

30 

EVAL internal 

N8SLP 

34 

EVAL INTFRNAI. 


SLPREC 

TSEC 

(SIZE 07C) 

0 

SECONDS FROM START OF THIS YEAR TO 

mopt 

PPOLY ( 

3,20s2) 

8 

THIS RECORD TIME INTERVAL 
POSITION polynomial COEFFICIENTS 


VPOLY < 

3^2092) 

3C8 

VELOCITY polynomial COEFFICIENTS 


APOLY { 

3 9 3 c 1 0 ) 

788 

polynomial COEFFICIENTS FOR MATRIX 

A 

CPOLY ( 

393(.i0 

ASH 

polvnomial coefficients for matrix 

C 

PDF.LH ( 

1 n ) 

D?8 

polvnomial coefficients for delta 


IDAY 


D78 

INITIAL DAY OF THIS RECORD 



! 
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/ 





DISP DEFINITION 



STM (SIZE ^CDfl) 


P ( 6 9 6 ) 

CXXS (6o 15) 
CXU(6{.8) 

CXy (69 15) 
PS(15»15) 
CXSU( 15«0) 
CXSV (15® 15) 
UO (8®8) 

VO ( 15« 15) 

PHI (6®6> 
TXXS(6®15) , 

TXU(698) 
Q(6«6) 


H ( 4 9 4 ) 

■AK (6®4) 

S ( 15s4 ) 

H (4»6) 

AM(49l5) 

6(4*8) 

AL (4® 15) 


HPHR 


PP (6?6) 
CXXSP(69l5) 
CXUP(698) 
CXVP (6* 15) 
PSP (15.15) 
CXSUPU5915) 
CXSVP( 15* 15) 


0 covariance matrix® state variables 

120 covariance matrix, state/sol ve-fok 

3F0 covariance MATRIX® NO LONGER USED 

570 COVARIANCE MATRIX® STaTE/MFAS. CONS, 

840 CONVAPIANCE MATRIX® SOLVE-FOR PARAMETERS 

F48 covariance MATRIX® NO LONGER USED 

1308 covariance MATRIX® SOLVE-FOR/MEAS, CONS, 

lAlO COVARIANCE MATRIX® CONTROL VARIABLES (BURN) 

ICIO COVARIANCE MATRIX® MEASUREMENT CONSIDERS 

2318 STATE-TO-STATF TRANSITION MATRIX 

2438 SOLVF-FOR-TO'STaTE TRANSITION MATRIX 

2708 CONTROL“TO“STATE TRANSITION MATRIX 

2888 dynamic NOISE COVARIANCE MATRIX 

29A8 MEASUREMENT NOISE OBSERVATION MATRIX 

2A28 FILTER GAIN MATRIX® STATE PARTITION 

2AE8 filter gain MATRIX® SOLVE-FOR PARTITION 

2CC8 OBSERVATION MATRIX PARTITION® STATE 

2DR8 OBSERVATION MATRIX PARTITION® SOLVF-FOR 

2F68 OBSERVATION MATRIX PARTITION® NOT USED 

3068 OBSERVATION MATRIX PARTITION* MFAS. CONS. 

3240 STORAGE FOR MATRIX AJ WHILE COMPUTING AK 
32C8 COV, MATRIX, STATE® JUST BEFORE MEAS, 

33E8 C0V„ mat,® STATE/SOLVE® BEFORE MEAS. 

36B8 COV, MAT,® NO LONGER USED 

3838 COV. MAT.® STATE/MFAS, CONS, BEFORE MEAS, 

3B08 COV. MAT,® SOLVE-FOR® BEFORE MEAS, 

4210 COV. MAT,® WO longer USED 

45D0 COV. MAT., SOLVE-FOR/ME AS, CONS. BEFORE MEAS. 


STVEC (SIZE 1B8) 


XI (6) 

XF (6) 

X8 (6) 

NDIMl 

NDIM2 

NDIM3 

IAUGIN(24) 


0 / state vector at TRTMl 

30 STATE VECTOR AT TRTM2 

60 NO LONGER USED 

90 DIMENSION OF SOLVE-FOR VECTOR 

94 =0 

98 dimension of MEAS, CONS, VECTOR 
9C INPUT ARRAY OF AUGMENTATION CODES 


TIM (SIZE 30) 


DATEJ 

0 

JULIAN DATE (REFERENCED TO 1950) OF TRTMB 

TRTMl 

8 

TIMF in days (SINCE ZERO) OF LAST 

deltm 

10 

TIME INTERVAL FOR PROPAGATION 

FNTM 

18 

final time 

UNI VT 

20 

UNIVERSAL time 
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\ 


NAMF (DIM) 

DISP 


DEFINITION 

VARMAT (SI^E ACl 





REXVO) 

0 

INPUT IMPULSIVE INSERTION DELTA-V VECTOR 

RADA (3fl6) 

IS 

INPUT variation matrix 

IFVMRI 

AS 

= 0, 

RADA NOT INPUT 



= 1 

RADA : 

INPUT 

VM (SIZE 10^) 





ALN(5TH 

0 

UNITS/AeU 

0 (VALUE SET FOR KM) 

TM 

8 

UNITS/DAV 

(VALUF SET FOR SEC) 

DELTP 

10 

NO 

LONGER 

USED 

RC (6) 

18 

NO 

LONGER 

USED 

DC 

68 

NO 

LONGER 

USED 

RSI (3) 

50 

NO 

LONGER 

USED 

VSK3) 

68 

NO 

LONGER 

USED 

OSI 

80 

NO 

(.0N6ER 

USED 

RVS (6) 

88 

NO 

LONGER 

USED 

VMU 

Be 

NO 

LONGER 

USED 

B 

CO 

NO 

LONGER 

USED 

BDT 

C8 

NO 

LONGER 

USED 

BDR 

DO 

NO 

LONGER 

USED 

delth 

D8 

NO 

LONGER 

USED 

timint 

EO 

NO 

LONGER 

USED 

RE (6) 

£8 

NO 

LONGER 

USED 

RTP (ft) 

lie 

NO 

LONGER 

USED 

CAINC 

168 

NO 

LONGER 

USED 

RCA 

150 

NO 

LONGER 

USED 

taca 

158 

NO 

LONGER 

USED 

SSS(3) 

160 

NO 

l.ONGEP 

USED 

NLP 

178 

LAUNCH PLANET ID (EARTH) 

NEP 

17C 

ephemeris 

Planet id (earth) 

NBOD 

180 

NO 

LONGER 

USED 

NB 

186 

NO 

LONGER 

USED 

NTP 

IBO 

TARGET planet ID (EARTH) 

INPR 

1B6 

NO 

longer 

USED 

IPROB 

188 

PROBLEM NUMBER (INPUT) 

ISPH 

1B6 

NO 

LONGER 

USED 

■INCMT 

ICO 

NO 

LONGER 

USED 

lEPMEM 

1C6 

NO 

LONGER 

USED 

ICL 

ica 

NO 

LONGER 

USED 

IPRINT 

ICC 

EACH IPRINT-TH measurement IS OUTPUT 

ICL2 

ICO 

NO 

LONGER 

USED 
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NAME (DIM) 


DISP 


definition 


XXXL {SIZE 

XSL (15) 
XU(8) 

XV( 155 

XLAB<6) 

XN«(24) 

KPRINT 


2245 


0 SOLVF-FOR PARAMETER NAMES 

78 NO LONGER USED 

B8 measurement CONSIDER PARAMETER NAMES 

130 state vector component NAMES 

160 augmentation parameter labels 

220 =0» PRINT ONLY PHI*P*PHI(T) 

sl» PRINT ALL covariance DATA 
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3«3o? COMMOfsl VfiMABLES IN ALPHABETICAL ORDER 


VARIABLE (DIM) 


A(3^3) 

ACC 

ADA (3 »b) 

ADOT (3«3) 
AEINT (6) 

AINC7 
AK (6v4) 

AL ( A « 1 5 ) 
ALNGTH 
AM(4,lb) 
AN0DE7 
AN ( A 0 15 ) 

APOLY (3o3« 10) 
B(3®3) 

B 

0DR 

RDT 

BIA (12) 

BRNTIM 
C (3o3) 

CAINC 

CECC(^99) 

CINC (4«9) 

CMEAN<499) 

CMIN(6« 15) 
CMPNMOO ) 
C0ME6(4,9) 
C0ME6T (4 »9) 
CPLU(69 15) 
CPOLY(3»39l0) 
CSAX ( 299 ) 
CXSU( 15s,B) 
CXSUG ( 15 sR) 
CXSUP(15»15) 
CXSV (IS, 15) 
CXSVG (15913) 
CXSVP( 15913) 
CXU(693) 
CXUG(6o8) 

CXUP (6 ) 

CXV(6, 15) 

CXVG (6» 15) 
CXVP (69 15) 
CXXS (69 15) 
CXXSG (Go 15) 
CttXSP (bsl5) 


BLOCK DEFINITION 


MATRIX CONVERSION* SELENOCENTRIC TO SEI EN06RAPHIC 

Wise NO LONGER USED 

OVERZ VARIATION MATRIX 

MATRIX DERIVATIVE OF A 

RLINK4 NOT USED THIS PROGRAM 

EVENT NO LONGER USED 

STM FILTER GAIN MATRIX* STATE PARTITION 

STM OBSERVATION MATRIX PARTITION* MEAS. CONS. 

VM UNITS/A. U. (VALUE SET FOR KM) 

STM OBSERVATION MATRIX PARTITION* SOLVR-FOR 

EVENT NO LONGER USED 

GENRL OBSERVATION MATRIX PARTITION ASSOCIATED WITH 
IGNORE Parameters 

SLPREC polynomial COEFFICIENTS FOR MATRIX A 
MATRIX CONVERSION? EARTH INERTIAL TO EARTH BODY 
FIXED 


VM 

NO LONGER USED 


VM 

NO longer used 


VM 

NO LONGER USED 


MISC 

NO LONGER USED 


EVENT 

DURATION (DAYS) OF FINAL INSERTION 

BURN 

matrix 

CONVERSION* MEAN 1950 TO TRUE OF DATE 

VM 

NO LONGER USED 


BLK 

NO LONGER USED 


8LK 

NO LONGER USED 


8LK 

NO LONGER USED 


GAINC 

CORRELATION MATRIX* POS/VEL AND SOLVE-FORS 

name 

NO LONGER USED 


BLK 

WO LONGER USED 


BLK 

NO LONGER USED 


GAiNC 

CORRELATION MATRIX* POS/VEL AND SOLVE-FORS 

SLPREC 

polynomial COEFFICIENTS FOR MATRIX 

C 

BLK 

NO LONGER USED 


STM 

covariance matrix* no longer used 


GUI 

CONTROL COVARIANCE* NO LONGER USED 


STM 

COV. MAT.* NO LONGER USED 


STM 

covariance matrix, solve-for/meas. 

CONS. 

GUI 

CONTROL covariance* SOLVE-FOR/MEAS, 

CONS, 

STM 

COVc MAT.* SOLVe-FOH/MEASo CONS. BEFORE MEAS 

STM 

covariance MATRIX* NO LONGER USED 


GUI 

CONTROL covariance* NO LONGER USED 


STM 

COV. MAT.* NO LONGER USED 



STM covariance matrix? STATE/MEAS. CONS. 

GUI CONTROL COVARIANCE* STATE/MEAS. CONS. 

STM COV, WAT.? STATE/MEAS. CONS. BEFORE MEAS. 

STM COVART aNCF ' waTPIxT ST 0 0 

GUI CONTROL covariance* STaTE/SGLV£-FOP 

STM COVc mat.? STATE/SOLVE* before MEAS. 
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VARIABLE (DIM) 

block 

DEFINITION 

DATEJ 

TIM 

JULIAN DATE (REFERENCED TO 1950) OF TRTM8 

DC 

VM 

NO LONGER USED 

delth 

VM 

NO longer used 

deltm 

T IM 

TIME INTERVAL FOR PROPAGATION 

OELTP 

VM 

NO LONGER USED 

OJ 

SLPOPT 

EVAL internal 

DNCNO) 

const 

dynamic noise CONSTANTS 

DSI 

VM 

NO LONGER USED 

DUMY (6> 

RLINK4 

NOT USED THIS PROGRAM 

DV8 (3) 

EVENT 

NO longer used 

DXQI18) 

matrix 

NOT USED 

ECC7 

EVENT 

NO LONGER USED 

EE(A) 

GEN6D 

ACTUAL MEANS OF EXECUTION ERROR PARAMETERS 

EEE (A) 

GENGD 

VARIANCES OF EXECUTION ERROR PARAMETERS 

EM(2«6) 

GUI 

NO LONGER USED 

EMN (15) 

BLK 

NO longer USED 

EMRES(4) 

GEWRL 

ACTUAL MEANS* MEASUREMENT RESIDUALS 

EMI 

DPNUM 

double-precision value OF UE-l 

EM2 

DPNUM 

double-precision value of 1oE-2 

EM3 

DPNUM 

double-precision value of UE-3 

EMA 

DPNUM 

double-precision value of ue-a 

EM5 

DPNUM 

DOUBLE-PRECISION VALUE Of 1 ,E-5 

EM6 

DPNUM 

DOUBLE-PRECISION VALUE OF 1 oE-6 

EM7 

DPNUM 

double-precision VALUE OF 1 <,£-7 

EM8 

DPNUM 

DOUBLE-PRECISION VALUE OF UE-8 

EMQ 

DPNUM 

double-precision VALUE OF l.E-9 

EM13 

DPNUM 

DOUBLE-PRECISION VALUE OF UE-13 

EPS 

CONST 

EARTH OBLIQUITY 

EP50 

DPNUM 

DOUBLE-PRECISION VALUE OF l.E*50 

EQEC(3«3) 

MEQEC 

equatorial to ECLIPTIC TRANSFORMATION MATRIX 

EU (B) 

GENRL 

NO LONGER USED 

EV (15) 

genrl 

ACTUAL MEANS OF INITIAL MEASeCONO® DEVIATIONS 

EVNM(ll) 

name 

■ EVENT NAMES 

EWI15) 

GENRL 

ACTUAL MEANS OF INITIAL IGNORE DEVIATIONS 

EXI (6) 

GENRL 

actual means OF INITIAL STATE DEVIATIONS 

EXSI ( 15) 

GENRL 

ACTUAL MEANS OF INITIAL SOLVE-FOR DEVIATIONS 

EXST (15) 

genrl 

ACTUAL MEANS* UPDATED ESTc ERRORS* SOLVE-FORS 

EXSTP (15) 

GENRL 

ACTUAL MEANS* PROPAGATED EST* ERRORS* SOLVE 

EXT(5) 

GENRL 

ACTUAL MEANS* UPDATEU ESTo ERRORS* STATE 

EXTP(6) 

GENRL 

ACTUAL MEANS* PROPAGATED EST= ERRORS* STATE 

FACP 

MISC 

NO longer USED 

FACV 

MISC 

NO LONGER USED 

EISAVE (6520) 

PHISAV 

TEMPORARY STORAGE FOR STM 

ENTM 

TIM 

FINAL TIME 

FOP 

C0NST2 

OFF-DIAGONAL ANNIHILATION VALUE (POSITION) 

FOV 

C0NST2 

OFF-DIAGONAL ANNIHILATION VALUE (VELOCITY) 

G ( A 5 8 ) 

STM 

observation MATRIX PARTITION* NOT USED 

GA(3«6) 

OVER2 

guidance MATRIX 
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VARIABLE (DIM) 


BLOCK 


DEEINITION 


GCXSUG ( 15»8) 

gengdi 

ACTUAL CONTROL SECOND MOMENT MATRICES 

NO longer used 

6CXSVG ( 15« 15) 

6ENGD1 

SOLVE-FOR/MEASUREMENT CONSIDERS 

GCX5WG ( 15-) 15) 

gengoi 

solvE-for/ignore Parameters 

GCXUG ( 698 ) 

GENGOl 

NO LONGER USED 

6CXVG(6.15) 

GENGDl 

STaTE/MEaSUREMENT C0NSI[)ERS 

GCXWG ( 6 « lb) 

GEWGDl 

state/ignore parameters 

GCXXSG( 6 *lb) 

GENGDl 

ST ATE/SOLVE-FOR VECTOR 

GPG ( 6 e fc ) 

GENGDl 

state 

GPSG( lb, 15) 

GENGDl 

SOLVE-FOR VECTOR 


GCUV (8. 15) 

GENRL 

NO longer used 

GCUW (B« 15) 

GENRL 

NO longer used 

GCVVK(15,15) 

GENRL 

ACTUAL 2ND MOMENT MATRIXs MEAS, CONS. /IGNORES 

GCXSU ( 15«8) 

GENRL 

NO LONGER USED 

6CXSV ( 16« lb) 

GENRL 

ACTUAL 2ND MOMENT MATRIXs SOL V F-KOR/M S , CN, 

GCXSW ( 15s 1 b) 

GENRL 

ACTUAL 2ND MOMENT MATRIXs SOLVE-FOR/IGNORES 

6CXSWP ( 15s lb) 

GENRL 

ACTUAL 2ND MOMENT MATRIXs SOLVE-FOR/ IGNORE 
BEFORE PROCESSING A MEASUREMENT 

6CXU(6s3) 

GENRL 

NO LONGER USED 

GCXV (6s 15) 

GENRL 

ACTUAL 2ND MOMENT MATRIXs STATF/MFAS. CONS. 

. GCX><^(6,15) 

GENRL 

ACTUAL 2ND MOMENT MATRIXs STATE/IGnORES 

GCXWP (6t 15) 

GENRL 

ACTUAL 2nd MOMENT MATRIXs S T A T E/ IGNORE 
BEFORE PROCESSING A MEASUREMENT 

GCXXS (b» 15) 

GENRL 

ACTUAL 2ND MOMENT maTRIas STaTF/SOI VE-FORS 

GDNCN(3) 

GENRL 

ACTUAL DYNAMIC NOISE CONSTANTS 

6 HA 

MATRIX 

GREENWICH HOUR ANGLE 

GM(1 1 ) 

RLINK9 


GMNCN(15) 

GENRL 

ACTUAL MFASUPEMENT N-ISE C-NSTAN7S 

GP ( 6 s 6 ) 

GENRL 

ACTUAL STATE 2ND M-MENT -♦T-*0 

GPS (I5sl5) 

GENRL 

ACTUAL 2ND moment MATRIXs SOLVE-FOP VECTOR 

6U (Bs8) 

GENRL 

NO LONGER USED 

GV(159 15) 

GENRL 

actual 2ND MOMENT MATRIX', M£AS. CONS, VECTOR 

GW(15s 15) 

GENRL 

ACTUAL 2ND MOMENT MATRIX, IGNORE PARAMETERS 

H 

LINK 11 

signed STEPSrZE (SEC) 

H(A, 6 ) 

STM 

OBSERVATION MATRIX PARTITION, STATE 

HALF 

DPNUM 

double-precision Value of half (O.s) 

HMSIC 

LINK39 

HOUR, minute, second in code (Same time) 

HPHR 

STM 

STORAGE FOR MATRIX AJ WHILE COMPUTING AK 

HPT 

EVENT 

NO longer used 

IAUGIN(P<>) 

STVEC 

INPUT ARRAY OF AUGMENTATION CODES 

lAUGW (E^) 

GCA 

IGNORE parameter AUGMENTATION VECTOR 

IC0Q3(E0) 

EVENT 

NO 1 ONGER USED 

ICOT3(EO) 

EVENT 

CODES FOR guidance POLICIES 

ICENTB 

FLAGSi? 

CENTRAL BODY NUMBER 

ICL 

VM 

NO longer used 

ICL? 

VM 

NO LONGER USED 

ICOOR 

Wise 

NO longer used 

IDAY 

slpopt 

DAY OF 1ST RECORD ON EPHEMFRIS FILF 
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VARIABLE 


BLOCK 


definition 


IDAY 

5LPREC 

INITIAL DAY OF THIS RECORD 

IDNF 

Wise 

FLAG FOR ASSUMED DYNAMIC NOISE 

lEIG 

EVENT 

NO LONGER USED 

lELVN 

LAGS15 

ACCUMULATED ACCELERATION POINTS TO BE WRITTEN 

lEPHEM 

VM U 

NO LONGER USED 

IEVNT(501 

EVENT 

COOFD EVENT TYPES CORRESPUNUiNG TO REV TIMES 

IFVMRI 

VARMAT 

=0» RADA NOT INPUT 

= U RADA INPUT , 

IGDNF 

GENRL 

actual DYNAMIC NOISE FLAG 

IGEN 

GCA 

=09 perform no generalized COVARIANCE ANALYSIS 
=1* PERFORM generalized COVARIANCE ANALYSIS 

IGHNF 

GENRL 

NO LONGER USED 

IGP 

OVERZ 

GUIDANCE POLICY CODE 

IHYPl 

EVENT 

NO longer used 

II 

NOVENT 

NO longer used 

I IPOL 

EVENT 

NO longer used 

IMNF 

MISC 

NO LONGER USED , 

INC 

COM 

NO LONGER USED 

INCMNT 

COM 

NO LONGER USED 

INCMT . 

VM 

NO LONGER USED 

INCPR 

COM 

NO longer used 

IND(25) 

FLAGS2 

CURRENT SECTION FLAGS 

INDSECaO»3) 

FLAGS2 

SECTION FLAGS 

INPB 

VM 

NO LONGER USED 

IPRINT 

VM 

EACH IPRINT-TH MEASUREMENT IS OUTPUT 

I0P7 

EVENT 

NO longer used 

IPOL 

EVENT 

NO LONGER USED 

IPR 

COM 

NO LONGER USED 

IPR08 

VM 

problem number (INPUT) 

IPRT (4) 

COM , 

NO longer used 

IPUN 

PUNK 

punch flag to correl 

IPUNE 

PUNK 

PUNCH FLAG FOR EIGENGECTOR EVENTS 

IPUNG 

PUNK 

punch flag FOR GUIDANCE EVENTS 

IPUNP 

PUNK 

PUNCH FLAG FOR PREDICTION EVENTS 

ISLP50 

SLPOPT 

FVAL INTERNAL 

ISPAN 

SLPOPT 

EVAL internal 

iSPH 

VM 

NO LONGER USED 

ISP2 

MISC 

NO longer used 

ITR 

Wise 

NO LONGER USED 

ITRAT 

COM 

NO longer USED 

I YEAR 

SLPOPT 

EVAL INTERNAL 

jPR(4siA) 

GENRL 

ACTUAL 2ND MOMENT MATRIX* MEASUREMENT RESIDUAL 

K 

LAGS12 

CURRENT ACCELERATION POINT 

KOUNT 

COM 

NO longer USED 

KPRINT 

XXXL 

=0, PRINT only PHI*P«PHI(T) 
=1» PRINT ALL COVARIANCE DATA 

MCNTR 

MEAS 

NUMBER OF MEASUREMENT TO BE PROCESSED NEXT 

MCOOE 1 1000) 

MEAS 

CORRESPONDING TYPES OF MEASUREMENTS SCHEDULED 

MMCODE 

OVERPR 

NEXT measurement TYPE ' 
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y/ARTABLF (DIM) 

BLOCK 

DEFINITION 

HNCN (12) 

CONST 

MEASUREMENT NOISE CONSTANTS 

MNNAMP ( 12 »3) 

name 

MEASUREMENT NAMES 

MUPLAN (11) 

BLK 

GRAVITATIONAL CONSTANTS ( KM ( ( 3 /SE C** 2 ) 

NAF 

event 

NO^LONGER USEt) 

NAF 6 ( 20 ) 

EVENT 

NO LONGER USED 

NB 

VM 

NO LONGER USED 

NROO 

VM 

NO LONGER USED 

NBSLP 

SLPOPT 

EVAL internal 

NCFDAY 

SLPOPT 

NOo OF DAYS PER CURVE FIT 

NOEGHE (3) 

SLPOPT 

DEGREE OF POLYNOMIALS 

NBEPM (3) 

SLPOPT 

BODIES FOR POLYNOMIAL COEFFICIENTS 

NROOY 

COM 

NO longer used 

NBODYI 

COM 

NO longer used 

NDIMl 

STVEC 

DIMENSION OF SOLVE-FOH VECTOR 

NDIM2 

STVEC 

-0 

ND1M3 

STVEC 

dimension OF MEASci CONS. VECTOR 

ND1M4 

GCA 

DIMENSION OF IGNORF PARAMETER STATE 

NEP 

VM 

EPHEMERIS planet id (EARTH) 

NEV 

EVENT 

number of events scheduled 

NEVFNT 

NOVENT 

number OF NEXT EVENT 

NEVl 

EVENT 

number of scheduled EIGENVECTOR EVENTS 

NEV2 

EVENT 

number of SCHEDULED PREDICTION EVENTS 

NEV3 

EVENT 

NUMBER OF scheduled GUIDANCE EVENTS 

NEVA 

EVENT 

NUMBER OF scheduled INSERTION EVENTS 

NEV5 

EVENT 

NO LONGER USED 

NFV 6 

EVENT 

NO longer USED 

NEV7 

EVENT 

NO LONGER USED 

NEVA 

event 

NO LONGER USED 

NEV9 

EVENT 

NO LONGER USED 

NEVl 0 

EVENT 

NO longer used 

NEVl 1 

EVENT 

NO longer used 

NGF 

EVENT 

COUNT OF GUIDANCE EVENTS 

NLP 

VM 

LAUNCH planet ID (EARTH) 

NMN 

MEA5 

number of measurements SCHEDULED 

NO( 1 1 ) 

BLK 

NO longer USED 

NOGEN 

OVEP7 

NO LONGER USED 

NPE 

EVENT 

count of predicti-n EVE-TCi 

NQE 

EVENT 

MO LONGER USED 

NR 

OVERPR 

NUMBER OF ROWS IN OBSERVATION MATRIX 

NSECTN 

FLAGS2 

CURRENT SECTION NUMBER 

NST 

CONST 

NUMBER OF STATIONS TO BE USED (MAXIMUM 

NTP 

VM 

target Planet id (earthj 

NTSFQS 

FLAGS2 

total number of sections 

OBLINT (14) 

RLINK4 

NOT USED this PROGRAM 

OMEGA 

CONST 

rotation RATE OF EARTh 

ONE 

DPNUM 

double-precision VAL(JE of one (1.0) 

P ( 696 ) 

STM 

covariance matrix, state variables 

PDFLH(IO) 

slprec 

polynomial coefficients for delta fi 

PERP? 

EVENT 

NO longer useo 
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VAPIflRLF (DIM) , BLOCK 


OLFINl TIO^ 


PHI (htb) 

PHIG (G»6) 
PHINFW(h.?U) 
PHI(M 0(B = ?CN 
PI 

P1.4MFT (11) 
PMftSS ( 11 ) 
PMIN(6*6) 

PP ((S«G) , 

PPLU (6«6) 
PPOLKO^^O^S) 
PS ( IS? IS) 

PSAV (H3:») 

PSG ( IS? IS) 
PSMTNf IS? IS) 

PSP ns, 15) 

PSPLil(l5?ls) 
PVINT (6) 

P7 

0(S) 

:) ( S , s ) 

OQ ( S) 

QPR (6,S) 

P ( 4 5 4 ) 

PAD 

p fi n A ( s ? s ) 
PAIITIIS (11) 

PC (^) 

PCA 
PF (S) 

PE«V (3) 

PF (S> 

PI (S> 

PMASS Ml) 

BPP(4,4) 

PSA'/F ( 6 ) 
PSI(3) 

PTP (S) 

PVS (S) 

S ( 1 5,4) 

SAL (3) 

S T G A I P 
SIGRFt 
SIGPPO 
SIGPF.S 

SKfll F 

SLAT ( 3) 

SL ON ( 3 ) 


GUI 

STM 

GUI 

PHISAV 

P H I 5 A V 

COM 

PPT 

BLK 

GAINC 

STM 

GA INC 

SLPRFC 

STM 
PSA VP 

GUI 
GAiNC 
STm 
GAINC 
PL I NK 4 
FVFNT 
PL IN<5 
ST'« 

PL Inks 
GFNHL 
STM 
COM 

V APM4 T 
HLK 

VM 

VM 

VM 

VAPMA T 

OVERZ 

OVFRZ 

BLK 

GFNPI. 

GAiNC 

VM 

VM 

VM 

STm 

CONST 

E'VFNr 

tVENl 

event 

EVENT 

scale 

const 

CONST 


CONTROL rOvAPlANCF, STaTE 

ST ATF-T,0-ST ATE TRANSITION MATpIX 

STm from Initial IImE on FILE TO TG 

STM F^OM initial FILL TTMt JQ TP (NEW TIME) 

STm FpoM initial FILE TIME TO II 

MATHEMATICAL constant PI 

NAMFS OF Pt. ANKTS 

gravitational constants ( AcO,^»»3/OAy<‘*2) 

POS/VFL COVAPIANCE BEFORE MEASUREMENT (*«LS) 
COv, maTPTX, STATr? Ji)ST BEPOPF. MEAS, 

POS/VEL covariance AFTER MEASUPEMEnT (WLS) 
POSITION POLVNomiaL coefficients 
COMVAPIANCF MATRIX, >OLvE-EON PARAMETERS 
STORAGt F-P covariances OUPINB PPEOICTl-N 
AND GUIDANCE EVENTS 

CONTROL COVARIflNCF? SOLVE-EOP VECTOR 
SOlvE-FOR covariance PEfOPF. measurement (WLS) 
COV, »«AT,, SOLVE-EOP? BEFORE MEAS. 

SOlvE-EOR COVARIAijCE after measurement (WLS) 
planet state, at givEn time 
NO LON(3ER USED 
S/C POSIT I DM 

□ VNftMiC NOISE covariance. MATRIX 
s/c velocity 

ACTUAi. ?ND MOMENT matrix? DYNAMIC NOISE 
measurement noise observation matrix 

N'JMMEP OF' degrees pEp RAOIAN 

INPltT variation matrix 

RAOM OF THE planets (A„U.) 

NO I OnGFP iJSE. O 
NO I, OnGFP used 
no longer USEO 

INPUT IMPULSIVE INSERTION IJtLTA-V VECTOR 
temporary storagf eop state* vector 
TFmdoraRY STOWAGE FOf^ bTATE VECTOR 
GRAVITATIONAL COMgTANTS (RELATIVE 10 SUN) 
actual PNO .moment MAIREX? MEASUREMENT NOISE 
state vector at TlAST 
NO longer used 
NO longer used 
NO longer USEO 

filter gain matrix? SOLVE-fOk partition 
station altitudes (ABOVE RADIUS OF EARTH) 
variance of error in pointing angle 1 

VARTANCE OF error IN POINTING ANCjLE 2 
VAPfANCF UE proportionality ERROR 
variance: of PESOLUTIOIi ERROR 

factor USE!.! TO SUHTRACI FRACTION OF KNOWLEDGE 

covariance from control covariance in guidm 

STAT ion I. ATTTUDES 
station longitudes 
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VARIABLF (DIM) 


BLOCK 


OEFir^ITIOM 


SMjP (?«q) 
SPHFWF (11) 
SPIN! (6) 
SSS(3) 

-SVl ( 3»E0) 
SV?(3,i!0) 
SXl ( 3) 
5X?(3) 

T 

T 

TACA 
TAU7 
TFV (50 ) 
TEVN 
TG 

THRFE 
T Imint 

TINJ 
TL AST 
TM 

TMN( 1 000) 

told 

TPT2 ( ?0 ) 

TKTMl 

TPTw? 

TSFC 


BLK 

RLK 

PLINK4 

VM 

LINrv3a 

LINK38 

LI.MK38 

BLK 

LlNKll 

VM 

event 

EVENT 

OVFRZ 

GUI 

DPNUM 

VM 

OVERZ 

GAINC 

VM 

ME AS 
PHISAV 
EVEMT 
TIM 

OVERPR 

SLPREC 


T olO 
TWOPI 
TWOPI 
TXU (6»B) 

TXW (6. IS) 
TXXS(6« 15) 

UNI VT 
UST (3) 

UO ( B»R) 

VMLI 

VPnLY(39?0f2) 
VST (3) 

VO ( IS 9 IS) 

VST (3) 

WST ( 3) 

XP (S) 
aF (6 ) 

XG (S) 

XI(4) 

XIG( IS) 
XLAR(6) 

XNM ( ) 

XP (S) 


OPNUM 

DPNUM 

RLINKB 

STM 

genrl 

STM 

TIM 

CON5T? 

STM 

VM 

SLPREC 

const? 

STM 

VM 

CONST? 

STVEC 

STVEC 

GUI 

STVEC 

GCA 

XXXL 

XXXL 

BLK 


NO LONGER USED 

SPHFPES of INFLUFNCE (A.U.) 

NOT USED THIS PROGRAM 

NO longer used ■ 

1ST SUM for variational eons 
?N o SUM for variational eons 

1ST SUM FOR FQNS OF MOTION 
?ND SUM FOR FONS OF MOTION 
trajectory time in UAYS 

time up TO which integration has progressed 
NO longer used 

NO LONGER USED 

SCHEDULED times OF EVENTS 

time of mex t event 

time of last guidance event 

DOURLE-PREC I SIDN VALUE OF THREE (3.0) 

NO I ONGE- USFD 
NO LONGER USED 

time WHEN measurement laST PROCESSED 
UNITS/DAY (VALUE SET FOR SEC) 

SCHEDULED times OF MEASUREMENTS 
T1 FOR PHIOLO 

timks “pedictfd to in prediction Events 

time in days (SINCE ZERO) OF LAST 

TIME OF NFXT MFASUREMENT OP EVFNT 

SECONDS FROM START UF THIS YEAR TO MDPT OF 

THIS RFCORD TIME INTERVAL 

DOUBLE-PPFCI S ION VALUE OF TWO {?,0) 

double-precision VALUE OF ?,**PI 

double precision VALUE OF ?.<*PI 

CONTROL-TO-STATF TRANSITION MATRIX 

STM PaRIITION ASSOCIATED WITH IGNORE PARAMETERS 

solvf-fop-to-state Transition matrix 
universal time 

X-ITIPFCTION cosine for stars 

covariance matrix* control variables (BURN) 

NO LONGER USFD 

velocity polynomial coefficients 
y^direction cosine for stars 
covariance matrix, measurement considers 

NO LONGER USED 

Z-DIRECTION COSINE FOR STARS 

NO tONGFR USFD 

STATE VECTOR AT ThTM? 

STATE vector at TG 
STATE VECTOR AT TRTmI 
IGNORE PARAMETER LAHLE S 
STATE vector COMPONENT NAMES 
AUGMENTATION PARAMETER LABELS 
STATE VECTOR OF PLANET 
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VartiARLF (0 IM) 


hL0Ct<. 


uEFINITIOM 


XP MCTkIX X P0L6W MOTIOf^l flNJ'iLt 

X5L(1‘^) XXXL SOLVf-'-FOW PARAMf:TtP NAV^S 

XU(fl) XXXL NO LONtiER 0'=^FI' 

XV(^,20) RLINKS '^/r POSITI'li'l P6PTIALS 

XV (IS) XXXL measurement consider parameter names 

xvn(3^?0) RLTNKS S/C VELOCITY PARTIALS 

YMOIC LINK3R YEAt>, MOMTS» DAY IN COf)E (INITIAL FILE TIME) 

YP matrix Y polar motion AN(.H.E 

ZnD(A(l»3) RLINKb ARRAY OF S/C ACCELERATION VECTORS 

ZERO DRNUM onURLE-RRFCiSION VALUE OF ZERO (0.0) 
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4. INDIVIDUAL SUBROUTINE DOCUMENTATION 

This chapter provides the detailed documentation of the subroutines 
comprising the STEAP-L programs. The subroutine hierarchy for NOMINAL 
and ERRAN are defined in Figures 2. 1 and 3. 1, An alphabetical listing 
of all subroutines appearing in STEAP-L is given in Table 4.1 with a 
description of the program(s) that use the subroutine. Subroutines 
that are part of the Goddard Trajectory Determination System (GTDS) 
and are unchanged are not documented in this report as existing doc- 
umentation is available at GSFC. Table 4.2 lists Che purposes of the 
documented subroutines (again in alphabetical order) for convenient 
cross-reference. 

The following pages then detail each subroutine in alphabetical 
order. The level of documentation of the subroutines is based on their 
complexity. Simple utility routines are described by defining their 
purpose, call sequence, and input and output arguments. The documenta- 
tion of more complicated routines defines local and common variables 
computed by the routines, mathematical analysis, and flowcharts. 



TABLE 4.1 ALPHABETICAL LISTING OF SUBROUTINES 
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GTDS ROUTINE (NOT DOCUMENTED) 



TABLE 4.2 PURPOSES 


NAME 

ACOSH 

ATCEGV 

BLOCK DATA 
ERRA.N 

BURN 

BURNV 

CALJUL 

CAREL 

C0RREL 

C0VMAT 

COWELL 

C START 

DABSV 

DANGMD 

DANGV2 

DATA 

DAVECT 

DDOT 

DDOTB 

DMADD 

DMATPY 

DMSUB 

DRD 


PURPOSE 

To Compute the Hyperbolic Arc:-Cosine 

To Compute Eigenvalues and Eigenvectors of Actual Target 
Condition and Moment Matrix 

To Load Constants into Common Locations used in Various 
other Parts of the Program 

To Compute the Current Acceleration due to a Finite Burn 
To Compute the Finite Burn Partials 

To Compute Julian Date from Calendar Date or Vice Versa 

To Transform Cartesian Coordinates to Conic Elements 

To Convert Covariance Matrix Partitions to Correlation Matrix 
Partitions and Standard Deviations and Write them out 

To Convert Standard Deviation/ Correlation input to Covariance 
Form j 

To Control the Integration Logic after the Integrator 
has been Initialized 

To Read the Header Record on the Sequential Orbit File 
To Calculate the Magnitude of a Vector 
To Modularize an Angle on Two PI 

Calculate a Directed Angle Between Two Vectors in 3 space 

To Read input Data, set Default Values, Initialize and Set 
Internal Parameters and Print Initial Conditions 

Vector Addition 

Vector DOT Product 

To Return the DOT (Inner) Product of Two 3-Vectors 
Matrix Addition 
Matrix Multiplication 
Matrix Subtraction 

To Compute Latitude and Longitude of a Vector 


/ 
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NAME 
DSHIFT 
' DSVECT 
DUNIT 
DUXV 
DVCOMB 
DVECRD 

DVSDIV 

DVSMLT 

DVMAG 

DVSTAT 

DXB 

DYN0 

DZERO 

ECOMP 

EIGHY 

EPHGT 

ERRAN 

EULMX 

GAINl 

GDATA 

GENGID 

GETGOW 

GHA 

GIDANS 


PURPOSE 

To Shift one Vector into Another 
Vector Subraction 
To Unitize a Vector 
To Form a Vector Cross Product 

To Combine (add) Two Vectors, Each Multiplied by Scalars 

To Compute a Unit Vector- from its Right Ascention and 
Declination 

Vector Scalar Division 

Vector Scalar Multiplication 

Calculate the Magnitude of a 3-Vector 

To Compute and Print Statistical Delta-V Parameters 

Calculate the Cross Product of Two 3-Vectors 

To Compute Assumed and Actual Dynamic Noise Covariance 
Matrix 

To Generate a Zero Vector 

To Compute Differential Transformation Relating Target 
Variables to State 

To Control the Computation of Eigenvalues and Eigenvectors 

To Retrieve from the Direct Access SLP File, the State Vector 
of a Planet with Respect to the Sun at an Arbitrary Julian Date 

To Control the Computational Flow of the Basic Cycle 

To Compute Transformation Matrices 

To Compute the Kalman Gain Matrices 

To Initialize Generalized Covariance Quantities 

To Generate Statistics Relating to Actual Guidance Events 

To Generate a State Vector at a Requested Time 

\ 

To Compute the Greenwich- hour Angle and Universal Time 
Dummy Link with Non- Halo Orbit Options 
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NAME 

GNAVM 

GPRINT 

GQC0MP 

GUID 

GUIDM 

HGIDNS 

HLAUCH 

HPRELM 

HTRJTY 

HZERIT 

INITLC 

JAC0BI 

LAMBRT 

LOOP 

MATIN 

MEAN 

MEN0 

MAIN 

M0MENT 

NTM 


PURPOSE 

To Propagate Covariances Between Measurements and Events 
and to Update them at Measurements 

To Print Actual Estimation Error Statistics 

To Confute Actual Execution Error Statistics 

To Compute the Variation, Guidance and Target Condition 
(BEF0RE) Covariance Matrices 

To Control the Execution of a Guidance Event 

To Compute the Change Required to the Control Variables 
for Targeting 

To Compute the Injection Time 

To Initialize Constants and Default Values, Read Input 
Data, and Calculate the Zero Iterate Guess 

To Control the Trajectory Generation Phase 

To Compute the Initial for Targeting when IZERO = 6 or 7 

To Initialize Constants 

To Calculate the Eigen-Values and-Vectors of a Real 
Symmetric Matrix 

To Solve Lamberts Problem for Transfers less than Two PI 
To Solve Lamberts Problem for Transfers Greater than Two PI 
To Compute the Inverse of a Matrix 

To Propagate and Update the Means of Actual State and 
Parameter Deviations and Estimation Errors 

To Compute Assumed and Actual Measurement Noise Covariances 
Entry Point to Program NOMNAL 

To Convert Covariance Matrix Partitions to Correlation 
Matrix Partitions, Calculate Eigen-Values and -Vectors and 
to Print 

To Read the Trajectory File and Manipulate State 
Transition Matrices 
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NAME 

ORBINT 

ORBEND 

0RBWRT 

PECEQ 

PRELIM 

PRED 

PRINT3 

PSIM 

P START 

SAVMAT 

SCHED 

SETEVN 

SHIFT 

SL0020 

SETl 

STAPRL 

steape' 

STMPR 

STVCPR 

SYMTRK 

SYMTRZ 

TIME 

TINE 


PURPOSE 

To Initialize the Sequential Orbit File with Partials 

To Write a ’Final’ Record to the Sequential Orbit File 

To Write Records to the Sequential Orbit File 

To Compute the Ecliptic to Equatorial Transformation Matrix 

Dummy Link with Non- Halo Orbit Options 

To Make Prediction Event Calculations 

To Print Measurement Information 

To Calculate the State Transition Matrix from T2 to T3 
using the T1 to T2 and T1 to T3 State Transition Matrices 

To Initialize the State Partial Matrix 

To Store a Vector P in a Vector PI 

To Order the Measurement Schedule 

To Control All Event Calculations 

To Shift a Double Precision Array to Another Location 
To Minimize F(X) 

To Initialize the Flags for use by Integration Routines 
To Compute Station Location ’Partials 

' I 

To Control the Error Analysis Mode of STEAP 
To Print the State Transition Matrices 
To Print the State Vector in Several Coordinate Systems 
To Symmetrize a Square Matrix 

To Fill the Upper-Right Triangle of a Symmetric Square 
Matrix whose Lower-Left Triangle was input 

To Convert a Time in Seconds to Days, Hours, Minutes and 
Seconds 

To Compute the Julian Date Relative to 1900 from the Calendar 
Date or Vice Versa 
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NAME 


PURPOSE 


TRAKM 

TRNSPS 

TRJTRY 

ZERMAT 


To Compute the Augmented Observation Matrix Partitions 
To Form the Transpose of a Matrix 
Dummy Link with Non- Halo Orbit Options 
To Zero a Matrix 


/ 
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ACOSH-A 


function ACOSH 

PURPOSES TO COMPUTE THE HYPERBOLIC ARC-COSINE 
calling SEOUENCE: RES=AC0SH<X) 

ARGUMENTS! 

X I VALUE OF HYPERBOLIC COSINE 

ACOSH 0 hyperbolic ARC-COSINE OF X 

SUBROUTINES REQUIRED: 

none 


SUBROUTINE ATCEGV 

PURPOSES TO COMPUTE EIGENVALUES AND EIGENVECTORS OF ACTUAL TARGET 
CONDITION 2ND MOMENT MATRIX 

CALLING SEQUENCES CALL ATCEGVail* ATC»EOT,FOV> 

ARGUMENTS* III I NUMBER OF ROWS IN ATC MATRIX 

ATC I ACTUAL TARGET CONOITIGN MATRIX 

EOT I ACTUAL TARGET STATE DEVIATION MEANS 

FOV I FINAL OFF-DIAGONAL ANNIHILATION TERM FOR JACOBI 

SUBROUTINES SUPPORTED* GENGIO 

subroutines required* eighy 

COMMON USED* 

local symbols* OUMt OUTPUT MATRIX FOR JACOBI 

EGVL EIGENVALUES 

PEIG INTERMEDIATE ARRAY 

ROW INTERMEDIATE VECTOR 

S ATC COVARIANCE ARRAY<3»3) 

SOUM ATC covariance ARRAY«2,2) FOR JACOBI 
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blkdat-a 


BLOCK DATA - ERRAN 

HURPOSt-I; TU load CONSTANTS INTO COMMON LOCATIONS USED IN VARIOUS 
OTHER HARTS OF' THE PROGRAM, 

CALLING SEUUFmCE: NONE 


argument: none 


SUBROUTiNEb SUPPORTED: 
SURROUTiNEb REOUIRED: 

HALF THE SUBROUTINES USE THE CONSTANTS 
STORED BY THIS BLOCK DATA 

NONE 

COMMON loaded: XLAB 

XNM 

EVNM 

MNNAME 

OMEGA 

PI 

HAD 

RMASS 

radius 

MUPLAN 

PMASS 

PLANET 

Tm 

SI GRES 

SIGPRO 

SIGALP 

bIGBET 

MNCN 

NLP 

NEP 

NTP 

ALNGTH 

TWOPI 

DNCN 

IGA In 

IGEN 

EM13 

EPSO 

IDNF 

IPR08 

IF VMRI 

ZERO 

ONE 

TWO 

HALF 

THREE 

, EMl 

EM2 

Em3 

EM4 

EM4 

EM5 

EM6 

SLON 

Em7 

EM8 

FM9 

SAL 

SLAT 



BLKDAT-1 


The subprogram BLOCK DATA loads constants into common blocks used by the sub- 
routines in the ERRAN program. 

The arrays loaded are: 


XLAB 

XNM 

EVNM 

MNNAME 

OMEGA 

PI 

RAD 

RMASS 

RADIUS 

MUPLAH 

PMASS 

PLANET 

TM 

SIGRES 

SIGPRO 

SIGALP 

SIGBET 

MNCN 

NLP 

NEP 

NTP 

ALNGTH 

TWOPI. 

DNCN 

IGAIN 

IGEN 

EMI 3 

EP50 

IDNF 

IPROB 

IFVMRI 

ZERO 

ONE 

TWO 

HALF 

THREE 

EMI - EM9 

SAL 

SLAT 

SLON 


1 


Hollerith names of the state-vector components 
Hollerith names of the augmentation parameters 
Hollerith names of the event types 
Hollerith names of the measurement types 
Mean sidereal rate (radians per day) 

IT (conversion factor) 

Degrees per radian (conversion factor) 

Mass ratios of the planets (mass of the sun = 1.0) 
Planet radii in A.U. 

Gravitational constant times planet mass (km^/sec^) 
Gravitational constant times planet mass (A.U.^/day^) 
Hollerith name of the planet 
Number of seconds per day (conversion factor) 

Execution error means of resolution, proportionality, 
and aiming angles Qt, B 


Noise constants for measurement types 

Planet number for launch, ephemeris, and target planets 
(nominally 4, for Earth) 

Number of km per A.U. (conversion factor) 

2ir (conversion factor) 

Dynamic noise constants (acceleration squared) 

Flag to choose Kallman-Schmidt filter (=1, nominal) 

Flag to use (=1) or not use (=0) generalized covariance 
capabilities (nominally 0) 

10-13 (constant) 

10+50 

(constant, "infinity") 

Flag to use dynamic noise constants 
Problem number 

Flag set if variation matrix (n) read in (nominally 0) 

Double precision constants 0., 1., 2., .5, 3. 

—1 9 — Q 

Double precision constants (lO”"*", lO”'^, . . . , 10 ^) 


Station location constants (altitudes in kilometers above 

mean radius of Earth, latitudes and longitudes in degrees) 
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BURN-A 


subroutine burn 

PURPOSE! TO CONFUTE THE CURRENT ACCELERATION DUE TO A FINITE BURN 


CALLING sequences CALL BURN(ACTH) 

ARGUMENTS! 

ACTH 0 ACCERATION VECTOR 

LOCAL SYMBOLS! 

DM3 

GRaV gravitational ACCELERATION 

amass initial SPACE CRAFT MASS 

SUBROUTINES REQUIRED! 

none 

COMMON USED: 

M TBURN 

T ALPHA 

SCMASS BETA 

THRMAG RPD 

XISP 

COMMON COMPUTED! 

OMASS COSB 

CURMAS SINB 

COSA 
SINA 
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SUBROUTINE SURNV 

PURPOSE! TO COMPUTE THE FINITE BURN PARTIALS 

CALLING sequence: CALL BURNV 

ARGUMENTS! 

none 

LOCAL SYMBOLS! 

none 

subroutines REQUIRED! 

none 

common USED! 

thrmag sinb 

SINA COSB 

COSA CURMAS 

COMMON COMPUTED! 


ACCPAR 



CALJUL-A 


subroutine caljul 

purpose; to compute julian date from calendar date or vice versa 
CALLING sequence; CALL CALUUL (DJ ♦ I Y *MO» ID* IH tMI *Si ICOOE ) 

ARGUMENTS; 


DJ 

JULIAN DATE 

lY 

CALENDAR YEAR 

MO 

MONTH OF YEAR 

ID 

DAY OF MONTH 

IH 

HOUR OF DAY 

MI 

MINUTE OF HOUR 

S 

SECONDS OF MINUTE 

ICODE 

I OPTION flag 


0=CONVERT FROM CALENDAR DATE TO JULIAN DATE 
laCONVERT FROM JULIAN DATE TO CALENDAR DATE 


SUBROUTINES REQUIRED: 

none 
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CAREL-A 


SUBROUTINE 6AREL 
PURPOSES 
CfiiLING SEGUEMGE8 
ARGUMENT 8 GN 

R(3J 

V«35 

TFP 


XI 




PP<3) 


CARTgSSAM COORDINATES TO CONIC ELEMENTS 
DALL G SREL SGMjRjyjTFPgAgEfWpKIpKKlpTfctiPPpQQpWMi 
2 GRAVITAf ZONAL C0MSTANT OF TME CENTRAL BODY 
I POSITION VECTOR RELATIVE TO CENTRAL BODY 

I VELOCITY VECTOR RELATIVE TO CENTRAL BODY 

0 Of FLIGMT FRON PERIAPSIS ON THE CONIC 

0 SEMl°MAJOR axis OF THE CONIC 

0 ECCENTRICITY OF THE CONIC 

0 ARGU50EMT OF PERIAPSSS OF THE CONIC 
0 IWeLINATIOM OF THE CONIC TO THE REFERENCE 


0 LONGITUOE OF THE ASCSNOICeS NODE OF' THE 


0 IS^STfiNTANEOUS ISSUE ANOMALY OF THE CONIC 
0 UeSIT VECTOR TOHASi PERIAPSIS ON CONIC 
0 U8sn VECTOR NORMAL TO PP IM ORBITAL PLANE 
0 UNIT VECTOR MORMAL TO ORBITAL PLANE 


SUBROUTINES REQUIRE08 
LOCAL SYMBOL SS AUXF 

fiva 

C®S0A 


NONE 

ECCENTRIC ANOMALY JMYPERBOLIC CASE) 

MEAM AMOHALY ^ELLIPTIC CASE) 

COSINE OF THE ECCEKITRIC ANOMALY ^ELLIPTIC 
CaSEH 


STA eOSIME OF TME TWi SMOMALT 

C WSGMITUOE OF THE aNGULAR MOMENTUM 

OIV INTERMEOIATE VARIABLE IN CALCULATION OF 

ECCE)3TRIC ANOMALY 
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CAR£L>A 


Eft 

P 

RftD 

RO 

RM 

SEMEft 

SSMHF 

STft 

TANG 

yM 

Z 


ECCENTRIC ANOHALf (ELLIPTIC CASE> 
SEHI-tftTUS RECTUM OF THE CONIC 

DEGREES TO RADIANS CONVERSION CONSTANT 

\ 

TIME DERIVATIVE OF RADIUS 

MAGNITUDE OF CARTESIAN POSITION VECTOR 

SINE OF THE ECCENTRIC ANOMALY (ELLIPTIC 
CASEJ 

HYPERBOLIC SINE OF AUXF 
SINE OF THE TRUE ANOMALY 

INTERMEDIATE VARIABLE USED TO CALCULATE 

SINHF 

MAGNITUDE OF THE CARTESIAN VELOCITY VECTOR 

INTERMEDIATE VECTOR USED TO CALCULATE 
PP, QQ VECTORS 
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GAR£L-1 




CAREL Analyeis 

CAREL Converts the cartesian state (position and velocity) of a massless 
point referenced to a gravitational body to the equivalent conic elements 
about that body. 

Lei the cartesian state be denoted r, v and let the gravitational con* 
etant of the central body be . 

The angular momentum constant c is 

c = I r X ^ I (1) 

The unit normal W to the orbital plane is 


W - tjLl 

c 


The semilatus rectum p Is 



The semi-major axis a Is 



( 2 ) 


(3) 


(4) 


Thus a ^0 for elliptical motion, a< 0 for hyperbolic motion. The 
eccentricity e is 


e “ 




1 - £ 
a 


(5) 


Thus e < 1 for elliptical motion, e > 1 for hyperbolic motion, 
inclination of the orbit i is computed from 


The 


cos 1 = W 


( 6 ) 


The longitude of the ascending node Q is defined by 


tan 0 


A 

w 

J 

/N 

-w 


(7) 
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CAREL-2 


The true anomaly f at the given state la computed from 


COB f “ P ~ ^ 

e r 


Bin f » 


cr 

MG 


( 8 ) 


Now define an auxiliary vector ? by 


/s 

z 


r -• 

— V 
c 


i ^ 
— r 
c 


(9) 


A A 

Then^ P , the unit vector to periapsls, and Q , the In-plane normal 
to P , are defined by 


A 

P 

= r cos 

f - 

z sin 

f 

(10) 

5 

.= r sin 

f + 

A 

z cos 

f 

(11) 


where r = — 
r 


The argument of periapsls w is then computed from 


tan (a) 



( 12 ) 


The conic time from periapsls t is computed from different formulae 
depending upon the sign of the simi-major axis. For a > 0 (elliptical 
motion) 


t 

P 



(E - e Bln E) 


cos E = e + COB £ E 

1 + e cos f 


^1 - e^ sin f 
1 + e cos f 


(13) 


For a 0 


(hyperbolic motion) the time from periapsls is 



(lA) 


Reference: Dattln, R. H. , Astronautlcal Guidance, McGraw-Hill Book Co., 

.New York, 1964. 


/ 
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CORREL-A 


SUOROUTl^g SORREL 

PURPOSES e@W^?ERT COVARiaMCE H&TRIX PARTITIOMS TO CORREtaiKON 
MfeT^SK Pi?iRTIYSOMS MD STSf^OSiRD DEVIATIOMS MQ WRITE 
THEM ©UT 


CALLING SEdUENCE® CALL CORREUPPcCKXSPijFSPsCXUPsUOsCKVP, V0#CXSUP. 

CHS¥PJ 


ARGUMENT 1 

PP 

I 

POSITION/VELOCITY COVARIANCE MATRIX 


CXXSP 

i 

CORRELATION BETMEEM SOLVE^FOR PARAMETERS 
AND POSITION/VELOCITY STATE 


PSP 

I 

SOLVE-FOR PARAMETER COVARIANCE MATRIX 


CXUP 

I 

CORRELATION 6ETMEEN POSITION/ VELOCITY STATE. 
AP3D DYNAMIC CONSIDER PARAMETERS 


u® 

I 

DYNAMIC CONSIDER PARAMETER COVARIANCE 
MATRIX 

' 

SKVP 

I 

CORRELATION BETWEEN POSITION/VELOCITY STATE 
mn MEASUREMENT CONSIDER PARAMETERS 

' 

VO 

I 

MEASUREMENT CONSIDER PARAMETER COVARIANCE 
MATRIX 


cxsup 

I 

CORRELATION BETWEEN SOLVE-FQR PARAMETERS 
AND DYNAMIC CONSIDER PARAMETERS 


cxsvp 

^ I 

CORRELATION BETWEEN SOLVE-FOR PARAMETERS 
AMD' MEASUREMENT CONSIDER PARAMETERS 


SUWOUTIMES SUPPORTED S 

prints SETEVIN GUIDf^ PREO 

local SVfSeOLS* ouw inverse of square root of .diagonal 

ELEWPNTS IM DVMANIC AND f^EASUREMENT 
eONSZDiER COVARIANCE PARTITIONS 

lEMD COUNTER INOICATSNG TOTAL NUMBER OF 
AUSKENTEO STATE VARIABLES 

ROM IKSTERMEOI ATE COMPUTATION AND OUTPUT VECTOR 


SQP INVEIRSE OF THE SQUARE ROOT OF OIEGOMAL 

elements in VEMICLE and SOtVE-FOR 
covaaiawcE partitions 


Z1 

lp[ IN 
XLA0 


STAKOARD OEVIATIOM 

SPRINT NDIMl NOIM2 NDIM3 ONE 

XSL XU XV 


COmON US£0!3 



COVMAT-A 


SUBROUTINE COV*^aT 


PURPOSE! TO CONVERT THE STANDARD DEVIATIONS AND CORRELATION INPUT 
TO COVARIANCES IN A SQUARE MATRIX 


CALLING sequence; CALL COVMaT ( P t NR » NACT UL ) 


arguments : 


P 

NR 

NACTUL 


array TO BE CONVERTED 
number Of rows USED 
ACTUAL DIMENSION OF MATRIX 


LOCAL SYMBOLS; K1 

JP 


SUBROUTINES REQUIRED 


INDEX 

INDEX 

: NONE 




COMMON USED! none 



COWELL-A 


subroutine COWELL 

PURPOSES TO CONTROL THE INTEGRATION LOGIC AFTER THE INTEGRATOR 
HAS BEEN INITIALIZED 

CALLING sequences CALL COWELL 5 TTO®XTOtVTO> 

ARGUMENTS: 

TTO I TIME {SEC) FROM EPOCH AT WHICH TRAJECTORY 
data is REQUIRE 

XTO 0 6 ELEMENT STATE VECTOR AT TIME TTO 

VTO 0 6 BY 20 ELEMENT MATRIX OF STATE PARTIALS 

LOCAL SYMBOLS: 

TINTEG TIME that HAS BEEN INTEGRATED TO 

STlME TIME aT INTERPOLATION 

STregT TIME ASSOCIATED WITH LAST ACCELERATION VECTOR 

IV FLAG FOR INTP TO INDICATE POSITION AND VELOCITY 

are REQUIRED 

ISW FLAG FROM TESTH TO INDICATE IF STEP SIZE MAS 

CHANGED 

lERR FLAG FROM CSTEP TO INDICATE THAT INTEGRATION 

STEP HAS FAILED 

SUBROUTINES REQUIRED; 

CSTEP 

testh 

INTP 

COMMON USED? 

? 

K 
H 

xoo 

NEQ 


XVDD SVl 

lOON SV2 

lELEVN 

SXl 

SX2 
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COWELL- 1 
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COWELL-2 


Interpolate acceleration arrays for state and state 
partials at requested time 
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COWELL- 3 


Increment point counter 


Have eleven points accumulated? 


Write last 11 acceleration points 
orbit file 


Check if step size should change 


Has stepsize decreased? 




COWELL-4 


Acceleration arrays still have space? 


Shift 20 old points out of acceleration arrays 


Reset acceleration index 


Skip writing an orbit file? 
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cstart-a 


subroutine CSTART (ORBINT entry POINT) 

PURPOSE: TO READ THE HEADER RECORD ON THE SEQUENTIAL ORBIT FILE 

CALLING sequences CALL CSTarT (NSEC « lERR ) 

ARGUMENTS: 

NSeC I DESIRED TRAJECTORY SECTION NUMBER 

lERR 0 ERROR FLAG 

=1 NORMAL RETURN 
=2 EOF detected 

=3, requested SECTION OUT OF RANGE 
=4 REQUESTED TIME OUT OF RANGE 

LOCAL SYMBOLS: 

IFRN LOGICAL FILE NUMBER 


COMMON computed: 

ymoic 

nstate 

xdd 

nsectn 

HMSIC 

kstate 

SXl 

NEQ 

AEINT 

IPART 

SX2 


SPINT 

GM 

xvod 


pvint 

T 

svi 


oblint 

H 

SV2 



FUNCTION DABSV 

PURPOSE: TO CALCULATE THE MAGNITUDE OF A VECTOR 

CALLING SEQUENCEI RES=DA8SV ( a t N ) 

arguments: 

A I INPUT VECTOR 

N I length of vector a 

DABSV 0 MAGNITUDE OF VECTOR A 

SUBROUTINES REQUIRED: 

NONE 
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SURPOUT:^i^f 


OAT A 




OATA-A 


PUPPOSt : TO Sh'T NECESSflPY IMUIAL VALUES AND DEFAULT VALUES FOR 

NAMFLIbT VARIARLES* TO READ INPUT DATA. TO TRANSLATE THESE 
INTO useable internal VALUF.Si TO COMPUTE DIMENSIONS OF 
state transition and covariance matrix PARTITIONS^TO ORDER 

measurement and event schedules, and to print the initial 

CONDITIONS. 

CALLING sequence: CALL DATA 


argumf:nts: 

NONE 



SURROUT INFS 

supported: 

STEAP- 

-F AMAIN 

subroutines 

required: 

COVMAT 

cstart 



SHIFT 

stvcrp 

LOCAL SYMBOLS: AR 

ARRAY 

USED To 


AMIN J 

minimum VALUE 


D 

TEMPORARY STOR 


DD 

TEMPORARY STOR 


DUM 

DUMMv 

VARI ABLF 


DUMl 

DUMMY 

VARI AULF 


GDATA 


PFCEO 

ZF.RMAT 


ECEU ROTATION matrix, ECLIPTIC TO EQUATORIAL 
KNOT JULIAN DATE OF FINAL TIME 

icnt counter used in scheduling 

IDAY DAY OF THE MONTH OF FINAL TIME 

lERPR FLAG TO PRINT NAMELIST *IF SET< 

lERR FLAG FOR FILE-READER INITIALIZATION SUCCESS 

IFCNRI flag SF T if control COVARIANCE INPUT 

IFGCOV FLAG SET IF CONTROL COVARIANCES INPUT IN 
FORM OF STANDARD DEVIATIONS AND 
CORRF.LAT IONS 

IFPCOV FLAG 'SET IF KNOWLEDGE COVARIANCES INPUT IN 
FORM OF standard DEVIATIONS AND 
■ CORRELATIONS 
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COMI^OM 


COMMON 


common 


DATA-B 


IFVMRI flag set if variation matrix ^araoa< input 
IHk hour of DAT OF FINAL TIME 

IMIN MINUTE OF HOUR OF FINAL TIME 
IMO MONTH OF year of final TIME 

lYR YEAH OF FINAL TIME 


IROW 

INDEX USED IN 

SCHEDULING 




LOAY 

DAY OF MONTH OF INITIAL T 

IME ON 

F 

ile 

LHP 

HOUR OF Day of 

INITIAL T1 

ME ON 

FI 

LE 

LMIN 

MINUTE OF HOUR 

OF initial 

TIME 

ON 

file 


LMO MONTH OF YEAH OF INITIAL TIME ON FILE 

LYR year of initial TIME ON FILE 

MEA5 TABLE OF MEASUREMENT TYPES PEUUESTED 

nDIMS number of columns OF AUGMENTED COVARIANCE 

nENT NUM8FH OF CAROS IN MEaSUkEMENT SCHEDULE INPUT 

SCriED ARRAY OF MEASUREMENT SCHEDULE TIMES 

SECI SECONDS OF MINUTE OF FINAL TIME 

SECL SECONDS OF MINUTE OF INITIAL TIME ON FILE 


computeo/used : OA.TFJ 

NOI M2 
PVINT 
XIG 


FNTM 

NOIM3 

TEV 

XSL 


lEVNT 
NO I m4 
TMN 
XU 


MCOOE 

NEV 

trtmb 

XV 


NOIMl 

NMN 

THTMl 


COMPUTFO : 

CPLU 

cxsu 

CXSUG 

cxsv 

CXSVG 

CXU 


CXUG 

CXV 

CXVG 

cxxs 

CXXSG 

EPS 


EuEC 

I AUG 

I AUGW 

lAUGDC 

1 AUGMC 

P 


PG 

PHI 

phig 

PPLU 

PS 

PSG 


HSPLU 

TG 

' UO 

VO 

XG 

XI 


USFU: OFLTM IftUGIN IGEN NEVl NEVi? NEV3 NEVA 

UNIVT SAL SLAT SLON 
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DAN(iMD-A 


FUNCTION DANGKiO 

PURPOSES TO MODULARIZE AN ANGLE ON TWO PI 
CALLING sequence: RES=DANGMD ( ANG) 

■v 

ARGUMENTSS 

ANG I THE ANGLE (RAO) TO BE MODULARIZED 

DANGMO 0 (ANG) MOD TWO PI 

subroutines required: 
none 


FUNCTION DANGV2 

PURPOSE: calculate A DIRECTED ANGLE BETWEEN TWO VECTORS IN 3 SPACE 

CALLING SEQUENCES RES=0ANGV2 ( A » B * REF ) 

ARGUMENTSS 

A I FIRST VECTOR 

8 I SECOND VECTOR 

REF I REFERENCE AXIS VECTOR 

0AN6V2 0 THE ANGLE FORMED BY ROTATING VECTOR A INTO 

VECTOR e CLOCKWISE* LOOKING IN THE DIRECTION 
OF THE vector ref 

SUBROUTINES PEQUIREDs 
DUXV 
oabsv 

DOOT 


subroutine DaVECT 

PURPOSES VECTOR ADDITION 

CALLING sequences CALL DA VECT ( A *8 *N* C ) 

ARGUMENTSS 

A I FIRST vector 

S I SECOND VECTOR ' 

N I length OF VECTORS A AND B 

C 0 VEC(A) ♦ VEC(B) 

SUBROUTINES REQUIREDs 
NONE 
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OOOT-A 


FUNCTION DDOT 


PURPOSE; WEcTOR dot PRODUCT 
CALLING sequence; R£S=D00T ( A cB »N) 


ARGUMENTS: 

A 

I 

FIRST VECTOR 

B 

I 

SECOND vector 

N 

I 

LENGTH OF VECTORS A AND B 

DDOT 

0 

VEC(A) DOT VEC(S) 


subroutines required: 
none 


function dootb 

PURPOSE: to return THE DOT (INNER) PRODUCT OF Tki/0 3-VECT0RS 

CALLING sequence; ADOTB = OdOTB(A»B) 

arguments; a first vector 
B second vector 

subroutines required; none 


SUBROUTINE DMADD 


PURPOSE; MATRIX ADDITION 


CALLING sequence; CALL OMADD ( A ♦ B ♦ C * I » J ) 
ARGUMENTS; 

A I FIRST MATRIX 

B I SECOND matrix 

C 0 RESULT matrix = (A) ♦ (B) 

I I NUMBER OF ROWS IN MATRIX 

U I NUMBER OF COLUMNS IN MATRIX 


SUBROUTINES REQUIRED: 

none 


« 
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DMADD 


>• 


subroutine OMADO 
PURPOSE: MATRIX ADDITION 


CALLING sequence? CALL DMADD ( A 9 0 ® C s I ? J ) 


arguments: 


B 

c 

I 

J 


I FIRST MATRIX 
I SECOND matrix 

0 result matrix = (AJ «B» 

1 NUMBER OF ROWS IN MATRIX 

I NUMBER OF COLUMNS IN MATRIX 


subroutines REQUIRED: ^ 
NONE 


SUBROUTINE DMATPY 

PURPOSE: MATRIX MULTIPLICATION 


CALLING sequence: CALL DMATPV ( A ® 0 »C »N »M » IP ) 


ARGUMENTS; 

A I FIRST MATRIX 

S I SECOND matrix 

c 0 result Matrix = ja) * ibs 

N I number of rows IN A 

M I number of columns IN A« NUMBER OF ROWS IN 8 

IP I number of columns in b 

subroutines required: 

NONE 


subroutine OmSUB (DMADO entry POINT) 

PURPOSE: MATRIX SUBTRACTION 

CALLING sequence: CALL DMSUB ( A « B »C • I » J) 

arguments: 

A I FIRST matrix 

B I SECOND matrix 

c 0 result Matrix = u> - (B) 

I I NUMBER OF ROWS IN MATRIX 

vi I number Of columns in MATRIX 

subroutines REQUIRED: 

NONE 


DRD 


subroutine dro 

PURPOSE? TO COMPUTE LATITUDE AND LONGITUDE OF A VECTOR 

CALLING sequences CALL DRD {SsRAkOECLJ 

ARGUMENTS? 

S 1 INPUT vector 

HA 0 RIGHT ascension OR LONGITUDE OF S VECTOR 

DECL 0 DECINATION OR LATITUDE OF S VECTOR 

subroutines REQUIRED! 

NONE 


subroutine dsmift 


PURPOSE! to shift ONE VECTOR INTO ANOTHER 
CALLING SEQUENCE! CALL DSHIFT (A»N»B( 


ARGUMENTS! 

A I ADDRESS OF FIRST DOUBLE WORD TO BE SHIFTED 

N I NUMBER OF DOUBLE WORDS TO BE SHIFTED 

B 0 ADDRESS OF FIRST DOUBLE WORD IN RECEIVING VECTOR 


SUBROUTINES REQUIRED! 

none 


subroutine DSVECT IOAVECT ENTRY POINT) 

PURPOSE! VECTOR SUBTRACTION 

CALLING SEQUENCE! CALL DSVECT ( A »B . N » C ) 

ARGUMENTS! 

A I FIRST VECTOR 

B I SECOND VECTOR 

N I length of VECTORS A AND 8, 

C 0 VECIA) - VEC(B) 

subroutines required: 
none 


DUNIT 


SUBROUTINE DUNIT 


PURPOSE; TO UNITIZE A VECTOR 


CALLING .sequence; CALL DUNIT 


ARGUMENTS; 

A 

I 

INPUT vector 

X 

N 

I 

LENGTH OF VECTOR 


B 

0 

UNIT VECTOR IN THE DIRECTION OF A 



subroutines required; 
dabsv 


subroutine duxv 

PURPOSE; TO FORM A VECTOR CROSS PRODUCT 
CALLING sequence; CALL UXVU.BiC) 
ARGUMENTS; 


A I FIRST vector 

B I SECOND VECTOR 

C 0 result vector = (A) CROSS (B) 


SUBROUTINES REQUIRED; 

none 


subroutine DvCOMB 

PURPOSE: TO COMBINE (AOOf TWO VECTORS* EACH MULTIPLIED BY SCALARS 

calling sequence; call nVCOMB{A*SA»B*SB»C) 

ARGUMENTS; 

A I FIRST VECTOR 

SA I SCALAR WHICH IS APPLIED TO VECTOR A 

B I SECOND VECTOR 

SB I SCALAR WHICH IS APPLIED TO VECTOR B 

C 0 RESULT VECTOR = ♦ S0*(B) 

SUBROUTINES REQUIRED; 

none 
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SUBROUTINE OVECRD 

PURPOSE; TO COMPUTE A UNIT VECTOR FROM ITS RIGHT ASCENSION AND 

declination 

CALLING sequence: CALL DVECRD ( R A . DECL » S ) 

arguments; 

HA I RIGHT ASCENSION OR LONGITUDE (HAD) 

OECL I DECLINATION OR LATITUDE (HAD) 

S 0 OUTPUT VECOTR 

SUBROUTINES REQUIRED; NONE 


FUNCTION OVMaG 

PURPOSE: calculate THE MAGNITUDE OF A 3-VECTOR 

CALLING SEUUENCE: AMAG = DVMAG(A) 

ARGUMENTS; A VECTOR 
subroutines REQUIRED; NONE 


subroutine DVSDIV (SAVECT entry POINT) 

PURPOSE: VECTOR SCALAR DIVISION 

CALLING sequence; CALL DVSOIV ( A»SA»N«C) 

arguments; 

A I INPUT VECTOR 

SA I SCALAR BY ii^HICH VECTOR A IS DIVIDED 

N I length of VECTOR A 

C 0 VEC(A) / SA 

subroutines required: none 


subroutine DVSMLT (SaVECT entry POINT) 
PURPOSE; VECTOR SCALAR MULTIPLICATION 
CALLING sequence: CALL DVSMLT ( A » SA » N» C ) 


ARGUMENTS 



A 

SA 

N 

C 


1 INPUT vector 

I SCALAR BY WHICH VECTOR A IS MULTIPLIED 
I LENGTH OF VECTOR A 
0 SA * VEC(A) 


subroutines REQUIRED: NONE 
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SUBROUTINE OVSTAT 

PURPOSE: compute AND PRINT STATISTICAL PARAMETERS OF A TRIM MANEUVER 


CALLING sequence: CALL OVST AT CS »E ♦ V ♦ T * TR 9 RHO i SMX ) 

I 


ARGUMENTS; 


S INPUT MATRIXc GAMMA*P*GAMMA (TRANSPOSE) ' 

E ARRAY OF EIGENVALUES OF S 

V array OF EIGENVECTORS CORRESPONDING TO E 

T ANNIHILATION LIMIT ( IF ( E (I ) oLT , T ) E (I ) =ZERO ) 

TP trace of S 

RHO MEAN VALUE OF DELTA-V 

SMX MAXIMUM EIGENVALUE IN E 


LOCAL SYMBOLS: 


D Table of percentile levels and oelta-v magnitudes 

DV ARRAY OF CONSTANTS USED TO COMPUTE MAGNITUDES 

IK * 

IKl *«* INDICES USED IN LAGRANGIAN 6-POINT INTER- 
IL * POLATION 

ILl * 

K2 Ratio* smd/smx 

L2 RATIO* SMN/SMX 

SO STANDARD DEVIATION OF DELTA-V VALUES 
SMD MIDDLE EIGENVALUE IN E 
SMN MINIMUM EIGENVALUE IN E 
STR SQUARE ROOT OF TR 
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DVSTAT Analysis 


The Lee-Boain analytic solution for V statistics involves a hypergometric 
function, and is described in detail in Reference 6 . ,For this subroutine, 
a table DV has been generated from their solution which is used to cal- 
culate the mean and standard deviation and values for the 90, 99, 99»9 
and 99.99 percentile levels in a much less time. The subroutine JAC0BI 
obtains the eigen values and eigenvectors of the input S matrix. The 
ratios of the middle and smallest eigenvalues to the largest eigenvalue, 

2 2 

k and 1 respectively, are determined and used in a Lagrangian 6 -point 
interpolation from table DV. The six points are located as follows 


. P . p 

4 5 


. , (iwidi [Mil) [xi - 

1 \10 ’ 10 r 


greatest integer less than or 


where P 
equal to X 

and the other points are at appropriate .1 intervals. 
Define 


p = [lOk^ + 1 ] - lOk^ 
q = [10^^ + 1 ] - loi^ 

Use Lagrangian 6 point interpolation: 

f(k^, P) = pq, f(p^) + q(q- 2 p+l)- f (P 3 ) / 2 

+ p(p- 2 q+l)' f (P 3 ) / 2 
+(l+pq-p^-q 2 ). f (p^) 

+p(p-l)* f(p^) / 2 
+q(q-l)- f(Pg) /2 

2 

This necessitates data points at .1 intervals over the range 04k ^ 1 . 1 , 

0^J?<k^, 
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Points outside this range, even though they might figure in the inter- 
polation theoretically, are unnecessary because either p or q or both 
become 1, and the corresponding terms drop out. 

After values have been so generated for Av Av Av 999, Av 

the values are multiplied by the square root of the trace, which is 
normalized to one in generating the data points. The standard deviation 
is calculated from the classic formula, tr= (tr-A'^) 2 , 
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S'JHP.UUTINt DYNO 

PUPPOSes COMPUTE ASSUMED ANU ACTUAL COVARIANCE 

rRIX IN THE ERROR ANALYSIS PROGRAM 


CALLING SEQUENCE^ CALL OYNOCICODEI 
AWGUWENT8 ICOOE 


I =0 ASSUMED 

=1 actual 


SU;^»OUTINES SUPPORTEDl ERRANN SETEVN 6UI0M P«EO 
local SYMBOLS 8 02 SQUARE OF {DELTM*TH) 

COMMON C0MPUT608 Q 

COMMON USED8 OELTM ONCM IDNF TM 

TrnKiF r.QHCN 


DYN0 Analysis 

SubrouLinu UYN 0 evaLuatuH Chu assumed dynamic cuvarlance matrix Q 
over the time interval t ■ t, , - t, if LCOOli “ 0 . If LC 0 DL * 1 

K+1 k. 

the actual dynamic noise covariance matrix Q' Is evaluated over the 
same Interval. In either case the dynamic noise covariance matrix 
Is assumed to have the form 

Q = dlag (Ja K] Ac*', k K2 At*', h K3 At**, Kj K2 A , K3 At^) 

where dynamic noise constants Ki , K2, and K3 have units of km^/a*'. 
To compute the actual dynamic noise covariance matrix Q', we simply 
replace K2, and K3 with the actual dynamic noise constants 
Kj , K2, and Kj, respectively. 
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SUBROUTINE DX8 

PURPOSE! calculate THE CROSS PRODUCT OE TWO 3-VECTORS 

CALLING SEUUENCE: CALL OXB ( A A 9 BB « VECPR ) 

ARGUMENTS: AA FIRST INPUT VECTOR 

RR SECOND INPUT VECTOR 

VECPR VECTOR PRODUCT 

SUBROUTINES REQUIRED: NONE 

SUBROUTINE DZERO 

I 

purpose: TO GENERATE A ZERO VECTOR 

CALLING sequence: CALL OZERO(A»N) 

I 

ARGUMENTS: 

A I address of first DOUBLE «IORD TO BE ZEROED 

N N number of double WORDS TO BE ZEROED 

SUBROUTINES REQUIRED: ^ 

none 



ECOMP-A 


SUBROUTINE ECOMP 

PURPOSE! TO COMPUTE DIFFERENTIAL TRANSFORMATION RELATING TARGET 

variables to state 

calling sequence? call ECOMP JSReAAsBB) 

ARGUMENTS 

XK I CURRENT STATE VECTOR 

AA 0 PARTIALS OF TARGETS WITH RESPECT TO POSITION 

B0 0 PARTIALS OF TARGETS WITH RESPECT TO VELOCITY 

LOCAL symbols? 


XPERT 

TEMPORARY PERTURBATION VALUE 



7FP 

TIME from PERIGEE 



El 

SEME-MAJOR AXIS 



E2 

ECCENTRICITY 



ATN 

3 VECTOR OF TARGETS CORRESPONDING 
PERTURBATION MADE TO XX 

TO 

NEGATIVE 

Alp 

3 VECTOR OF TARGETS CORRESPONDING 
PERTURBATION MADE TO XX 

TO 

POSITIVE 


subroutines REQUIRED! 

carel 


COMMON USE!>! 

PERT 

TM 

SMU il ^ 5 
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SUBROUTINE EIGHT 

PURPOSES TO CONTROL THE COMPUTATION OF EIGENVALUES, EIGENVECTORS, 
ANO MYPERELLIPSOSOS® 

CALLING SEaUE^CEs CALL EIGHT^VEIG ,f OK,HARG,IFMT) 

ARGUMENTS ^EZG I MATRIX TO BE QIAGONALIZED 

FOX S FINAL OFF-DIAGONAL ANNIHILATION VALUE 

HARG I MATRIX FOR WHICH THE HYPERELLIPSOID IS TO 
BE COMPUTED 

IFMT I FORMAT FLAG 

si, PRINT POSITION EIGENVALUE TITLE 

=2, PRINT VELOCITY EIGENVALUE TITLE 

=3, PRINT EIGENVALUE TITLE 


SUBROUTINES SUPPORTED! ATCEGV MOMENT PRED ' GENGID GUIOM 

SUBROUTINES REQUIREOI JACOBI 

LOCAL SYM80LS8 EGVCT EIGENVECTOR MATRIX 

I 

EGVL EIGENVALUE MATRIX 

OUT SQUARE ROOTS OF EIGENVALUES 

EIGHT Flow Chart 
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EPHGT 


SUBROUTINE EPHGT 

PURPOSE! TO RETRIEVE FROM THE DIRECT ACCESS SLP FILE* THE STATE 

VECTOR OF A PLANET WITH RESPECT TO THE SUN AT AN ARBITRARY 

JULIAN DATE 

CALLING sequence; CALL EPHGT aP»DJ»R»V ) 

ARGUMENTS? 

IP I PLANET NUMBER (1=5UN* 2=MERCURY» ETC) 

DJ I JULIAN DATE 

R 0 RADIUS VECTOR FROM SUN TO PLANET IP 

V 0 VELOCITY VECTOR OF PLANET IP WITH RESPECT To THE SUN 

LOCAL SYMBOLS! 

IPGSFC VECTOR CORRELATING NOMNAL PLANET NUMBERING SYSTEM 
with GTDS PLANET NUMBERING CONVENTION 
array TEMPORARY TRANSMISSION ARRAY BETWEEN EPHGT AND 

subroutine EVAL 

IARRAY TEMPORARY TRANSMISSION ARRAY BETWEEN EPHGT AND 
subroutine EVAL 

COMMON USEI>; 

ymdic 

HMSIC 
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EPHGT Analysis 

Subroutine EPHGT is used by ERRAN and NOMNAL to retrieve heliocentric- 
ecliptic planetary state vectors from the solar/lunar/planetary direct 
access epheraeris file. This is accomplished by saving and resetting the 
year, month, day and hours, minute, seconds variables (for initial con- 
ditions) and a dummy gravitating bodies vector with the sun as the 
central body and the planet of interest as the only non-central body 
is then generated. Subroutine EVAL is then called to generate these 
vectors , 
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PROGRAH £RR6W 

PURPOSES TO COf^TROL THE COMPUTATIONAL FLOW THROUGH THE BASIC 
CYCLE (MEASUREMENT PROCESSING) AND ALL EVENTS IN THE 
ERROR ANALYSIS MODE. 

subroutines supported* error 

SUBROUTINES REQUIRED® SCHED NTM PSIM DYNO TRAKM 

MEMO GNAVM PRINTS SETtVN GUIOM 

MEAN GPRINT 

LOCAL SYMBOLS® 


ICOOE EVENT CODE 

IPRN MEASUREMENT COUNTER FOR PRINTING 



NEVENT EVENT COUNTER 


TRTM2 TIME OF THE MEASUREMENT 


COMMON 

computed/useos 

ICOOE 

XF 

MCNTR 

XI 

RI 

TEVN 

TRTMl 

COMMON 

COMPUTED® 

OELTM 





COMMON 

USED* 

FNTM 

NMN 

lEVNT 

NR 

IPRINT 

NTMC 

ISTMC 

RF 

NEV 

TEV 


kprint 
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ERRAN Analysis 


Subroutine liRRAN controls the computational flow through the basic 
cycle (measurement processing) and all events In the error analysis/ 
generalized covariance analysis program. 


In the basic cycle the first task of ERRAN is to_control the gen- 
eration of the targeted nominal spacecraft state X, at time 


^k+l' S^ven the state Xj^ at time tj^. Then calling PSIM, DYNjJ, 

TRAKM, and MEN?), successively, ERRAN controls the computation of 
all matrix information required by subroutine GNAVM to compute the 
actual and assumed knowledge covariance matrix partitions at time 



immediately 


following the measurement. 


At an event, ERRAN simply calls the proper event subroutine or 
overlay where all required computations are performed. 
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ERRAN-3 


Call MEN0 to compute assumed 


Call GNAVM to compute assumed covariance 
matrix partitions at and tj^^^ 


Increnrent print counter IPRN 

" ' ' i — ~ 

Is it time to print? 


IGENj 0? 

Call DYN0 to compute actual Q 


k+l,k 


Call MEN0 to compute actual 


Call GNAVM to compute actual 
2nd-moment ratrix partitions 

n ard 


Call MEAN to compute actual 
estimation error means at 












4-47 


ERRAN-4 






ERRAN-5 
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SUBROUTINE EULMX 

PURPOSES TO COMPUTE THE MATRIX REQUIRED TO DEFINE TRANSFORMATIONS 

FROM ONE COORDINATE SVSTEM TO ANOTHER. 

CALLING SEQUENCES CALL EULMX «ALP«NN, BET >MMfGAN»LL * P ) 

ARGUMENTS ALP I FIRST ROTATION ANGLE SRAOIANS) 

MM I FIRST AXIS OF ROTATION 

BET I SECOND ROTATION ANGLE (RADIANS) 

MM I SECOND AXIS OF ROTATION 

GAM I THIRD ROTATION ANGLE (RADIANS) 

U I THIRD AXIS OF ROTATION 

P(3,3) 0 TRANSFORMATION MATRIX 


SUBROUTINES REQUIREDS NONE 

LOCAL SYMBOLS S A INTERMEDIATE ROTATION MATRIX 


ALPHA temporary LOCATION FOR EACH OF THE 
ROTATION ANGLESS ALP» BET. AND GAH 

D INTERMEDIATE PRODUCT MATRIX 


F TRANSFORMATION MATRIX FOR ANGLE ALP 

G TRANSFORMATION MATRIX FOR ANGLE BET 

H TRANSFORMATION MATRIX FOR ANGLE CAM 

N , COUNTER SHORING NUMBER OF COORDINATE AXES 

FOR WHICH CALCULATIONS REMAIN 


MAXIS TEMPORARY LOCATION FOR EACH OF THE AXES 
OF ROTATIONS NN,MM, AND LL 





GAINl-A 


SUBROUTINE CftINl 

PURPOSES TO COMPUTE THE KALMSN GftlM MATRICES 

CALLINC, SEQUENCES CALL GAI N1 (NR , A J s AKH » S M , IE NO J 

ARGUMENTS8 NR 1 NUMBER OF ROWS IN THE OBSERVATION NATRIK 

AJ I MEASUREMENT RESIDUAL COVARIANCE ANO ITS 

INVERSE 


AKH I 

INTERMEDIATE 

ARRAY 

SW I 

INTERMEDIATE 

ARRAY 

I END I 

NR-1 


SUBROUTINES SUPPORTED 1 

GNAVM 


SUBROUTINES REQUIRE08 

MATIN 


LOCAL SYMBOLS? OUM 

INTERMEDIATE 

VECTOR 

X J 

INTERMEDIATE 

ARRAY 

SUM 

INTERMEDIATE 

VARIABLE 

COMMON C0MPUTE08 

AK S 


COMMON USEOI 

ONE MALE 

ZERO 
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6AIN1 Analysis 

Subroutine GAINI computes the Kalman-Schmldt filter gain matrices 

and that are used in subroutines GNAVM and NAVM to update 

estimation error covariance matrices after a measureir.ent has been 
processed. 


The 

matr 


measurement residual covariance matrix J, , and the auxiliary 

k+1 

ices A^^j^ «nd are assumed to be available (from GNAVM or 


NAVM) when GAINI is called. Subroutine GAINI then evaluates the 
following equations to determine the filter gain matrices: 


\+l " \+l “^k+l 


( 1 ) 



'k+1 


k+1 



( 2 ) 



/ 
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SUBROUTINE GDATfl 

PURPOSE* TO INITIALIZE GENERALIZED COVARIANCE QUANTITIES 


CALLING SEQUENCES CALL GDATA 


SUBROUTINES SUPPORTEOI 

DATA 





COMMON COMPUTED/USEOI 

EU 

EV 

EVA 

EVB 

EVK 


EVS 

EH 

exi 

EXSI 

exsT 


EXT 

GCUV 

GCUW 

GCVH 

GCXSU 


GCXSUG 

GCXSV 

GCXSVG 

GCXSH 

GCXSWG 


GCXU 

GCXUG 

GCXV 

GCXVG 

GCXH 


GCXWC 

GCXXS 

GCXXSG 

gqncn 

GMNCN 


GP 

GPG 

GPS 

GPSG 

GU 


GV 

GH 

lONF 

VARA 

VARO 


VARK 

VARS 




COMMON USEOI 

CXSU 

CXSV 

CXU 

CXV 

CXXS 


DNCN 

lONF 

MNCN 

NOIHl 

NDIM2 


NDIM3 

N0IM4 

P 

PS 

SIGALP 


SIGBET 

SIGPRO 

SIGRES 

TG 

UO 


\J Q 

ZERO 
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SU8R0UT IN£ 


GEMGIO 


PURPOSES FO GENERATE THE ENSEMBLE STATISTICS OF THE ACTUAL 

COMMANUED VELOCITY CORRECTIONS THE ACTUAL EXECUTION 
ERROR AND THE ACTUAL TARGET MISS 

CALLING SEQUENCES CALL 

GENGIO 

SUBROUTINES SUPPORTED* 

ERRAN 

SUBROUTINES REQUIREU* 

SAVMAT OYNO GNAVM MEAN MOMENT 

EIGHT GQCOMP ATCEGV JACOBI DVSlAT 

local symbols* AMAX 

INTERMEDIATE VARIABLE 

ATC 

actual target CONDITION 2ND MOMENT MATRIX 

B 

INTERMEDIATE VARIABLE 

BBB3 

BLANK label ARRAY 

C 

INTERMEDIATE VARIABLE 

OELTM 

TIME DIFFERENCE 

EBOVB 

MAGNITUDE OF ACTUAL STATISTICAL OELTA-V 

EOVN 

MEAN OF ACTUAL COMMANDED VELOCITY 

correction 

EGM 

MAGNITUDE OF EIGENVECTOR CORRESPONDING TO 
MAXIMUM EIGENVALUE 

1 

EGVCT 

EIGENVECTOR ARRAY 

EGVL 

eigenvalue vector 

elab 

LABEL 

EX IS 

STORAGE FOR EXI 

EXSIS 

STORAGE FOR EXSI 

EXTS 

STORAGE FOR EXT 


ACTUAL STATISTICAL OELTA-V 

; GAP 

ACTUAL VELOCITY CORRECTION 2ND MOMENT 

MATRIX 

GPSAVE 

STORAGE FOR GP 
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GTG TIME OF ACTUAL GUIDANCE EVENT 

IFLAG SI BEFORE GUIDANCE EVENT 

s2 AFTER GUIDANCE EVENT 

III INDEX DEPENDING ON GUIDANCE EVENT TYPE 

HAP INDEX OF HAXXHUH EIGENVALUE 

PEIG INTERMEDIATE ARRAY 

Q actual execution error 2NO MOMENT MATRIX 

ROM INTERMEDIATE VECTOR 

S INTERMEDIATE ARRAY 

SUM INTERMEDIATE VARIABLE 

U actual commanded VELOCITY CORRECTION 

VEIG INTERMEDIATE VECTOR 

ZLAB LABEL 

ZV ACTUAL EXECUTION ERROR MEANS 


ZZ INTERMEDIATE VARIABLE 


COMMON 

C0KPUTE0/USE08 

DUMMY Q 

EXl 

EXMEAN 

EXSI 

EXT 


' 

GCXSUG 

GCXSVG 

GCXSVG 

GCXUG 

GCXVG 



GCXHG 

XLAB 

GCXXSG 

GP 

GPG 

GPSG 

COMMON 

USEOI 

ADA 

DVUP 

EE 

EEE 

EU 



EV 

Em 

EXST 

FOP 

FOV 



GA 

GCXSU 

GCXSV 

GGXSR 

GCXU 



GCXV 

GCXN 

GGXXS 

GPS 

GU 



GV 

GM. 

IGP 

I GUI a 

11 



NOIMi 

NDIMZ 

NDIMS 

NOIMi^ 

PI 



OPR 

RPR 

TEVN 

TG 

T INJ 



XIC 

XSL 

XU 

HV 
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( 


GENGID Analysis 

Subroutine GENGIU controls the execution of generalized guidance 
events. Generalized guidance has been extended to all guidance 
options defined for subroutine GUIUM except for final insertion. 


Unlike GUIDM, which computes target dispersions and fuel budgets 
based on filter-generated statistics, subroutine GENGID computes 
target dispersions and fuel budgets based on actual statistics. 
In other words, the generalized covariance technique as applied 
to the guidance process is programmed in GENGID. The required 
equations are summarized below. 


Before the guidance event at time t^ can be executed, it Is neces- 
sary to propagate the actual control mean and control 2nd-moment 
matrix partitions forward to tj from the previous guidance event 

at time t 


J-1‘ 


The control mean propagates according to 


Y.', " , + 0 

1 J-1 XX 


8 


c + u" + 8 w' 
S . XU o xw o 
0 


( 1 ) 


where , 8 


XX 


s 


XU 


and 


8 are stats transition matrix partitions 
xw 


over the Interval 


fj-l» ^j]* 


and X 


* 8 * 


u , and w denote actual 


position/velocity and solve-for, dynamic-conslder, and Ignore param- 
eter deviation means. The notation ( ) ’' Indicates actual values 
as opposed to the unprimed assumed values, while ( ) and ( )'^ 
Indicate values Immediately before and f,fter the execution of the 
guidance event, respectively. The actual control position/velocity 
2nd'moment matrix is defined by 


- H [x- . 

j 


( 2 ) 


The remaining control 2nd-moment matrix partitions are defined 
similarly. Th® propagation equations appearing in subroutine 
GNAVM are used to propagate the control 2nd-moment matrix parti- 
tions over the interval pj t”|. 
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The actual target state deviation la related to tlie actual 


state deviation Xj at time tj according to 


“i ■ “j 


(3) 


where Hj Is the variation matrix for the appropriate midcourse 
guidance policy. The mean of 6 t' is given by 


E [«,'] . nj E [.'] 


(4) 


The statistical target dispersions are represented by the actual 

(5) 


target condition 2nd-moment matrix W', which is defined as 


W' = E 6Tj . 


Substitution of equation (3) into equation (5) yields 


( 6 ) 


Equations (4) and (6) are evaluated Immediately before and after 
the guidance correction to detemine how much the target errors 
have actually been reduced by the velocity correction at t^ . 

The actual commanded velocity correction 2nd-moment matrix is de- 
fined by 


Sj -E 


AVj AV'- 


<7) 


where the actual commanded velocity correction is glvan by 

The guidance matrix corresponds to the appropriate linear mid- 
course guidance policy. The equation used to evaluate S" is given 


by 


J 


■j-'i 




(9) 


where s is a scalar input by the analyst, generally 0. s ^ 1., 


and 
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where all E 
of equation 



terms have been neglected in the derivation 


The mean of 
by applying 


the actual commanded velocity correction Is obtained 
the expectation operator to equation (8) : 


h'] ■ ' j I" [-j'] ^ " [“)]|- 


< 10 ) 


Since this equation gives no useful information 
studies, the Hoffman-Young formula will be used 


USIX CO « 


for fuel-sizing 
to evaluate 


E 




where 


A = 


^race 



B = A' X' + >^3 + ^2 


( 11 ) 


and X', and X' axe the eigenvalues of the 2nd-moment matrix 

Sr* If this mean is vanishingly small, the Lee-Boain analysis is used 
to obtain the statistical parameters, including the effective ^V. 
Otherwise the actual effective or statistical is defined as 



( 12 ) 


where a 


J 


denotes a unit vector In the most likely direction of 


the velocity correction. The most likely direction Is assumed 
to be aligned with the eigenvector associated with the maximum 
eigenvalue of ST. 


With 


"E 


available, the actual execution error statistics 


can be computed (by calling subroutine GQCf^MP). Those are the 
actual execution error mean E 2nd-moment maxtrix 


Qj defined as 


L 


" E 




5AV 


JJ 


(13) 
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It remains to summarize the equations which are used to update 
all actual control and knowledge means and 2nd-momeut matrix par- 
titions immediately following the execution of a guidance event. 
The actual estimation error means and 2nd-moment matrix partitions 
are updated using the following equations: 



i T 

where A = [0 i I] . The actual deviation means are updated using 
the following equations: 


E 

E 


4* 


“ 


= - E 





-h' 


'Vf-" 

X 

= - E 

L'jJ 


(18) 


(19)' 


The entire set of actual control 2nd-moment matrix partitions is 
updated by equating them to the corresponding actual knowledge 

2nd-moment matrix partitions at tT. 
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GENGID Flow Chart 
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4-6G 




I 


GENGID-7 


[AvpU' 0. 


Use Hoffman-Young 
empirical' eqns. to 
obtain statistical 
parameters and size of 
actual statistical 
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Compute actual target dispersions 

EUtp and W^, immediately afterj 

executing the desired midcourse 
guidance event. Write out, along 
with eigenvalues and eigenvectors. 


t 


All actual statistics have been 
updated immediately following the 
execution of the guidance event. 



RETURN 
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SUBROUTINE GETCOW (ORBINT ENTRY POINT) 

PURPOSES TO generate A STATE VECTOR AT A REQUESTED TIME BY 
interpolating DATA ON THE SEQUENTIAL ORBIT FILE 

CALLING sequence; CALL GETCOR ( NSEC » TREQ ♦ 1 ERR ♦ X ♦ STM ) 


ARGUMENTS; 

nsec 

I 

number of desired trajectory section 

TREQ 

I 

REQUEST TIME OF DATA (SEC) 

lERR 

0 

ERROR FLAG 



=1 NORMAL RETURN 
=2 EOF DETECTED 

=3 REQUESTED SECTION OUT OF RANGE 
=4 REQUESTED TIME OUT OF RANGE 

A 

0 

STATE VECTOR 

stm 

0 

STATE PARTIALS 

local symbols: 

ifrn 


LOGICAL FILE NUMBER 


SUBROUTINEJi REQUIRED: 
INTP 


COMMON USED/COMPUTED: 
NEO 


COMMON computed: 


T 

XVDO 

H 

SVl 

XDD 

SV2 

SXl 

NSECTN 

5X2 



« 
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SU8R0UT£ME 

PURPOSES 70 eOl^PyTE THE GPEEMWZCH MQLE AND THE UNIVERSAL 

7If1E iM QAVS» NMICH IS USED IN THE TRACKING HOOULE TO 
ORIENT THE TRACKING STATIONS ON A SPERICAL ROTATING 
EARTH« 

CALLING SEQUENCES CALL GHA 
ARGUMENTS 8 ^ONE 

SU8R0UTINES SUPPORTED 8 OATAIS OATAl 

NUMBER OF OATS IN TSTAR 
EARTH ROTATION RATE 
GREENWICH HOUR ANGLE 
INTERMEDIATE VARIABLE 
JULIAN DATE OF JAN. 0, 1950 
FRACTION OF DAY IN TSTAR 

/ 

JULIAN DATE, EPOCH JAN. 0, 1950» OF 
INITIAL TRAJECTORY TIME 

UMIVT 

CA7EJ EM13 


LOCAL SYMBOLS 8 0 

EQMEG 

6H 

ID 

REF JO 
TFRAC 
TSTAR 

» 

f 

COMMON CONPUTEOe 
COMMON USfg)S 
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Subroutine GHA computes the Greenwich hour .angle In degrees and days at 
some epoch T* referenced to 1950 January 1 0 . £poch T* is computed from 

T* - J.D.q + 2415020,0 - 

where 

J.D. “ Julian date at launch time t referenced to 1900 
January 0 12 . ° 

“ Reference Julian date 2433282.5 

KXLr 

” 1950 January I'^O*' referenced to January 0^12^ of 
the year 4713 B.C. 

j u H h 

and 2415020i0 ■ 1900 January 0 12 referenced to January 0 12 

of the year 4713 B.C. 

Then ia the Julian data at launch time C referenced to 1950 January 

ld0\ ® 

The Greenwich hour angle correapondlng to T* is given by ' 

GHA(T*) - 100.0755426 + 0,985647346d + 2.9015 x lO"^^ d^ +ut 

where 0 ^ GHA(T*) < 360? 

and d » integer part of T*, t “ fractional part of T*, 

and Earth's rotation rate is degrees/day. 

GHA 

The Greenwich hour angle la days ia given by 
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GHA Flow chsrt 


ENTER 



YES / * V NO 

I C GH 4 0 7 ) 


GH GH + 360. 


yes y ^ NO 

-C GH - 360. ^ 0 ? > 


GH “ GH - 


Compute Greeawich hour angle 
in days. 


RETURN 
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GIDANS-A 


SUBROUTINE GIOANS (PRELIM ENTRY POINT) 
purpose: dummy link with non halo ORBIT OPTIONS 

CALLING sequence: CALL GIDANS 
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SUTRQUTINE CNAyM 


PURPOSEt 10 PROPA&aTE ASSUMED COVARIAnCE MATRIX PARTITIO^S P, 

OXXS,CXUjCX V»PS ,CXSU,CkSV pOR ACTUAL SECOND MOMENT MATRIX 
PARTITIONS GP,&CXXS,GCXJpGCX\/,GCXW,GPSiGCXSUp GCXSVp 
GCXSW FROM THE TIME OF THE LAST MEASUREMENT OR EVENT TO 
TME PRESENT TIME AND TO UPDATE THESE MATRIX PARTITIONS 
IF A MEASUREMENT IS 10 OF PROCESSED 


calling sequence* call 


ARGUMENTS* NR I 

IFLAGl I 

ICODE I 

UO I 

VO I 

i 

GCXW I 

GCXSW I 

' P I 

CXXS I 

CXU I 

, cxv 1 

. PS I 

cxsu I 


GNAVMINR, IFLAGl p I CODE ,U 0, VO p GCXW , GCXSW, 
P,CXXS,CXU,CXV,PS,CXSU,CXSV,Q,R) 

NUMBER OF ROWS IN THE OBSERVATION MATRIX 

-1 FOR assumed COVARIANCE PROCESSING 
=2 FOR ACTUAL SECOND MOMENT PROCESSING 

=0 FOR UPDATE 
= 1 FOR PROPAGATION 


ACTUAL OR ASSUMED DYNAMIC CONSIDER 
PARAMETER 2ND MOMENT MATRIX 


ACTUAL OR ASSUMED MEASUREMENT CONSIDER 
PARAMETER 2ND MOMENT MATRIX 


ACTUAL POSITICN-VELOCIT Y STATE /, IGNORE 
PARAMETER 2ND MOMENT MATRIX 


ACTUAL SOLVE-FOR PARAMETER / IGNORE 
PARAMETER 2ND MOMENT MATRIX 


ACTUAL OR ASSUMED POS I T I ON- V E LOC I TY 2ND 
MOMENT MATRIX 
2ND MOMENT MATRIX 

ASSUMED OR ACTUAL POS IT ION- VELOC I TY STATE 

/ solve-for parameter end moment matrix 

ASSUMED OR ACTUAL PQ S I T T 0 N- VF LO C I T Y STATE 
/ DYNAMIC D0N5I0ER PARAMETER 2NU MOMENT 
MATRIX 

ASSUMED OR actual POS IT ION- VE LO C I TY STATE 
/ measurement CONSIDER PARAMETER 2ND 
MOMENT MATRIX 


assumed or actual solve-for PARAMtTER 

COVARIANCE OR 2ND MOMENT MATRIX 


assumed or actual 

/ DYNAMIC CONSIOFR 
MATRIX 


SOL VF-FOR 
PA RARE ir R 


PAR 

'N 


A M 1. T t R 
) NT 
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CXSV I ASSUMED OR ACTUAL SOLVE-FOR PARAMETER 

f MEASUREMENT CONSIDER PARAMETER 2ND 
MOMENT MATRIX 

Q I ASSUMED OR ACTUAL DYNAMIC NOISE 2 NO MOMENT 

MATRIX 

R I ASSUMED OR ACTUAL MEASUREMENT NOISE 

2ND MOMENT MATRIX 


SUBROUTINES SUPPORTED* 

ERRANN 

SETEVN 

GUIDM 

PREO GENGID PROBE 

SUBROUTINES REQUIREO* 

GAINI 

GAIN2 




LOCAL SYMBOLS! AKW 

INTERMEOIAFE 

ARRAY 



OS 

INTERMEDIATE 

ARRAY 



ES 

INTERMEDIATE 

ARRAY 



FS 

INTERMEDIATE 

ARRAY 



I END 

NRrl 





NOIM4S 

NDIM4 

VALUE STORAGE 



N1 

NOIMl 

-1 




SUM 

INTEKMEOI ATE 

VARIABLE 



SW 

INTERMEDIATE 

ARRAY 



COMMON COMPUTEO/USEDI 

AK 

AL 

AM 

AN 

CXSUP 


CXSVP 

CXUP 

CXVP 

CXXSP 

G 


CCUV 

GCUR 

GCVM 

GCXSHP 

GCXMP 


HPHR 

IGAIN 

GH 

H 

HALF 


JPR 

ZERO 

PP 

PSP 

S 

COMMON USED! 

NDIMl 

N0IM2 

NDIM3 

NOIM4 

PHI 


TXU 

TXW 

TXXS 
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GNAVM Analysis 

Subroutine GNAVM propagates and updates (at a Q>eeBurement) both 
assumed (or filter) covariance matrix partitions and actual 2nd 
moment matrix partitions. The equations programmed in GNAVM are 
independent of the filter algorithm employed to generate gain 
matrices. 


The covariance and 2nd moment matrix partitions manipulated by 
Q1AVH are defined as follows: 


P =■ E(3t S ] 


C ° E 




P - E 
8 


r* *'i 

L S 8 


C “ E(x Q 1 

XU 8 

s 


C “ E[x u ] 

XU 


C - E[X v’^] 

XV ■’ 


C = E[x w^] 
xw ^ ^ 


C - E[x V ] 

XV ‘ 8 

8 


C - E[x * 1 . 

X w ‘■8 ■* 

3 


( 1 ) 


The following matrix partitions are tised in GNAVM, but are not 
changed in GNAVM: " 


C - E[u v^] 
uv 

C = E|u 
uw ■* 

- E[v 9 ^] (2) 

U = E[u u^] 

V « E[v 
W “ Efft . 
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GNAVM-2 


In these definitions S, u, and 0 representc respectively » 

the estimation errors in position/velocity state, eolve-for param- 
eters, dynamic consider parameters, measurement consider parameters, 
and Ignore parameters. Ignore parameteza, of course, are not de- 
fined when assumed (or filter) covariance matrix partitlona are 
being propagated or updated. Furthermore, the assumed C has been 
set to zero. 


The equations used to propagate covariances or 2nd moment 
from time to are summarlzedt 


P. , - I 4 >P^ + e C'^’^ +0 c'*'^ + e c'^ ^ 

k+1 \ k XX XX XU XU. xw xw. 

8 8« IC tC I 

k 




. T - T “ T 

+ c 8 ^ + c +c , 

k+l 


- + + +T +T 

C » 4>C + e P + 6 C + 0 C 

XX XX XX S, XU X U xw X V 

\+l “k “ “ k ‘ k 


c" ' - +0 c"*" +e u+e 

XUj^^l XUj^ XX^ X^U XU o xw uw^ 


c“ = + e c'*’ + 0 C + 6 

XV, , , XV, XX X V XU uv xw w 

k+l k s a , o o 

\ k 


c “ +0 + 6 C +0W 

xw, . , xw, XX X w XU uw xw o 

k+l k as, o 

k 


®k+l ®k 


C “ C 

x_u x„u 

k+l k 


matrices 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 
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C - C 

XV XV 

® k+1 ® k 


( 10 ) 


c ■ c 

XV XV 

® k+1 ® k 


( 11 ) 


In these equations ( ) indicates immediately prior to processing 

a measurement; ( )^» Immediately after. The state transition ma- 
trices over the interval [t^, are indicated by 4>, 0^ , 0^^, 

s 

and dynamic noise covariance or 2nd moment matrix is de- 

noted by 


Before covariance (or 2nd moment) matrix partitions can be updated 
at a measurement, the measurement residual covariance (or 2nd moment) 
matrix, defined by 


“^k+l " ^[^k+1 ^+l] 

must be computed. The required equations are summarized 


( 12 ) 


“^k+l ■’ “®k+l ^\+l ^^k+1 **\+l ®k+l 


(13) 


Vi - "k+1 « c 


XX 


T . - 
M + C 


k+1 


‘\+l 


T 

G + C 


XV, 


+ C" (14) 


k+1 


XV 


k+1 


B « P” + C“^ + C“ + C” + C” 

k+1 8, XX XU XV XV 

Vi “ k+1 “ k+j * k+1 


II ” + C"^ + U 0^ + c" + C" (16) 

^+1 ° '"^O 


k+1 


E. , , = C” + C M + C N + V L + C G 

k+1 XV, XV vw o uv 

k+1 ® k+1 ° 


(17) 
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W HP « T T' T 

^ = W N + C H + C M + C L + C 

K+1 O XW, , , X W VW UW 

k+1 8 1 , , o o 

k+1 


In these equations H, M, G, L, and N represent observation matrix 
partitions, and represents the measurement noise covariance 

(or 2nd moment) matrix. 


Gain matrices and are also required before covariance 

(or 2nd moment) matrix partitions can be updated. These are not 
computed in GNAVM but are obtained by calling either subroutine 
GAINI or GAIN2, depending on which recursive estimation algorithm 
is desired. 


With ^k+1 

used in the updating process: 


the following equations are 


“ P" 

k+1 k+1 


\+l ^ “ ^+1 '*■ *St+l “^k+l \+l 


XX 


* C 


XX 


Vi® - ^k+1 ^ \+i Vi ®k+i 


k+1 


k+1 


(19) 

( 20 ) 




XV, , , XV, , , 

k+1 k+1 


’^k+l ™k+l 


P ■= F 
®k+l ®k+l 


k+1 


x^u 


X u 


- \+l 

(21) 

/ 

- 

(22) 

T 

- \+l " 

(23) 

. ®k+l \+l ®k+l 

(24) 

T 

- ®k+l " 

(25) 


k+1 


k+1 


(18) 
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+ - T 

C - C - , B 

XV XV k+1 

® k+1 ® k+1 


+ - T 

C = C - 

X w X w k+1 

® k+1 ® k+1 


(26) 

(27.) 


It should be noted that propagation equations (3) through (11) are 
also used to propagate both assumed control covariance and actual 
2nd moment matrix partitions over the time Interval separating two 
successive guidance events. The update equations, of course, are 
not used in this situation. 
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ENTER 


Are assumed or actual 
covariances to be propagated? 


Actual 


Assumed 


Set ND4 to zero 


Set ND4 = NDIM4 


Propagate all covariance matrix 
partitions forward from t|^ to 

tj^+;^ using equations (3) through 

( 11) . Summarize Pj^+p 


Is a measurement time? 


Evaluate equations (13) 
through (18) to compute 
measurement residual co- 
variance matrix 



Equate covariance matrix par- 
titions at t|^^j to the cor- 
responding partitions at 


RETURN 


6NAVM Flow Chart 
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GPRINT-a 


S 

P 

c 

A 


UBROUT 

INE 

GPO TNT 


URPOSE 

8 TO 

PRINT 

ACTUAL 

ESTIMATION E 

alling 

SEQUENCE * 

GALL 

GPRINT (IFLAG* 

RGUMF.NTS 8 

I FLAG 

I 

=3 PRINT ACT 
GUIDANCE 
=10 PRINT act 
STATISTIC 
=2 PRINT ACTU 
ST ATISTICS 



T IMM 

I 

TIME TO BE PR 


RRQt< STATISfICS 
TIMM) 

UAL STATISIICS AT A 
EVENT 

UAL ESTIMATION ERROR! 

S 

al estimation error 
AT A prediction EVENT 

INTEG 


SUBROUTINES SUPPORTEOI PREO ERRANN SETEVN 
subroutines REQUIRED* MOMENT 


LOCAL SYMBOLS* A 

0 

OUM 


HOLLERITH WORD -AFTER- 
HOLLERITH WORD -BEFORE 
INTEKHEOIATE VECTOR 


EXSTSV TEMPORARY 
EXTSV TEMPORARY 


STORAGE for EXST 
STORAGE FOR EXT 


ROW INTERMEDIATE VECTOR 

ZZ INTERMEDIATE VARIABLE 


COMMON USED* CXSUP 

EMRES 

EXSI 

GCXSV 

GCXH 

GU 

N0IH2 

PSP 

XSL 


CXSVP 

CXUP 

CXVP 

CXXSP 

EU, 

EV 

EH 

EXI 

EXST 

EXSTP 

EX TP 

GCXSU 

GCXSH 

GCXSHP 

GCXU 

GCXV 

GCXWP 

GCXXS 

GP 

GPS 

GV 

GW 

JPR 

NDIMl 

NOIM3 

NDIM4 

NR 

PP 

RPR 

TRTM2 

XIG 

XL A8 

XU 

XV 




I 
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SUnHOUTINE GQOOMP 

PURPOSEI TO COMPUTE ACTUAL EXECUTION ERROR STATISTICS 
CALLING S£(aUENCE8 CALL GQCOMP « V V , EE e EEE » F V , QJ 


ARCUHFNTS8 

FE 


I ACTUAL COMMANOEO yfLOCITT CORRECTION 
I MEANS OF ACTUAL EXECUTION ERRORS 
EEF I 2ND moments OF ACTUAL EXECUTION ERRORS 
EV 0 EXPECTED VALUE OF ACTUAL EXECUTION ERROR 
Q 0 ACTUAL EXECUTION ERROR 2ND MOMENT MATRIX 


ORTEDI 

GENGIO 


F ACTR 

INTERMEDIATE 

VARIABLE 

9H0P 

MAGNITUDE OF 

VV VECTOR 

RHOP2 

RHOP»*2 


VI 

V(l)'»^2 


M2 

M12>**2 


V3 

V<3) 


M>: 

V1*V2'^V3 


XI 

INTERMEDIATE 

VARIABLE 

XMUP 

INTERMEDIATE 

VARIABLE 

ZETA 

INTERMEDIATE 

VARIABLE 
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GQC0MP Analysis 

Subroutine GQC0MP computes the actual execution error mean anti 
2nd moment matrix for use in the generalized covariance analysis 
of a guidance event. The actual execution error 6&V' Is assumed 
to have the form 


AV’ 


6AV’ = k’ AV' + s' 

j ' i 


AV’I 

j' 


( 1 ) 


where k' denotes the actual proportionality error; s', the actual 
resolution error; actual pointing error; and 

AVj, the actual commanded velocity correction. 

The means of the three ecliptic components of 6AVj are given as: 

( 2 ) 


E[6AV] = 

(it' 

+ 


AV' + 

X 



p / 

X 

E[6AV] = 

(k' 


K) 

AV’ + 

y 



p' / 

y 

E[6AV] = 
z 

(k' 

+ 


AV' - 

z 


o'AV' 6a' + AV' AV' 6B' 
y X z 


AV' AV* 66' - p'AV’ 6a' 
y z X 


(3) 

( 4 ) 


where p' * UV’i, u' = (AV'^ + AV'^]*^, and 6a' and 66' are the 

X ' y ^ 

actual pointing angle errors, and both E( ) and (~) indicate mean 
values, 

The actual execution error 2nd moment matrix Is defined by 


Q’ " E 


6AVj’^ 


(5) 


the elements of matrix are given as: 

q" « AV'2 + 5^’ + aV'2 AV'2 66' 66' + 

11 X \ y X z ' 


\ 2AV' / 

2p' AV AV' AV' 3Tv'j c' fPAV' AV' AV' 

xyz /u xz 


Ci) 
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Q’ “ AV'^ + -tV 66' 66' + p'2 AV'2 6 u' 6a' - 

22 y M ^ \ y z X 


2p' AV’ At^' AV* 6a' 66') + — C' ( 
X y 2 / g \ 


AV AV 

y 2 


: 6B’ - p' AV' 6a’V 

z X / 


( 6 ) 


OU ■“ AV’^ + g’2 66' 66’ - 2A^’ g' c' 66' 


‘33 


^l2 


(7) 


Q' - C A^' AV' + V 
21 X y u 


2AV' AV’ AV' 66' - p'/aV’2 - a^’^'Iao' 
X y z \ X y / 


1 

rrz-' 


- p'7 6V' AV' 6a' 6a' + p' AV' (aV' 2 - AV'^] 6a' 
X y . 2 \ y X / 


AV' AV' AV’^ 6’6’ 66' 
X y 2 


66 ' + 

( 8 ) 


Q[3 - = C' AV; AV; 4. 


AV' / _ _ X __ 

— T P'AV’ 6a' + AV' AV' 66’)- g' A^' 66’ 

g’ Vy xz J ^ - X 


- p'AV' 6a' 66' - AV’ AV' 66' 66' 
y X z 


(9) 


Q'3 - Q^2 “ ^ 


AV' 

z 


(aV* AV' 66' - p'AV' 6a'|- g’ AV' 66' 

\ y 2 / y 


+ p'AV 6a' 66' - AV' aV 66' 66' 
X y 2 


( 10 ) 


where 


and 


c = k' k' + -I- k' s' + 

p P 


( 11 ) 




( 12 ) 
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SUBROUTIME guio 

PURPOSE; COMPUTE THE GUIDANCE MATRIXi THE VARIATION MATRIX, AND THE 

TARGET condition (JUST BEFORE) COVARIANCE MATRIX AT A 
MIOCOURSe GUIDANCE EVENT 

calling SE(3UENCEJ call GUI0(RF»IGP«TEXN,GA,A0A) 

arguments; rf state at time of event 

, IGP guidance policy (=1 for FTA* =2 FOR VTA) 
tevn time of event 
GA GUIDANCE matrix 
ADA variation MATRIX 

SUBROUTINE SUPPORTED: GUIDM 

subroutines REQUIRED: CSTART EIGHY MATIN NTM RSIM SHIFT 

STMPR zermat 

^LOCAL SYMBOLS: ALF ARRAY OF HOLLERITH CONSTANTS FOR PRINT 

OUHl INTERMEDIATE 2X2 MATRIX 

EGVL VECTOR OF EIGENVALUES 

lERR ERROR FLAG RETURNED BY FILE READER 

PHI3 INTERMEDIATE 3X3 MATRIX 

ROT ROTATION MATRIX FOR VTA COMPUTATIONS 

ROW TEMPORARY STORAGE VECTOR 

SQP TEMPORARY STORAGE VECTOR 

YPM INTERMEDIATE VALUE 

ZPM INTERMEDIATE VALUE 

ZTM INTERMEDIATE VALUE 

COMMON COMPUTEO/USED: DELTM "PHI TRTMl 

COMMON USED: FISAVE FISAVE 
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GUID Analysis 


Subroutine GUID is called at a midcourse guidance event at t . in the 

J 

error analysis mode to compute three primary quantities for the selected 
midcourse guidance policy. These three quantities are the variation 
matrix 7]^ the target condition covariance matrix prior to the velocity 

correction , and the guidance matric Two midcourse guidance 

policies are available: fixed- time- of- arrival (FTA), and variable- time- 
of-arrlval (VTA). Both are linear impulsive guidance policies having 
form 

Av. = r 8 X. 

J j J 


where AVj is the commanded velocity correction, and SX^ is the estimate 


of the spacecraft position/velocity deviation from the targeted nominal. 
The relevant equations for each guidance policy will be summarized below. 


The variation matrix 77 . for FTA guidance relates deviations in space- 
craft state at t. to position deviations at the final time t„, and 
j F 

is given by 


where 


K i ^7] 


t.). 


1 i ”21 

.J 


is the upper half of the state transition matrix 


The guidance matrix for FTA guidance is given by 

-1 


r. = 

J 


0 I -1 

I 


The variation matrix for VTA guidance relates deviations in state at 
tj to deviations from the nominal normal to the impulsive insertion vel- 
ocity vector. Consequently, the state deviations at the nominal target 

time, t j.are rotated to a coordinate system whose z-axis is along the 
£ 

input delta-V, REXV, and only the upper 2x6 partition of the rotated 
tji(tj,, tj) matrix becomes 7] 

"via ' ; “] 

where 

A = (r 0 1 B = R0nl 

^ J 2x3 L J2x3 
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Next the guidance matrix is obtained, by 


^ * 
VIA 


[ T T-1 * T T-1 

) A I -B (BB ) B 


= -B^(BBT)"^ 

' VTA 


Whichever guidance policy is used to obtain Tf, and F., the target 
condition covariance matrix is confuted as 


W" = P " 

J J c , j 

J : 

Where P is the control covariance matrix immediately prior to the 
to the guidance event. 

Finally, saved values are restored and the trajectory file is re- 
initialized for use by the file reader when called from other sub- 
routines. 



GUID-3 
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110 



Compute and print Wj , 
its std . dev . 's and 
correlations 



Partition 0^ 


± 

q - [- 0/1 01-1] 


GUID-4 




I 
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SUBROUTIME GUIOM 

PURPOSE CONTROL EXECUTION OF A GUIDANCE EVENT IN THE ERROR 
ANALYSIS PROGRAM 

CALLING SEQUENCE! CALL GUIOM 

SUBROUTINES SUPPORTED! ERRANN 

SUBROUTINES RFOUIRED: COPREL OTnO GUIl) EIGHY JACOBI GNAVM 

I'VSTaT NTm PSIM STMPR CSTART SAVMAT 


LOCAL SYM0OLSS 


ADA VARIATION MATRIX 

ftMAK INTERMEDIATE VARIABLE USED TO FIND MAXIMUM 
eigenvalue of VELOCITY CORRECTION 
COVARIANCE MATRIX IS MATRIX) 

CXSUl STORAGE FOR CXSU KNOWLEDGE COVARIANCE 

CXSVl STORAGE FOR CXSV KNOWLEDGE COVARIANCE 

CXUl STORAGE FOR CXU KNOHLEOGE COVARIANCE 

CXVl STORAGE FOR CXV KNOWLEDGE COVARIANCE 

CXXSl STORAGE FOR CXXS KNOWLEDGE COVARIANCE 

DUMl INTERMEDIATE VARIABLE 

DUM VECTOR SUM OF UPDATE AND STATISTICAL 
VELOCITY CORRECTIONS 

EGM maximum EIGENVALUE OF S MATRIX 

EGVCT ARRAY OF EIGENVECTORS 

EGVL ARRAY OF EIGENVALUES 

FXEC f.KECUTIOM ERROR COVARIANCE MATRIX 

txv EXPECTED VALUE OF VELOCITY CORRECTION 

GA GUIDANCE MATRIX 

GAP intermediate ARRAY EQUAL TO GA TIMES P 

ICODE INTERNAL CONTROL FLAG 

rcooE? internal control flag 

'^GP ’/'^T.DCQURSE GUIDANCE POLICY CODE 
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SQP 

TRS 


ISPHC TEMPORftRV STORAGE FOR ISPH 

nkP OF IKIAXZHUH EIGEf^VALUE OF S 

@yT SPACECRAFT VELOCITY RELATIVE TO TARGET 

PLANET IN PLANETO'CENTRIC EQUATORIAL 
COORDINATES 

P2 S70R4GE FOR P CONTROL COVARIANCE 

Pi STORAGE FOR P KNOWLEDGE COVARIANCE 

RF NOMINAL TRAJECTORY STATE AT GUIDANCE EVENT 

RHO MAGNITUDE OF STATISTICAL OELTA^V 

ROW INTERMEDIATE VECTOR 

INTERMEDIATE VECTOR 
TRACE OF S MATRIX 

VEIG MATRIX TO BE DIAGONALIZED 

Z INTERMEDIATE ARRAY 


COMMON 

COHPUT ED/USED S 

CXSU6 

CXU 

ISPH 

P 

COMMON 

COMPUTEOS 

OELTM 

COMMON 

USED: FOP 

FOV 


Q 

SKALE 


uo 

VO 


cxsu 

CXSVG 

CXSV 

CXUG 

CXVG 

CXV 

CXXSG 

CXXS 

NGE , 

PG 

PSC 

PS 

TG 

XG 




TRTMi 

XI 




ICDT3 

NDIMI 

NDIM2 

NDIM3 

ONE 

SIGALP 

XF 

SIGBET 

ZERO 

SIGPRO 

SIGRES 

TWO 
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GUIDM Analysis 


Subroutine GUIDM is the executive guidance subroutine in the error analysis 
program. In addition to controlling the computational flow for all types of 
guidance events, GUIDM also performs many of the required guidance computations 
itself. 

Before considering each type of guidance event, the treatment of a general 
guidance event will be discussed. Let t. be the time at which the guidance 
event occurs. _Before any guidance evenf^can be executed, the targeted 
nominal state X., knowledge covariance , and control covariance P must all 

^ s- j j 

be available, where ( ) indicates values immediately before the event. The 
first two quantities are available prior to entering GUIDM. However, GUIDM 
controls the propagation of the control covariance over 'the interval [^tj.j^, tj], 
where tj_^ denotes the time of the previous guidance event. 

The next step in the treatment of a general guidance event is concerned with 
the computation of the effective velocity correction and the execution error 
covariance. In the error analysis program, only a statistical velocity 
correction can be computed. The effective velocity correction AVj is then 
used to compute the execution error covariance matrix ^4. A summary of the 
execution error model and the equations used to compute Qj can be found in 
the subroutine GQC0MP analysis section. 

The last step is concerned with the updating of required quantities prior to 
returning to the basic cycle. An assumption underlying the modeled guidance 
process is that the targeted nominal remains unchanged at a guidance event. 

The knowledge covariance is updated using the equation 

0 9 0 

. -I- - 

0 B TJ. 

if an impulsive thrust model is assumed. 

In the impulsive case, the control covariance is updated simply by setting 



This equation is a direct consequence of the assumption that the targeted 
nominal state is always updated at a guidance event. 

Each specific type of guidance event involves the computation of other 
quantities not discussed above. These will be covered in the following 
discussion of specific guidance events. 
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1. Midcourse Guidance 

Linear midcourse guidance policies have form 


Av. 


N. 


= r, 


SX. 


where the -subscr ipt N indicates that this is the velocity correction required 
to null out deviations from the nominal target state. This notation is used 
to differentiate between this type of velocity correction and velocity 
corrections required to achieve final insertion. Linear midcourse guidance 
policies are discussed in more detail in the subroutine GUID analysis section. 


Subroutine GUIDM calls GUID to compute the guidance matrix, p-, aud the 
target condition covariance immediately prior to the guidance ^vent, Wj , and 
then uses f- to compute the velocity correction covariance Sj, which is 
defined as: 

s. . E , 

and is given by the equation 

S-. = r. (p ~ - s p ~ ) rl 


where s is a (real) scalar, input by the analyst; generally, 0. $ s 1 . 

This equation assumes that an optimal estimation algorithm is employed in the 
navigation process, since the derivation of this equation requires the 
orthogonality of the estimate and the estimation error. 


In the error analysis program AVj^ _ is never available since no estimates & Xj 

are ever generated. Only the ensemble statistics of are available which 

means only a statistical or effective velocity correction "E [AVn J can be 

computed. In the STEAP error analysis program, this effective velocity 
correction is assumed to have form: 


The 

the 


magnitude is given by the 

analysis of subroutine DUSTAT. 


gj, 

ail 

Lee -Boa in 


ana lytic 


solut ion as 


descr ibed 


in 
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The direction of the effective velocity correction is assumed to coincide with 
the eigenvector corresponding to the maximum eigenvalue of S . . This eigenvector 
is denoted by a . • ^ 


After the updated control covariance has been computed, the 

+ , ^ 

condition covariance matrix W. following the guidance correction 
using the equation: ^ 


target 
is computed 


w"!" = »?. 

J J •= 


T 

j j 


where variation matrix ’Ij has been previously computed in subroutine GUID. 


2. Final Insertion 

A final insertion event describes an insertion (into halo orbit) which may be 
accomplished by an impuslive or by a finite burn. If the burn is impulsive, 
the expected correction has been input as REXV. This vector is used to ^ 
compute the executive error matrix, Qp, just as AVjj . is used to compute Qj . 

If the burn is finite, GUIDM calls NTM and PSIM to compute the state to state 
transition matrix, <£(tp, tg), and the control to final state transition 
matrix 0 (tp, tg) , where the control parameters are the pointing angles a and/3 
and the thrust magnitude T. Both the control and the knowledge covariances 


are propagated by 

(tp, 

tg) 


II 

'J’ (tp , 


d)(tp, tg)^ 

II 

4>(tp, 


cj)(t t )^ 

tg; 


The execution error matrix for the finite burn is the 6x6 matrix 

Qp “ 6 (tp, tg) Up 0 (tp, tg) 

„ n 2 

where is the diagonal matrix whose diagonal elements are ® , snd 

a T J Uq is input to the program. In either, the impulsive or the finite case, 
both the knowledge and the control covariances are updated by adding the 
execution errors. These results are printed, and the end of the final 
insertion event is the end of the ERRAN run. 
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GUIDM FLOWCHART 




■3 7 0 /. 


TYPE OF^ 
GUIDANCE EVENT 


/. Ann 





GUIDM-5 
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HGIONS-A 


SUBROUTINE HGIDNS 

PURPOSE: TO COMPUTE THE CHANGE REQUIRED TO THE CONTROL VARIABLES 

FOR TARGETING 

CALLING sequence: CALL HGIDNS 

arguments: 

NONE 


LOCAL SYMBOLS: 
TFP 
El 
E2 
E4 
DF 

AETA 

BETA 

ETaI 

STmB 

STmd 

temxi 

PXU 

TXU 

atar 

aer 

DElTAV 

xfeq 

tchang 


TIME FROM PERIAPSIS (SECS) 

SEMI-MAJOR AXIS (KM) 

eccentricity 

argument OF PERIAPSIS (DEG) 
computed injection time (DAYS) 

SENSITIVITY PARTIALS OF TARGETS l»RT POSITION CHANGES 
SENSITIVITY PARTIALS OF TARGETS WRT VELOCITY CHANGES 
TARGETING MATRIX GAMMA 

PARTION OF STATE TRANSITION MATRIX OF POSITION 
CHANGES WRT VELOCITY CHANGES 

PARTION OF STATE TRANSITION MATRIX OF VELOCITY 
CHANGES WRT VELOCITY CHANGES 
SENSITIVITY matrix 

PARTIALS OF STATE WRT FINITE BURN CONTROLS 
AT THE END OF THE BURN 

PARTIALS OF STATE WRT FINITE BURN CONTROLS 

at target time 

ACTUAL TARGET VECTOR ON CURRENT NOMINAL TRAJECTORY 
target error VECTOR (DESIRED-ACTUAL) 

CONTROL UPDATE VECTOR 

INJECTION state VECTOR IN EARTH EQUATORIAL- 
GEOCENTRIC COORDINATES 

ACCUMULATED TOTAL CHANGE IN CONTROL VECTOR 


COMMON used: 




DL 

STM 

DTOL 

SPD 

XB 

ACCTH 

PERT 

SMU(4) 

TB 

IBURN 

LIBR 

H 

XF 

UMAX 

PI 

ECEQ 

ASTM 

OTAR 

RPO 


COMMON computed: 




XL 

ALPHAO 



ITOL 

betao 



iter 

TBURN 



KWIT 





« 
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HGIDNS FLOW CHART 


^ START ^ 


Zero out accumulated change 
first time through 


Compute actual target from 
nominal trajectory 


Compute target error 




AER = 
DTAR - 
A TAR 



Error greater than tolerance 







HGIDNS-2 


Terminate case if maximum 
iterations are exceeded 


Is this a finite burn case 


Generate update control variables 




RETURN 




HLAUCH-A 


subroutine hlauch 


PURPOSE; TO COMPUTE THE INJECTION TIME 

f 

CALLING sequence: CALL HLAUnCH 1 X « D J ) 

arguments; 

X I INJECTION STATE VECTOR IN ECLIPTIC-GEOCENTRIC 

COORDINATES (X» Y» Z (KM) AND XDOT , YDOT* ZOOT 
(KM/SEC) ) 

DJ I/O INPUT AS DESIRED INJECTION JULIAN DATE* OUTPUT 
AS actual injection JULIAN DATE 

LOCAL SYMBOLS; 

REEJD JULIAN date JAN 1* 1950.0 

OLA EARTH EQUATORIAL DECLINATION OF ANGULAR MOMENTUM 

VECTOR (RAD) 

SOLA SIN OF *DLA» 

RhilSR LATITUDE OF LAUNCH SITE (RAD) 

SIGR LAUNCH AZIMUTH RADIANS 

SSIG SINE OF *SI6R* 

SPhI sine OF *PHILSR« 

CPHI COSINE OF *PHILSR» 

CSiG COSINE 0 *SI6R‘ 

TD DAYS FROM REFJD TO INJECTION DAY 

SHa GREENWICH HOUR ANGLE (DEG) 

CThE cosine OF »ThE» 

SThE sine of *THE» 

the RIGHT ascension AT LAUNCH (RAD) 

TL launch TIME ON DAY OF LAUNCH (DAYS) 

FL true anomaly of LAUNCH SITE (RAD) 

PSIB angle BETWEEN LAUNCH AND INJECTION (RAD) 

TC COAST TIME (SEC) 

T0 TIME BETWEEN LAUNCH AND INJECTION (DAYS) 

SIG launch AZIMUTH (DEG) 

TI injection TIME OF DAY (DAYS) 

X£0 INJECTION STATE VECTOR IN EQUATORIAL-GEOCENTRIC 

coordinates (X. Y* Z (KM) AND XDOT, YDOT, ZDOT 
(KM/SEO) 

ELEMS 6 vector OF ORBITAL ELEMENTS (SEMI-MAJOR AXIS 

(KM), eccentricity, INCL INAT ION ( DEG) * LONGITUDE 
OF ASCENDING NOOE(OEG), ARGUMENT OF PERIAPSIS 
(DEG), TRUE ANOMALY (DEG) > 

k(HAT 3 element ANGULAR MOMENTUM VECTOR IN EQUATORIAL- 

GEOCENTRIC COORDINATES 
ZAXIS PSEUDO POLE VECTOR (0., 0., 1.) 


SUBROUTINES REQUIRED; 

OM^TPY CALJUL 

OUXV DUNIT 

OVFCRO DANGV2 

CAREL 
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subroutine HlAUCH (CONTINUED) 
COMMON USED: 


Fi 

THELS 

ECEQ 

PS 11 

PHILS 

PI 

PSI2 

THEDOT 

TWOPI 

TIMI 

RPRAT 

RPD 

Tim? 

SIGMAL 

SPD 
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HLAUCH Analysis 


HLAUNCH computes the injection time from the injection state and the 
launch profile parameters input by the user. 


The injection state is first rotated from the ecliptic coordinate system, 
which was input, to the earth equatorial system. 


RcA = <J>Eceq Rec 
VcA = cfJEceq Vec 


( 1 ) 


The unit normal to the launch/orbit plane is then calculated in earth 
equatorial coordinates as 



X V , 
CA CA 


X V 
CA 


ca| 


( 2 ) 


The inclination of the orbit plane i (= arc cos Wz) should equal the 
desired input value. The orbit plane inclination must equal or exceed 
the latitude of the launch site to permit a coplanar parking orbit 

and transfer orbit as Indicated in Figure 19. In the case that 
sin i sin ^ the launch azimuth is defined by 


sinlj^ 


cos i 
cos 


(3) 


and the solution with 0^ 2^ ^90 degress is selected. In this case 

the parking orbit nominal is identical to that of the transfer plane 
given by (2). 




Pa rking 
Orb it 

T ransfer 
Plane 


Figure 1 Transfer Plane/Parking Orbit Geometry 


4-99 



HLAUCH-2 


If Isin ll < the parking orbit and the transfer orbit 

cannot be coplanar (Figure 16). In this case the parking orbit is de- 
fined to be in the plane have a launch a'zimuth of L = 90 deg. , 

jj 

containing the closest approach radius vector Rp., and nearest the 

Lin 

transfer plane, (Note the alternate parking orbit plane in Figure lb 
which also satisfies the first two of these requirements). The unit 
normal to the parking orbit plane is given by 


W = 
P 


CA 


I^CA^ 


X V 
V 


(A) 


where is the velocity vector at the injection point in the parking 

orbit, V is given by 
P 


V 

P 




-cos 0 
P 

sin 5 

P 

cos £ 

p 

-sin 0 

P 

sin E 

-Sin 

sin 6 

P 

cos 2 ; 

P 

+COS Q 

p 

sin E 


cos 5 

cosl 




p 


where (0,5) are the equatorial right ascension and declination of the 
P P 

periapsis position R^. For the specific parking orbit plane having 
= 90 deg, including R„. , and nearest the transfer plane must 

Ll La P 

satisfy 

cos 0 

sin = (6) 


cos 


sgn (cos X ) = sgn 
P 


V 


CA 


V (0) 
P 


where 0<X 180 deg and where the equation (5) is used, 

P 

Thus the unit normal to the parking orbit plane may be computed by 
either (2) or (4) and the launch azimuth is either given by (2.47) or 

= 90 deg. In either case the remaining calculations proceed as follows. 

the right ascession at launch is defined by. 
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cos ^ 


W sinO- sin£ 

X i-i 1 j 



+ W cos X 

1 


L 


sin 


W 


sin X 



L 

1 


W cos X T 
X L 


(7) 


The launch date input by the user is recalculated as the integer day 

(0^ ut) closest^to the inital date input by the user. The Greenwich 
hour angle at 0 ut of the launch date is then 


GHA = 100°. 07554260 + 0°, 9856473460 T^ (8) 

0 

2°. 9015 X lO'^^ 

d 

The launch time on the day of launch is 

, <61 ,gt . 2. 

w 


where w is the rotation rate of the launch planet and is the 
longitude of the launch site, both being read in as input. 


The unit vector toward the launch position is the 

= (coscj5^ cos cos sin(^j^|, sin't'^^j) (10) 

The true anomaly of the launch site f.^ calculated as: 

L) 


^L~ \ ' ^GA 
sin = \ ■ 


( 11 ) 


The angle between launch and injection is 




( 12 ) 


The coast time t^ may now be computed 



(13) 


where ^ and ^ s.re the angle of ° the first and second burns and 
ktj) is the inverse parking orbit coast rate, all of which are input. 
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The time between launch and injection is therefore 


S - h '^2 'c , 

where tj^ and t^ are the input time durations of the first ' 
and second burn 

The injection time is then 

tj = tj^ + tg (15). 


I 
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SUBROUTINE HPRELN 

PURPOSE: TO INITIALIZE CONSTANT SAND DEFAULT VALUES. READ INPUT 

DATA. AND CALCULATE THE ZERO ITERATE GUESS= 

CALLING sequence: CALL HPRELM 


arguments: 

NONE 

LOCAL symbols: 

BTIME DURAION of FINITE BURN (DAYS) 

I0IAS INPUT FLAG INDICATING THAT A BIAS VECTOR 

IS TO BE ADDED TO THE LIBRATION STATE PRIOR 
TO TARGETING 

IZERO INPUT FLAG INDICATION THE SOURCE OF THE 

INITIAL CONDITIONS 

XInT injection STATE VECTOR TO BE INTEGRATE IF 

ITMAX=0 

DTaR array of target VALUES 

ZOAT INITIAL GUESS OF SPACECRAFT VELOCITY AT LIBRATION 

POINT. IF IZER0=5 

RE RADIUS VECTOR FROM THE SUN TO THE EARTH 

SUBROUTINES REQUIRED: 

DZERO TRNSPS 

CALJUL DSHIFT 

PECEQ EPHGT 

DVSMLT DAVECT 

DVECRD hzerit 

common computed: 


ALPHA 

FI 

itmax 

NB 

PERT 

atry 

IBTYPE 

, LAUNCH 

NBOO 

PHILS 

beta 

lOISK 

LIBR 

NFR 

PSIl 

otol 

ipre 

msglvl 

NPOINT 

PSI2 

RLIBR 

TOUR 

TIMl 

XISP 

ITOL 

rprat 

THEDOT 

TIM2 

Z8IAS 

KWIT 

scmass 

thels 

tmpr 

IDON 

RETRO 

SIGMAL 

THRMAG 

TP 

ITER 

LNON 

IBURN 

ECEO 




tburn 

XF 




TB 

volp 




dl 

XL 
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HPRELM Flowchart 




HPRELM-2 


! 


Initial velocity of S/C 'at libration point 
comes from input 


Initial velocity comes from tables or 
Lamberts solution 
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htrjty-a 


SUBROUTINE HTRJTY 

PURPOSE; TO CONTROL THE TRAJECTORY GENERATION PHASE 


CALLING sequence: CALL HTRJTY 

arguments; 

none 


LOCAL SYMBOLS; 

OJ 

TBIAS 

TCOAST 

SH 

TTo 

TACT 

ttopre 

ST 

TTP 

DJp 

RMaG 

vmag 

RA 

DEC 

delv 

KSWB 

IRONG 

KPOINT 

ISTOP 

kshp 

X 


JULAIN date of INJECTION 

BIAS date added TO PRINT OUT TIME 

COAST time 

SIGN OF INTEGRATION STEP 
PRINT time (DAYS) 

print time of next special PRINT POINT (DAYS) 
PRINT time of previous PRINT POINT (DAYS) 

SIGN OF PRINT POINT 
DAYS FROM INJECTION 
JULIAN DATE OF PRINT POINT 
MAGNITUDE OF RADIUS VECTOR (KM) 

MAGNITUDE OF VELOCITY VECTOR (KM/SEC) 

RIGHT ASCENSION (DEG) 
declination (DEG) 

MAGNITUDE OF IMPULSIVE BURN DELTA V VECTOR 
FLAG TO INDICATE INITIATION OF FINITE BURN 
NUMBER OF ERRORS ENCOUNTERED DURING TRAJECTOR 
PRINT phase 

CURRENT SPECIAL PRINT POINT NUMBER 
FLAG TO INDICATE ARRIVAL OF STOPPING CONDITIONS 
FLAG TO INDICATE THAT CURRENT PRINT POINT IS 
A SPECIAL PRINT POINT 

state vector 





subroutines REOUIREDs 


HLAUCH 

EPHGT 

DVSMLT 

davect 

seti 

PSTART 

MSTART 

OR0INT 

COrtELL 

□shift 

C'alJUL 

eval 

dsvect 

DABSV 

DRD 

DMa’TPY 

ORBEND 

BURN 



COMMON USED; 





NBOD 

TMPR 

thrmag 



NB 

PI 

XISP 



iarray 

RPD 

scmass 



itol 

SPO 




iter 

H 




COTAR 

LAUNCH 




lOlSK 

ECEQ 




NPOINT 

MS6LVL 




TP 

planet 





« 
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SUBROUTINE HTRJTY (CONTINUED) 
COMMON COMPUTED/USED: 


XL 

VOLP 

OL 

IBURN 

XB 

lOON 

TB 

LNON 

XF 


TOUR 



COMMON computed: 
BSTM 
ASTm 

stm 

accth 



KWIT 

ALPHA 

beta 

tburn 
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HTRJTRY FLOWCHART 


Check if targeting has converged 


Check if finite burn is being targeted from 
impulsive burn 


Should run target based on actual launch 
profile 


Check if actual launch profile has been 
generated 
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Generate actual Launch profile 


Reset convergence flag &; iteration counter 


Generate projected launch profile 


Initialize integrator for COAST phase of 
trajectory 



0 


i 

r 

CALL HLAUCH & 

INITIALIZE 

FLAG 


r 

CALL SI 
CALL P£ 
CALL Mf 

;ti 

START 

.TART 
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Check if orbit file is being 
written 


Skip special print logic if 
first point 
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Check if there are special print points 


Has last special print been used 


Is special print beyond nominal print point 


Is special point beyond last point 






HTRJTY-5 


Set special print flag, set current print 
point to special print point jdecrement 
index for NOMNAL point 


Set prev ious ■ pr int point to current print 
point 


Check to see if current print point is 
within section 


Reset print point to terminal value & set 
termination flag 


Check if current print point is outside 
integration range 
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HTRJTY-6 


Are we in the burn phase 


Is this an impulsive burn 
case 


Is current print point 
within coast phase 


Integrate to burn 
initiation 


Reset flags 



Integrate to 
print point 
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HTRJTY-8 
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Have 10 bad print point 
accumulated 


Set burn flag 


Determine finite burn 
parameters from impulsive- 
burn data 








HTRJTY-10 



Is this an impulsive burn 


Initialize integrator 


Integrate burn phase of trajectory 


Turn burn flag off. & restart integrator 


Integrate coast phase of trajectory 
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V 


HTRJTY-11 
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SUBROUTINE HZERIT 

PURPOSE: TO COMPUTE THE INITIAL FOR TARGETING WHEN IZERO = 6 OR 7 

CALLING sequences CALL HZER I T ( I ZERO , RE ♦ 01 e RL » TUUR e VL ) 


arguments: 

IZERO 

RE 

01 

RL 

TOUR 

VL 


I INITIAL CONDITION FLAG 

1 injection RADIUS VECTOR (KM) 

2 JULIAN DATE OF INJECTION 

I RADIUS VECTOR TO LIBRATION POINT (KM) 
I FLIGHT TIME (DAYS) 

0 velocity VECTOR AT LIBRATION POINT 


LOCAL SYMBOLS: 

alp ANGLE BETWEEN RL AND VL (RAO) 

VLM MAGNITUDE OF VL 

Th transfer angle (RAD) 

Sma SEMI-MAJOR AXIS OF TRANSFER ORBIT 

E ECCENTRICITY OF TRANSFER ORBIT 

P SEMI-LATUD RECTUM OF TRANSFER ORBIT 

NTaB NUMBER OF LOOKUP TABLES 

KTaB VALID LOOKUP TABLE 

NP NUMBER OF POINTS IN LOOKUP TABLE MINUS ONE 

ZAXIS PSEUDO-POLE VECTOR / 

FTaB LOOKUP tables 

' FLIM LOOKUP table LIMITS 

NLIM NUMBER OF POINTS IN LOOKUP TABLES 


SUBROUTINES REQUIRED: 

OUXV LAMBRT 

OANGMO LOOP 

OVCOMB DANGV2 

COMMON USED: 

PI retro 

HPO ATHY 

SPo NFH 


SMU 




\ 





HZERIT-1 


HZERIT FLOWCHART 


Are initial conditions generated from 
Lamberts solution 


Find correct table for this flight, time 


Does flight time lie within any table 


Interpolate table for velocity 
magnitude and rotation angle 


Compute initial guess of libration point 
velocity 
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0 


HZERIT-2 


Compute transfer angle 


Check for transfer greater 
than one full revolution 


Compute velocity at libration 
point from Lamberts solution 





■ 

THETA - 
f (RETR0, 

Rl, R2> 

1 

r 


c 



RETURN 


) 


YES 
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SUBROUTINE InITLC 

PURPOSE: TO initialize CONSTANTS 

ARGUMENTS: 

NONE 


LOCAL SYMBOLS; 

rln character string of, planet Names 

subroutines reouired 

OShIFT 

COMMON computed: 

pi smu 

TWOPI RSOI 

RPO RP 

SPn MONTH 

AKMPAU planet 
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PUR^OSSS OF a REat SVMMETRIC 70 OiaCOMfiS. 

F0R51 BY a SUCCESSION OF PtaNE ROT&TIONS TO ANNIHILATE 
THE ®FF°DSA60NAL ELEMENTS AND SUBSEQUENT COMPUTATION 
OF THE EIGENVALUES ANO EIGENVECTORS Of THAT MATRIX 


OALLING SEaUEMCEi 
ARGUMENTS k 


SUBROUTINES SUPPORTED I EIGHT 

GUIDM 


CALL JAC0e!(A9U2tV»N«>F00S 

I MATRIX TO BE OlAGONALl ZED (HILL BE 
DESTROYEOJ 

0 VECTOR OF EIGENVALUES (LENGTH N$ 

0 MATRIX OF EIGENVECTORS (N BY N OIHENSXON) 

1 DIMENSION OF SQUARE MATRIX A 

I FINAL OFF-OIAGONAL ANNIHILATION VALUE 

PRESIM SETEVN 


GUISin GUISS 
GUIO PRED 


LOCAL SYMBOLS S 


DEL 


INTERMEDIATE VARIABLE 
INTERMEDIATE VARIABLE-A UPSP) 
INTERMEDIATE VARIABLE**AUPJP» 
INTERMEDIATE VARI ASLE-A ( JPJP) 
INTERMEOIATE VARIABLE 
DIFFERENCE IN ELEMENTS OF A 


KR 


DIMENSION OF A 
KR -8- 1 


IMTEKMEOIATE 

INTERMEDIATE 


?a 

VHP 


LARGEST OFF-DIAGONAL ELEMENT 

INTERMEDIATE 
ONE TWO 
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ZERO 



JACOBI- i 


JACOBI Analysle 


The Jacohi ineCliod tubjectK a real, symmetric matrix A to a sequence ol 
tran; forma tion s l>aeecl rjii a rotuLioii matrix: 


cus (b 
K 

sin <b 

K 


- 8 In ^ 

K 

cos (f> 

K J 


where all other elements of the rotation matrix are Identical with the 
unit iTiatrlx. After n multiplications, A Is transformed Into: 


A' 



o;\a 0^,... 0^ 


If is chosen at each step to make a pair of off-diagonal elements 

xero, then A' will approach diagonal form with the eigenvalues on the 
diagonal. The columns of 0^^ ^2‘‘’^N to the eigenvectors of A. 


The angle of rotation 0 is chosen in the following way. If the 
four entries of 0- are in (l,i), (j,l) and (j,j) then the corres- 

K . 1 \ 

ponding elements of Oj^ A Oj^ are 


If 


(3 


2 2 
b , = a cos + 2a sin 0 cos p + a sin {3 
li i-i Ij JJ 

^ij “ ■ ®il^ ^ ^ aj^j(cos2(J - sin20) 

b . = a . . sln^(3 - 2a . sin ^ cos (J + a . . cos^lJ 

JJ 11 1| JJ 


is chosen so that tan 2(3 “ 2a, /(a., - a .) then 

ij 11 JJ 


ij ji 


Each multiplication creates a new pair of zeros but will introduce a non- 
zero contribution to positions zeroed out on previous steps. However, 
.successive matrices of the form ^2*^ *^l”^ ^ 1 ^2 approach the 

required diagonal form. 

Reference: Scheld, Frances; Theory and Problems of Numerical 

Analysis, McGraw-Hill Book Company, Inc., New York, 
1968. 
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Set Initial V iiiatriK tc> imlLy. 
Set W2(l) - A(l) 


Set Tl - AKS(A(2)). 

Scati upper 

off- diagonal elements 

of 

matrix A 

hy rows to find ^*reateet 

element In 

a1)uoLute value. Set 
tills e Lement . 

Tl 

equal to 


Set IREDO ■ U. Scan upper off- 
dlagonal eleraenti! of matrix A by rows 
until n value greater than Tl lb 
found. Pivot on this element. 


Compute rotation angle . 
Set IREDO - 1 


Compute eigenvectors and 
diagonalize matrix A, 


Tl .LE. F0D7 


Place diagonal from 
A into W2 


Tl - Tl ★ O.OUl 


RETURN 













LAMBRT-A 


SUBROUTINE LamBRT 

PURPOSE: TO SOLVE LAMBERTS PROBLEM FOR TRANSFERS LESS THAN TWO PI 

CALLING sequence; CALL LA0RT (R1 »R2f TS»THETA»XMU» VI » » V2» AsE»P) 
ARGUMENTS; 

R1 I POSITION VECTOR AT DEPARTURE POINT 

R2 I POSITION VECTOR AT ARRIVAL POINT 

TS I TRANSFER TIME 

theta I TRANSFER ANGLE (RAO) 

XMU I GRAVITATIONAL CONSTANT OF CENTRAL BODY 

VI 0 VELOCITY VECTOR AT DEPARTURE POINT 

V2 0 VELOCITY VECTOR AT ARRIVAL POINT 

A 0 SEMI-MAJOR AXIS OF TRANSFER ORBIT 

E 0 ECCENTRICITY OF TRANSFER ORBIT 

P 0 SEMI-LATUS RECTUM OF TRANSFER ORBIT 

LOCAL SYMBOLS; 

RIM 
R2M 
URl 
UR2 
C 

UC 
CM 

tpara 

SUBROUTINES REQUIRED; 

dabsv dangmd 

ACOSH DVCOMB 

dsvect 

OVSDIV 

COMMON USED; 

PI 

TWOPI 
RPD 

COMMON computed: 

KUIT 


MAGNITUDE OF R1 
MAGNITUDE OF R2 

UNIT vector IN DIRECTION OF R1 
UNIT VECTOR IN DIRECTION OF H2 
R2 - R1 

UNIT VECTOR IN DIRECTION OF C 

magnitude of C 

PARABOLIC TIME OF FLIGHT 
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LAMBRT Analysis 


LAMBRT solves Lambert's problem for transfer angles less than 360°. 

Given: 1) Initial and final position vectors 

(R^ and Rj) 

2) The time of flight (t£) 

3) The transfer angles (9) 

The initial and final velocity vectors (V^ and V 2 ), for a two-body conic 
trajectory connecting these points are found. 

The problem is solved using the following steps: 


S = B>2 - 


( R. 


i- R^ + C )/2.0 


= iL (s 3/2 - g (s - C )3/2) 
3 


2 2 

where g = sign ( - 9 ) 


and 


0<9 < 2-n- 


For tj > t , the conic is an ellipse and the following transcendental 
equation mBst be solved for A 


t£ = ^ 


, ( A- sin A ) - g (B - sin B) 

) 1 


1-cos A ^ 1 

where : 

s ( 1 - cos B) = (s - c ) (1 - cos A ) 

and 


0 4 ^ 


2 7T 


0 A ^ 

0 4 B 4 A 

M = gravitational constant of central body 

The semi major axis of the ellipse may now be calculated as: 
a = s / (1 - cos A ) 


Also calculate: 


; g 2 = sign ( 7T 


a2) 


For t£ < t^, the conic is an hyperbola, and r must be solved for in the 
transcendental equation: 
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LAMBRT-2 


tf =/ - s ^ (siah r-r) - g (sinh S - 5 ) 

'' cosh r-1 - 

where : 

s (cosh S - 1) = (s - C ) (cosh ■ r - 1) 

and 

0 "v S $■ r ^ ijo 

The semi major axis of the transfer hyperbola is 
a = s / (1 - cosh r ) 


and set = +1 


For both elliptical and hyperbolic transfer orbits the initial and final 
velocities are now calculated from the following equations: 




= yTZ: (A + B) 

^ 2 



(A - B) 


V, 

= V 


+ V 


-1 

c 

= V 


p 

- V 

1 

-2 

c 


p 

2 


For the derivation of these equations, the reader is referred to Astronautlcal 
Guidance by R. H. Batten. 
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SUBROUTINE LOOP 

PURPOSE: TO SOLVE LAMBERTS PROBLEM FOR TRANSFERS GREATER THAN TWO PI 

CALLING sequence: CALL LOOP ( R 1 « R2 . TS » TmET A ♦ XMU f ATRY t NER i V I t V2 » A » E » P ) 

ARGUMENTS: 

Rl I POSITION VECTOR AT DEPARTURE POINT 


ner 

I NUMBER OF FULL REVOLUTIONS BEFORE ENCOUNTER 

nfr 

I NUMBER OF FULL REVOLUTIONS BEFORE ENCOUNTER 

ATpy 

I INITIAL GUESS 

LOCAL SYMBOLS: 


RMl 

MAGNITUDE OF Rl 

WM2 

MAGNITUDE OF R2 

CM 

MAGNITUDE OF C 

C 

R2 - Rl 

cthi 

COSINE OF TRUE ANOMALY AT POSITION 1 

CTh2 

COSINE OF TRUE ANOMALY AT POSITION 2 

sthi 

SINE OF TRUE ANOMALY AT POSITION 1 

STh2 

SINE OF TRUE ANOMALY AT POSITION 2 

THSCI 

TRUE anomaly AT POSITION 1 

THSC2 

TRUE anomaly AT POSITION 2 

VM2 

MAGNITUDE OF VI 

VMl 

MAGNITUDE OF V2 


subroutines REQUIRED: 

DSVECT OABSV SL0020 

OUNIT DUXV DVSMLT 

DVCOMB DANDMO 

COMMON USED: 

PI 

HPO 


& 


'4-129 



LOOP-1 


LOOP Analysis 


LOOP solves Lambert's problem for transfer angles greater than 360°. 

Given: 1) Initial and final position vectors 

and R 2 

2) The time of flight (t^) 

3) The transfer angle (9^^) 

The initial and final velocity vectors (V^ and V 2 ) for a two-body elliptical 
conic connecting these points are found. 


The problem is solved using the following step. 

C = ^ - 

s = ( R^ + R 2 + c' )/2.0 

Q = ~ cos (0/2) 

s 

where 9 = (8jj) nod 2 •r- 

The following sets of equations must be solved for X by interating 



E = - 1 

y = (-E)^ 

K = 

z = (1 + KE)^ 
f = y (z-Qx) 
g = xz-QE 
X = tan'^ (f/g) 
d = m + X 

The orbital elements are then 
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LOOP -2 



2 

e 


(2^3)^^ Qz (s- Ri ) - X (s - R, ) 


( 1 - 


' ■ ■ ' } + ■ 


For the derivation of the previous equations, the reader is refered to NASA 
Technical Note D -5368 (A Unified Form of Lambert's Theorem, by E. R. Lancaster 
and R. C, Blanchard), 

The initial and final velocities are now calculated using the method; 


= a 


(1 - e2) 


cos 


cos 


= (P - R^ )/(e ) 

N2 = (P - R2 )/(e R2 ) 


sin = (cos N^ cos 9 - cos, N2)/sin 0 


N2 = (Nj^ + 9 ) mod 2 tt 


cos ^ 


aP 


( 2 a - R. ) 


2 sh 


sin Fj^ = e * sin / (1 + 2e cos N^ + e ) 


VI 


P(- 


cos = ( “^ 


aP 


\ 


(2a - R2) 


sin = e + sin N2/(l + 2e cos N2 + e^)^ 


2 2 
V2" = ( — 


- ^ 


w = 

§1 ( 

^1 

X R 2 

where 

®1 

= 

+ 1.0 

and 

®1 

= 

- 1.0 

A 



/S 


w 

X 


Y T 

0 




0 ^ 0 


7T e <.2 r 
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= VI (sin + cos F^ ) 

= V2 (sin F2 ^^2 ■'■ ^2 "^2 ^ 
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ROUTINE MAIN 

PURPOSES ENTRY POINT TO PROGRAM nOMNAL 
LOCAL symbols 5 

OPTION CHARACTER STRING READ FROM INPUT 

halo CHARACTER STRING 'HALO* 

SUBROUTINES REQUIRED: 

INITLC 

HPRELM 

HTRJTY 

PRELIM 

HGIONS 

trjtry 

GIOANS 

COMMON USED: 

KWIT 

COMMON COMPUTED: 

IPRE 
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MAIN-1 


MAIN FLOWCHART 


Check if halo orbit mission 





MAIN-2 



4 - 13 ^ 


I 






MAIN-3 


Initialize data for halo orbit mission 


Case termination? 


Generate halo orbit trajectory 


Case termination? 


Target halo orbit trajectory 


Case termination? 
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MATIN-A 


SUTROUTINt MATIN 

PURPOSES TO COMPUTE TME INVERSE OF A MATRIX. 
calling StQUENCEl CALL MATIN(A,R»N> 

ARGUMENTS A(N,N) I MATRIX TO 0E INVERTED 

R(N,N> 0 RESULTANT INVERSE OF MATRIX A 
N I DIMENSION OF A AND R 


SUHROUTINES REQUlREQl 
LOCAL SYMBOLS* AL 

ALDAR 

B 

DETR 

G 

IX 
KR 
MIXI 
HIXJ 
MIXL 
S 

X 

XOFF 


NONE 

ACLL) ♦ S tINTERMEDIA T£ 

intermediate variable 
intekmeoiate vector 
intermediate vector 
intermediate vector 
intermediate vector 
dimension of a 
intermediate variable 
intermediate variable 
intermediate variable 

INTERMEDtATE VARIABLE 

intermediate variable 

INTERMEDIATE VARIABLE 


VARIABL E) 
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MEAN -A 


SU0ROUTIN£ HEftN 

PURPOSE* TO PROPft&ATE ANO UPOATE MEANS OF ACTUAL STATE ANO 

PARAMETER OEtflATIONS ANO ACTUAL STATE ANO PARAMETER 
estimation errors 

calling sequence* CALL ME A N ( EX T P » EXS TP , I FL AG ,I FL AG 1 , NR ) 

ARGUMENTS* EXTP I STATE OEVTATIONS OR ESTIMATION ERRORS 

EXSTP I 30LVE-FOR PARAMETER DEVIATIONS OR 

ESTIMATION ERRORS 

t 

IFLAG I =1 FOR UPDATE 

=2 FOR PROPAGATION 

I 

f 

IFLAGl I =1 FOR DEVIATION MEANS 

=2 FOR ESTIMATION ERROR MEANS 



NR I NUMBER OF ROWS IN THE OBSERVATION MATRIX 


SUBROUTINES SUPPORTED* 

ERRANN 

SETEVN 

PROBE 

GENGIO 

PR ED 

LOCAL SYMBOLS* IGO 

INTERNALLY SET FLAG 



SUM 

intermediate 

STORAGE 



ZERO 

value 

0.0 




COMMON COMPUTED/USED* 

DUME 

EU 

EV 

EH 

EXIF 


EXSIP 





COMMON USED* 

AK 

AL 

AM 

AN 

G 

■ 

H 

NDIMi 

N0IM2 

N0IM3 

NDIM4 


PHI 

S 

TXU 

TXW 

TXXS 
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MEAN Analysis 

Subroutine MEAN propagates and updates actual estimation error 
means over the time interval [t. . t. separating twc aucceaaive 

measurements or events. The equation^ programmed In MEAN are In- 
dependent of the filter algorithm employed to generate gain ma- 
trices. Gain matrices are assumed to have been computed during a 
prior call to subroutine GNAVM. The propagation equations pro- 
grammed in MEAN are also used to propagate actual deviation means 
over the time interval separating two successive guidance events. 

The update equations, of course, are not used in this situation. 

The actual estimation errors for position/velocity state, solve-for 
parameters, dynamic consider parameters, measurement consider param- 
eters, and ignore parameters are defined, respectively, by the 
following: 


\+l ” \+l ■ \-t-l 

it =• it - X 

®k+l ®k+l \+l 

^+1 ^+1 

\+l " ^k+1 ■ ^'k+l “ 

O ^ tJ a 

k+l k+1 k+1 



( 1 ) 

( 2 ) 

( 3 ) 

(A) 

( 5 ) 


where (*) indicates estimated values, and x, x , u , v , and w 

SCO o 

are the actual deviations from nominal. 


Only the means of ic and x are paopagated and updated since the 

s 

means of u, v, and w are constant, The propagation equations are 
aummkrized 



where $,0 ,6 , 

XX^ XU 


and 0 are state transition maCrlcea over 
xw 
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MEAN-2 


Before the means of x and x can be updated at a meaaureoient , the 

B 

mean of the measurement residual must first be computed using 

X 1 + Gu + Lv + Nw (8) 

. ®k+lj o o o 
where H, M, G, L, and H are observation matrix partitions. 

The update equations are summarized as: 

^[*k+i] ■ ®[\+i] * \+i ■ '[“^k+i] 

"[vi] " ^[\^] * \*x • 

where K, and S, are the filter gain matrices. 
k+1 k+1 

To propagate actual deviation means requires that x and x be 

s 

replaced by ft and x^, respectively, in equations (6) and (7), and 
that the minus signs in equations (6) be replaced with plus signs. 


(9) 

( 10 ) 
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MEAN Plow Chart 


ENTER 
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MEN0-A 


SUBROUTINE MEMO 

PURPOSE I COMPUTE ASSUMED AND ACTUAL MEASUREMENT NOISE COVARIANCE 
MATRICES IN THE ERROR ANALYSIS PROGRAM 

CALLING SEUUENCEI CALL MENO ( MMCODE » I CODE ) 

ARGUMENT I ICOOE I INTERNAL COOE USED TO DISTINGUISH BETWEEN 

THE TWO alternatives LISTED ABOVE 

MMCODE I MEASUREMENT MODEL CODE 

SUBROUTINES SUPPORTEOI ERRANN 

COMMON C0MPUTE08 R RPR 

COMMON USED* IMNF MNCN IGMNF GMNCN 
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MEN0 Analysis 

The linearized observation equation employed by the navigation 
process la given by 


where 6Y, la the measurement deviation from the nominal meaaure- 

A 

mento ia the augmented observation matrix, is the augmented 

state deviation from the nominal augmented state, and is the 
assumed measurement noise. 

The function of subroutine M£N0 is to compute the assumed measure- 
ment noise covariance matrix 




„ T 
^ \ \ 


If IC0DE ° 0. The constant measurement noise variances associated 
with all available meaaurement types are stored in the vector MMCN. 
Subroutine MEN0 selects the appropriate element from this vector 
to construct 


If IC0DE 0 the actual measurement noise covariance matrix 

\ 




where is the actual measurement noise, is computed instead. In 

this case subroutine MEN0 selects the appropriate actual measure- 
ment noise variances from the vector GMNCN to construct R^. 

The accompanying flow chart Indicates the computational flow for 
computing An identical procedure is used to compute R^. 
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MENO Flow Chart 


8ETUKN 


) 


BETUKH ^ 



ENTER 


) 


Only constant noise 

NO 

y 

1 available 

at present. 



L 

‘ 



k 

Zero out matrix. 

\ 

! 


Is measurement 
noise constant? 


YES 


Compute the 

R matrix for 

3 8tAr-plAtiet angle meas- 

urements . 





MMCPDE=7 


NO 


10 


,11,12,13 




1,2, 3, 4,5, 6. 7, 8 
YES 


MMC0DE even? 




Compjte the matrix 

for a range- rate meas- 
urement from the 
appropriate station. 


Compute the matrix 

for a range and range- 
rate measurement from 
the appropriate station. 


> i *C ' 


10 


^ RETURN ^ 


Compute the matrix 

for the appropriate 
star-plane angle meas- 
urement . 


Compute the matrix 

for an apparent planet 
diameter measurement. 
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MOMENT-A 


SUHROUTINr MOMENT 

PUPPOSEI TO CONVERT AN ARBITRARY NON-SQUARE 2ND MOMENT MATRIX TO 
THE ASSOCIATED CORRELATION MATRIX PARTITION AND PRINT IT 
ALSO COMPUTE AND PRINT EIGENVALUES, EIGENVECTORS , AND 
HYPERELLIPSOIOS 

calling sequence* call MOMENT(NI,N2,EXYT,EX,EY,CORW,CORM1,ABL,I1, 

I2,IFLAG,IF2) 


ARGUMENTS* 

N1 

I NUMBER OF ROWS IN 2ND MOMENT MATRIX 


N2 

I number OF COLS IN 2ND MOMENT MATRIX 


EXYT 

I N1 BY N2 2ND MOMENT MATRIX OF X AND Y 


EX 

I N1 VECTOR MEAN OF X 


EY 

I N2 VECTOR MEAN OF Y 


CORN 

I 2ND MOMENT MATRIX CORRESPONDING TC VECTOR 
X OF DIMENSION N1 


CORWl 

I 2ND MOMENT MATRIX CORRESPONDING TO VECTOR 
Y OF DIMENSION N2 


ABL 

L VECTOR OF ROW LABELS CORRESPONDING TO 
CORWl 


I 1 

I ROW INDEX MAXIMUM 


12 

I COL INDEX MAXIMUM 

■■■ .,!■ 

IFLA6 

I =0 DO NOT COMPUTE EIGENVECTORS, ETC. 


IF2 

I =0 DO NOT COMPUTE STD. DEV. 

subroutines 

SUPPORTED* 

CPRINT GENCIO 

SUBROUTINES 

REQUIRED* 

EIGHY 

LOCAL SYMBOLS* OUT 

INTERMEDIATE ARRAY 


PEIG 

INTERMEDIATE VECTOR 


ROW 

INTERMEDIATE VECTOR 


SQP 

INTERMEDIATE VECTOR 


SQPl 

INTERMEDIATE VECTOR 


VEIG 

INTERMEDIATE VECTOR 

COMMON USED* 

FOP FOV 
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/ 


M0MENT Analysis 

T 

Subroutine M0MHNT transforms an arbitrary 2nd moment matrix E[xy ] 
Into a correlation matrix and. If x ■■ y. Into a vector of standard 
duvlattons. The transformation consists of two steps; 

T 

1) Transform E[xy ] into the covariance matrix 


cov (x,y) - E[xy*^] - E[x] ♦ E(y^]; 


2) Transform cov (x,y) into the correlation matrix having 
correlation coefficients 


where 




iJ Oj 


1 i i 


'iJ 


E[x^ yj] 


E{x2] 


2ih 


2 1^ 


E[y^l 


Subroutine M0MENT writes out the correlation matrix and. If they 
exist, the standard deviations. Subroutine M0MENT can also compute 
and write out the eigenvalues, eigenvectors, end hyperellipsold 
of cov (x,y) if X ■ y. 
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SUBROUTINE NTM 

PURPOSE: TO SHIFT THE LAST STATE TRANSITION MATRIX OBTAINED FROM THE 

FILE INTO PHIOLDf TO CALL THE FILE READER ANO OBTAIN 
. THE NEW state VECTOR ANO STATE TRANSITION MATRIX* PHiNEW* 
(I.E., TO SET UP COMMON BLOCK PHISAV FOR USE BY SUB- 
ROUTINE PSIM). OR TO FLAG THE ERROR ENCOUNTERED WHILE 
TRYING TO READ THE FILE. 

CALLING SEOUenCE: CALL NTM ( RF ♦ PH I » K SECT ) 

ARGUMENTS: RF STATE VECTOR OBTAINED FROM FILE 

PHI state transition matrix 

KSFCT INDEX OF SECTION OF FILE TO BE READ 
=1 FOR COAST SECTION 
=2 FOR FINITE BURN SECTION 

LOCAL SYMBOLS: lERR ERROR FLAG RETURNED BY FILE READER 

TSEC TIME IN SECONDS PAST THE INITIAL TIME ON THE FILE 

subroutines REQUIRED: GETCOW 

COMMON computed: PHINEW TOLD 

COMMON USED: TRTMI TRTMB DELTM PhIOLO TM 
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subroutine ORBENO (ORBINT entry POINT) 

PURPOSE: TO WRITE A «FINAL* RECORD TO THE SEQUENTIAL ORBIT ERE 


CALLING sequence: CALL ORBENO 

ARGUMENTS: 

none 


local symbols: 

ifrn 

LOGICAL 

FILE NUMBER 

common USED; 

T 

SKI 

NSECTN 

H 

XVDO 


XDD 

SVl 


SKI 

SV2 

s 


COMMON computed: 

ielevn 



0RBINT«A 


subroutine ORRINT 

PURPOSES TO initialize THE SEQUENTIAL ORBIT FILE WITH PARTIALS 
CALLING sequence; CALL ORBINT 

arguments; 


NONE 

LOCAL SYMBOLS; 

IFRN 

COMMON USED; 


LOGICAL FILE NUMBER 


VMDIC 

nstate 

H 

HMSIC 

kstate 

XDO 

AEINT 

IPART 

sxi 

SPINT 

GM 

SX2 ' 

pvint 

nsectn 

XVOD 

ObLINT 

T 

svi 


SV2 
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ORBINT Analysis 


Subroutine ORBINT contains five entry points; ORBINT, ORBWRT, ORBEND, 
CSTART and GETC0W. The purpose of these entry points is to trans- 
mit and retrieve information on a sequential orbit that contains partials. 
Entry point orbint writes a 1048 byte logical header record followed 

by a 1056 ^q byte second logical header record. This routine is called 

only once at the start of trajectory generation. After the header re- 
cords have been written the ORBWRT logic block is executed. This block 
writes a data record containing the latest 11 acceleration 

vectors, first and second cowell runs, current integration timcj step size 
and section number. This logic block is repeatedly used during trajectory 
integration by a call through entry point ORBWRT. After the^ trajectory 
has been completely integrated a call to entry point ORBEND will write 
a data record of any remaining acceleration vectors and set an end-of- 
file mark. 

Entry points '■CSTART and GETCOW are used to retrieve information from the 
sequential trajectory data file written by ORBINT/ORBWRT/ORBEND, An 
initial call to entry point CSTART prior to trajectory generation will 
read the two headon records and the first data record. Each successive 
call to entry point GETCOW will read successive data records (if required) 
and call subroutine INTP to interpolate the trajectory data for the 
requested time period. 
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SUBROUTINE ORBWRT (OHBINT ENTRY POINT) 

PURPOSE; TO WRITE RECORDS TO THE SEQUENTIAL ORBIT FILE 
CALLING sequence: CALL ORBWRT 

arguments: 

NONE 

LOCAL SYMBOLS: 

IFRN LOGICAL FILE NUMBER 

COMMON USED; , 

T SX2 NSECTN 

H XVDD 

XDD SVl 

SXl SV2 

COMMON computed; 

IELEVN 
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SUHROUTINE PtCEQ 

PURPOSE « TO COMPUTE THE MftTRIX DEFINING THE TF? ANSF ORMft T I OM FROM 
PLANET CENTERED ECLIPTIC COORDINATES TO PLANET CENTERED 
EQUATORIAL COORDINATES AS A FUNCTION OF THE PARTICULAR 
PLANET AND TIME. 

calling sequences call PECEQ (NP,D,ECEQ> 


ARGUMENT NP I CODE OF PLANET 

0 I JULIAN OATEs EPOCH 1900, OF REFERENCE TIME 

ECEQ<3,3> 0 COORDINATE TRANSFORMATION MATRIX FROM 

PLANETOCENTRIC ECLIPTIC TO PL ANETOCENTRIC 
EQUATORIAL COCRDINATES 


SUBROUTINES REQUIRED? EULMX ORB 
LOCAL SYMBOLS* 


AHCGC COORDINATE TRANSFORMATION MATRIX FROM 
GEOCENTRIC ECLIPTIC TO GEOCENTRIC 
EQUATORIAL COORDINATES FOR EARTH - FROM 



ECLIPTIC 

MOON 

TO orbital plane COORDINATES FOR 

CSDECL 

COSINE OF 

DECL 

CSE09L 

COSINE OF 

EO0L 

CSINM 

COSINE OF 

INM 

CSNOM 

COSINE OF 

NOOEM 

CSRASC 

COSINE OF 

RASC 

DECL 

DECLINATI 

ON OF TARGET PLANET POLE 

DGTR 

CONVERSION FACTOR FROM DEGREES TO RADIANS 

ED 

Julian date , epoch 4713 b.c. 

E OBL 

OBLIQUITY 

OF ECLIPTIC 

I t4H 

INDEX 
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NODCM 

INDEX 

^8 0H'M 

-yNIT VECTOR normal TO TARGET PLANET 

ORBITAL PLANE 

\ 

P0AK 

CROSS PRODUCT OF POLE AND NORN 

PMflG 

MAGNI TUDE OF P6AR 

com 

UNIT VECTOR aligned WITH TARGET PLANET 
POLAR AXIS 

PCLM&G 

magnitude OF POLE 

Q0ARP 

CROSS PRODUCT OF POLE AND PBAR 

OMAT, 

MAGNITUDE OF QBARP 

RASC 

RIGHT ASCENSION OF TARGET PLANET POLE 

SNOECL 

SINE OF OeCL 

SNEOBL 

SINE OF EOBL 

SNINM 

SINE OF INCLINATION INM 

SNNOM 

SINE OF NODE NOM 

SNRASC 

SINE OF RASC 

TPRIM 

BESSELIAN DATE 

X I 

INTERMEDIATE VALUE 

X IQ 

INTERMEDIATE VALUE 

XL 

INTERMEDIATE VALUE 

XLQ 

INTERMEDIATE VALUE 
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PECEQ Analysis 

Subroutine PECEQ computee the coordinate transformation maxtrix 
A from planetocentrlc ecliptic to planetocentric equatorial co- 
ordinates for an arbitrary planet. 

The derivation of A for a planet other than the earth or moon 
will be summarized. Matrix A is defined by 

A = [x j Y 1 z]^ (1) ■ 

where X, Y, and Z are unit vectors aligned with the planetocentric 
equatorial coordinate axes and referenced to the planetocentric 
ecliptic coordinate system. Unit vector Z is aligned with the 
planet pole. Unit vector X lies along the Intersection of the 
of the planet equatorial and orbital planes and points at the 
planet vernal equinox. Unit vector Y completes the orthogonal 
triad and is given by 

Y = Z X X. (2) 

It remains to obtain expressions for X and Z. Let N denote the 
unit vector normal to the planet orbital plane, and let P denote 
the unit vector aligned with the planet pole Then 

Z = P (3) 


and 


P X 

Ip X n| 


( 4 ) 


The unit vector N, referred to the ecliptic coordinate system, is 
given by 


N 


sin , 1 sin Q 
-sin i cos Q 
cos i 


(5) 


where 1 and are the inclination and longitude of the ascending 
node, respectively, of the planet orbital plane. The unit vector 
P, referred to the ecliptic system is given by 
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I cos a cos 6 

P “ I cos L sin u cos 6 + sin t. sin 6 
I - sin c sin u CDS 6 + cos t sin 6 


( 6 ) 


where n and <‘> are Che righc ascension and declination, respec^ 
Cively, of the planet pole relative to the geocentric equatorial 
coordinate system, and c is the obliquity of the ecliptic. Ex- 
pressions for >1 and 6 for each planet were obtained from JPL TR 
32-1306, Cans tan and Related Information for Aatrodynamio Cal- 
culations t 1968, by Melbourne, et at. 

For the earth and the moon, the transformation matrix A is writ- 
ten as Che produce of two transformation matrices 


A = A 2 A^. 


(7) 


For the earth A^ is 


the identity matrix and A^ is given by 



0 

cos G 
sin e 


T 

0 

-sin c . 
cos e 


( 8 ) 


The following figure defines the transformations A^ and A„, 
using the definitions given. 


XYZ Ecliptic coordinate axes 

X Y Z Orbital plane coordinate axes 

000 

X^YqZq Moon's equatorial coordinate axes 

1 Inclination of moon's orbital plane Co ecliptic 

plane 

Right ascension of moon’s orbital plane to 
ecliptic plane 

6 Inclination of moon’s equatorial to orbital 

p lane 

6 Right ascension of moon's equatorial to orbital 

plane 


4-155 



PECEQ-3 


Z 



^ ^0 


The transformation from ecliptic to orbital plane coordinates 

is performed by rotating about the z-axis through an angle- Q and 
then about the resulting x-axis through an angle 1. Symbolically » 

= (Q about 3, 1 about 1). (9) 

The transformation A^ from orbital plane to equatorial coordinates 
can be written similarly as 


A^ « (6 about -3, 6 about -1) . 


( 10 ) 



PECEQ-4 


PECEQ Flow Chart 


ENTER ^ 




r: 


Call 0RB to compute the 
Orbital elements of the 
planet 


i 

’ 


Compute the obliquity of 
the ecliptic e 



: earth? ^ 



i No 


Compute sine 
inclination i 
of planet. 

and cosine of 
and node n 


ICompute the right ascension a 
land declination 6 of the 
iplahet pole. 



Compute the coordinate 
transformation matrix A. 


RETURN 3 


Compute transformation 
matrix A^. Set 

A2 ” I . 


L5 







PRED-A 


SUopnUTlMF PRE^ 


PURPOSE; CONTROL CALCULATIONS FOR A PREDICTION EVENT IN ER«AN 


CALLING sequence: CALL PPED 

subroutines REQUIRED: CORREL CSTART DYNO 

MEAN NTM °SIM 


EIGHY GNAVM GPRINT 

5AVMAT SHIFT STMPR 


LOCAL SYMBOLS: DUM2 

DUM3 


ARPAt OE eigenvectors 
ARRAY Of- EIGENVALUES 


DUM intermediate array 

EGVCT ARRAY OE EIGENVECTORS 

EGVL ARRAY OF EIGENVALUES- 

EXSTS TEMPORARY STORAGE FOR EXST 

tXST temporary storage FOP EXT 

ICOOE INTERNAL CONTROL FLAG 

IPR temporary storage FOR IPRINT 

peig tntermeoiate array 

PSAVE temporary storage FOR KNOWLEDGE COVARIANCES 



RE 

state 

VECTOR AT 

TPT 





PI 

state 

VECTOR AT 

TIME OF 

EVENT 




TPT 

time predicted 

TO 





VEIG 

N A T h 1 X 

Tu BE OT AGONAL I ZED 



COMMON 

computfo/used: 

CXSU 

cxsv 

cxu 

CxV 

cxxs 

OELTM 



GCXSR- 

GCXW 

IPrInT 

NPE 

P 

PS 



TRTmI 

XI 





COMMON 

USED: EH 

Exsr 

EXT 

FOP 

FOV 

GP 

ISTMC 


NO I M 1 

NO I 

N r.H M 3 

N T M C 

ONE 

Q 

TPT2 


' liO 

VO 

XF 






\ 
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PRED-1 


PRED Analysis 


Subroutine PRED executes a prediction event in the error analysis 
program ERRAN. At a prediction event, the knowledge covariance par- 
titions, and the estimated position/velocity deviations from the most 
recent nominal trajectory are propagated forward to t , the time to 

which the prediction is to be made. The knowledge covariance partitions 
are propagated using the prediction equations found in the GNAVM 
Analysis section. The estimate is propagated using the equation 


= «(t , t) <5x + e (t , t ) dx 

p p J XX p J s 


where ^ and 



are the state transition matrix partitions over the 


time interval 



The position and velocity partitions of the propagated knowledge co- 
variance matrix are diagonalized at time t^ and the eigenvalues, and 


eigenvectors are computed. 



PRED-2 


PRED FLOW CHART 


^ , ,E^ER ^ 


Increment prediction event counter and 
obtain time tp to which prediction is to 
be made. 


Save all knowledge covariance matrix 
partitions at t j . 


I 


Call NTM to compute the nominal 
trajectory at time tp. Call PSIM to 
compute the state transition matrix 
partitions over the time interval 


Si 


jcall DYN0 to compute the dynamic noise 
jcovariance matrix for the interval 

Iwrlte out the state transition matrix 
partitions and eigenvalues. 


I 


E li GNAVM to propagate knowledge 

variance partitions forward to time tp. 
Write out the knowledge correlation 
matrix partitions and standard deviations 
at time t„. 




PRED-3 


Restore knowledge covariance matrix 
par tic ions at t - . 


IGEN = L 


Save actual knowledge covariance matrix 
partitions at t . . 



RETURN 







PRELIM-A 


SUBROUTINE PRELIM 

PURPOSE; dummy link WITH NON HALO ORBIT OPTIONS 
CALLING sequence: CALL PRELIM 



PRINT3-A 


subroutine PR1NT3 

I 

PURPOSE: print PERTINENT INFORMATION AT SPECIFIED MEASUREMENTS 

CALLING sequence: CALL PR INT3%MMCODE » NR • 

arguments; MMCODE CODE FOR TYPE OF MEASUREMENT 

NR NUMBER OF ROWS IN OBSERVATION MATRIX 

subroutines required: correl STMPR STVCPR 

LOCAL SYMBOLS; 01 HOLLERITH CONSTANT 

02 HOLLERITH CONSTANT 

lA STATION number 

IB STAR NUMBER 

lONE =1 

ITHREE =3 

I TWO =2 

M INTERNAL MEASUREMENT CODE 

TRTM2 TIME OF MEASUREMENT 


COMMON USED; 

AK 

AL 

AM 

AN 

oeltm 

H 

mcntr 


Q 

R 

S 

TRTMl 

XIG 

XLAR 

XSL 


XU 

XV 
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subroutine PSIM 

PURPOSE; TO invert A 6X6 STATE TRANSITION MATRIX (USING ITS SYMPLECTIC 
CHARACTER) FROM T1 TO T2» AND MULTIPLY A STATE TRANSITION 
MATRIX FROM T1 TO T3 BY THAT INVERSE* AND THEREBY OBTAIN 
THE STATE TRANSITION MATRIX (STM) FROM T2 TO T3. 

CALLING sequence: CALL PS I M ( P3 1 . P2 1 . P32 ) 

ARGUMENTS: P31 STM FROM Tl TO T3 

P21 STM FROM Tl TO T2 

P3? STM FROM T2 TO T3 . 

subroutines REQUIRED: NONE 

t, 

COMMON USEO; none 
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PSTART-A 


SUBROUTINE PSTART 

PURPOSE; TO initialize THE STaTE PARTIAL HATRIX 
CALLING sequences CALL PSTART 

arguments: 

none 

LOCAL symbols: 

aB mean RADIUS OF CENTRAL BODY 

COMMON USED: 

INO(i) IBURN 

ICENT IND(16) 

X 

xo 

GM 


SUBROUTINES REQUIRED: 


ELEM 


POLAR 


PARTE 


COMMON computed: 


elemq 

NEQ 

orbel 

XV 

sphcoq 

XVD 

pvint 

PXPD 

AEINT 


spint 


OBLINT 
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savmat-a 


SUBROUTINE SAVMAT 

PURPOSE; TO STORE THE 034 VALUES OF ARRAY P IN ARRAY PI 

CALLING sequence; CALL SAVMAT (P«PU 

ARGUMENTS; P ARRAY TO BE SAVED 
PI STORAGE ARRAY 
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SCHED-A 


SUBROUTINE SCHEO 

PURPOSES TO DETERMINE RHAT TYPE' OF MEASUREMENT IS TO BE TAKEN 
NEXT AND AT NHAT TIME IT WILL OCCUR. 

CALLING SCOUENCEt CALL SCHEO ( T 1 , T 2 fMMCOOE) 

ARGUMENTS MMCOOE 0 MEASUREMENT MODEL CODE 

T1 I PRESENT TRAJECTORY TIME 

T2 0 trajectory TIME AT WHICH THE NEXT 

MEASUREMENT OCCURS 


SUBROUTINES SUPPORTEOI 
LOCAL SYMBOLS I M 

INDEX 

ERRAN 




COMMON USEOS 

MCNTR 

MCDOE 

NMN 

TMN 


SCHED Flow Chart 

MCNTRP 

MC00E2 

NMNP 

TMNl 

MCOOEl 

TMN2 


C enteQ 


^ MCNTR 

X 




loop on M 
greater or 

to find TMN(M) 
equal to T1 


' 

T2 - T> 
MMC0DE 

IN(M) 

- MC0OE(H) 


^ RETURN ^ 


4-167 




setevn-a 


SUPROUTINP bFTEVN 

PURPOSE: COMTROL COMPUT f, T I QK.S COMMQM TO ALL cVEnTS IN ERRAM 

CALLING SE'JUFMCE: CaLL St TtVN ( p I , TE NCOUt ) 


arguments: RI STATf OF TRIMI 

TF\/'J TTMt OF E\/ENT 
NCOOe TYPE OF EVENT 


SUhROUIINFS REQUIRED: 

CORREI. nvNO 

E JGhY 

GNA VM 

GPPINT 

MEAN 

' 

NjM RSIM 

STMPR 

STVCPR 



LOCAL symbols: EGVCT 

array of eigenvectors 




EGYL 

CORRESPONDING 

ARRAY OF 

eigenvalues 


E X T I J 

INTERmEOI ATE 

variable 




OUT 

SQUARct rtOOTS 

OF eigenvalues 



PEIG 

INTEPPEOl ATE 

ARRAY 




RF 

STATE VECTOP 

AT event 

T 1 ME 



VE IG 

matrix to PE 

NT AGONAL I4FO 



COMMON rOMRUTED/USED : 

1,'FLTm TpTmI 

XF 




COMMON USED: CXSU 

xcsv xcu 

' xcv 

CXXS 

FOP 

FOV 

TpRT 

TSTMC PS 

P 

0 

UO 

VO 

XLAH 

IGFN GO 

GV 

GCXW 

GCXSW- 

GP 

GCXXS 

GCXU GCXV 

GPS 

GCXSU 

GCXSV 

QPR 

RPR 

“^xT EXST 
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SETEVN Analysis 

Before executing any ev^t in the error analysis/generalized co- 
variance analysis program subroutine SETEVN is called to perform 
a series of computations that are common to all events. Subroutine 
SETEVN computes the targeted nominal trajectory at t^ and propagates 

the assumed and actual knowledge covariance partitions at tj^ ^ — 

the time of the previous event or measurement — forward to time 
t. using the propagation equations found in subroutine GNAVM. The 


actual estimation error means are also propagated forward to t 
using the propagation , equations found in subroutine MEAN. 




For any event other than a prediction event, subroutine SETEVN also 
computes eigenvalues, eigenvectors . of the 

position and velocity partitions of the assumed and actual know- 
ledge covariance at t . 

J 
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SETEVN Flow Chare 


( ENTER j 
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seti-a 


SUBROUTINE SETl 

PURPOSE; TO initialize FLAGS FOR USE BY INTEGRATION ROUTINES 
CALLING SEUUENCF: CALL SETl(Y»OJ) 

ARGUMENTS; 


Y I INITIAL 6 ELEMENT STATE VECTOR 

UJ I JULIAN date CORRESPONDING TO STATE VECTOR Y 


subroutines REQUIRED: 
CALJUL 
DShIFT 
OZERO 


COMMON USED: 

NbM 

N80PTH 

SMU 

IBURN 

COMMON COMPUTED; 

YMDIC 

HMSIC 

T 

X 

xo 

iters 

NCONOT 

nsectn 

SECTIM 


retro 


XDD 

ICENT 

XV 

6M 

xvo 

INO 

XVOO 

H 

NB 

I SUN 

NEQ 

NOFC 

NOCOWL 

NTSEQS 

NSTR 

CETOL 


« 
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subroutine Shift 

PURPOSE! TO SHIFT A DOUBLE PRECISION ARRAY TO ANOTHER LOCATION 

CALLING SEQUENCE! CALL SHIFT (AiNtK. BOUT) 

ARGUMENTS! AIN ARRAY TO BE SHIFTED 

K . NUMBER OF CONTIGUOUS VALUES TO BE SHIFTED 
BOUT RECEIVING ARRAY 


SUBROUTINE SL0020 
PURPOSE: TO MINIMIZE F(X) 

CALLING SEQUENCE! CALL SLO 020 ( X 9 Y . TS I 
ARGUMENTS! 

X I 3 VECTOR OF VALUES OF X 

Y 13 VECTOR OF F(X) VALUES CORRESPONDING TO VECTOR X 

TS 0 PROJECTED VALUE OF X FOR WHICH F(X) IS MINIMUM 

SUBROUTINES REQUIRED: 

NONE 


« 
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SURR0UTIt>3E S7APRL 

PURPOSE!! ¥0 COMPUTE THE PftRTEAL DERIVATIVES OF STATION LOCATION 

ERRORS s 


calling seoucnces 

CALL 

STAPRL(ALf AL0N,ALAT9P'AT2,VEC«PAr 

ARGUHENT9 AL 

I 

ALTITUDE OF THE STATION 

ALAT 

I 

LATITUDE OF THE STATION 

ALON 

I 

LONGITUDE OF THE STATION 

PA 

0 

PARTIAL OF STATION POSITION AND VELOCIT 
WITH RESPECT TO ALTITUDE, LATITUDE AND 
LONGITUDE 

PAT2 

1 

LONGITUDE ♦ ONEGA* CCURRENT TIME-LAUNCH 
TIME) 

VEC 


UNUSED 

SUBROUTINES SUPPORTED! 

TRAKS TRAKM 

LOCAt. SYMBOLS S Ci 


SINE OF latitude 

G2' 


COSINE OF LATITUDE 

G3 


SINECPHI ♦OMEGAa-UNIWT) ) 

G4 


COSINECPHI ♦OHEGAJT-UNIVT) S 
WHERE PHI *LOMC1TUOE 

OHEGAaEARTH NOTATION RATE 
T *TIME 

UNIVT=UNIVERSAL TIME 

65 


SINE OF OBLIQUITY OF EARTH 

G6 


COSINE OF OBLIQUITY OF EARTH 

OMEG 

OMEGA IN PROPER UNITS 

COMMON IISEOfi 


EPS OMEGA TM 
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STAPRL Analysis 

The ecliptic components of the position and velocity of a tracking station 
relative to the Earth are related to station location parameters R, 0, and 
^ through the following sat of equations; 

X a R COB 0 cos G 
e 

Yg R cos 0 cos sin G 4 R eln 0 sin € 

-R cos 0 sin € sin G -• R sin 0 cos £ 

X “ - tJR cos 0 sin G 
8 

Y ■ R cos 0 cos € COB 0 
z m COS 0 sin € COB V, 

B 

where G ■ 0 + W (t - T) , and T Is the universal time at some epoch(usually 
launch time) . 

Subroutine SZA?RL computes the negative of the partlals of the previous 
quantities with respect to the station location parameters R, 0, and , 

These partlals are summarized below: 

- cos 0 cos G 


R Bln 0 cos G 


R cos 0 eln G 


- 1^8 In € sin 0 + cos c cos 0 Bin ^ 


R cos € sln0 slnG -RsinC cos 0 


°R COB £ coo Q cos G 


sin i COS ® Bit* 0 " COB € sin 0 


aXg 

«R 

iia 

ae 

dX 

a 

3^ 


©R 

5Y. 

T© 

6Y 


i . 


3 


8 


dZ 
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STAPFL-2 






d'Y 

e 

§Y 

8 

b 0 



dR 


dZ 

a 

©e 



» - [r sin € aln 6 aln G + R coa € coe 9^ 


R Bin € cos 6 coa G 


“ (J coe 9 sin G 


“ - cjR aln 9 sin G 


■■ cob 9 cos G 


■> o»coa 0 coa e cos G 


(jR coa € sin 9 cos G 


B ^^fR cos £ cos 9 sin G 


Uf Bln € coe 9 cos G 


« -u'Sl sin £ sin 9 cos G 


■■ - wR Bln £ cos 9 sin G 
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STEAPE-A 


MAIN PKOGRAM STFAPE 

PURPOSE.; TO CONTROL TmE ERROR ANALVSIS MODE OF STEAP 
CALLING sequence: NONE 

ARGUMENTS: NONE 

SUBROUTINES REQUIRED; DATA ERRAN 
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STMPR-A 


SUaSOUTlNE SfMPR 

PUffPOSE® TO PRINT OUT THE TRANSPOSES OF THE STATE TRANSITION 
HATRIS PARTITIONS PHI 9 TXXS, TXW, ANO TXU OVER AN 
ARBITRARY INTERVAL OF TINE. 


CALLING SEQUENCE t 

CALL 

SIMPR(TRTMljTRTM2} 


ARGUMENT 8 

TRIMl 

I 

TIME AT BEGINNING OF INTERVAL OVER 
STATE TRANSITION MATRIX PARTITIONS 
BEEN COMPUTED 

WHICH 

HAVE 


TRIM 2 

I 

TIME AT END OF INTERVAL OVER WHICH 
transition MATRIX PARTITIONS HAVE 
COMPUTED 

state 

BEEN 


SUBROUTINES SUPPORTEOt PRINT 4 

SETEVS 

GUISIM 

GUISS 

PRESIM 

PRINT3 

SETEVN 

GUIOH 

GUIO 

PRED 

PROBE 

PROBES 





COMMON USED* 

NOIMl 

N0IM2 

PHI 

TXU 

TXXS 


XLAB 

XSL 

XU 




NDIM4 

TXW 
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STVCPR-A 


subroutine stvcpr 

PURPOSES print the EPOCH? T» AND PRINT THE STATE VECTOR XI IN SEVERAL 
coordinate systems (GEOCENTRIC ECLIPTIC? HELIOCENTRIC ECLIP- 
TIC? ROTATING 8ARYCENTRIC? ROTATING L1=CENTRIC« AND ROTATING 
La-CENTRIC) 

CALLING SE(SUeNCES CALL STVCPR(Xl,T) 


ARGUMENTSs XI STATE VECTOR 

T EPOCH OF state VECTOR 


LOCAL 

VARIABLES; 

AUX 

intermediate vector 



EL 

ROTATION VECTOR 



R 

VECTOR TO BE PRINTED 



RB 

SUN-TO-BARYCENTER DISTANCE 



RM 

RADIUS magnitude TO BE PRINTED 



RX 

intermediate value 



theta 

rotation matrix 



VM 

VELOCITY magnitude TO BE PRINTED 



VP 

velocity due to rotating frame of reference 



XB 

HELIOCENTRIC STATE VECTOR OF 0ARYCENTER 



XM 

HELIOCENTRIC STATE VECTOR OF MOON 



XP 

HELIOCENTRIC STATE VECTOR OF EARTH 

COMMON 

USED: 

DATEJ 

MUPLAN TRTMB 
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STVCPR Analysis 


Subroutine STVCPR prints the state vector XI at epoch T in several coordinate 
systems: geocentric ecliptic, heliocentric ecliptic, rotating barycentric, 
rotating Ll-centric and rotating L2-centric. The state vector XI is read 
from the GTDS (Cowell) file in the geocentric ecliptic coordinate system, and 
so the first print involves no transformation. The earth's heliocentric 
ecliptic state is then added to XI for the second print, XI in heliocentric 
ecliptic state. 



Next the moon's heliocentric ecliptic state is obtained, and the barycenter's 
state is computed 

E M 

The rotating barycentric frame has the same Z-axis as the ecliptic, but 
rotates so that the X'-axis is always the barycenter-sun line, A matrix 
is defined to effect the rotation of the position and velocity components. 




Then the velocity change due to the rotation of the frame and to the 
spacecraft position in the frame is added to the velocity, and the third 
print is the rotating barycentric representation of XI. 
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For the last two representations, the origin of the coordinate system is 
simply translated back and forth along the X' -axis to the Lagrangian LI 
and L2 points respectively. 


SUN 


MOON 



EARTH 
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SYmTRK-A 


SUBROUTINE SYMTRK 

PURPOSE: Symmetrize a square matrix 

CALLING sequence: CALL S YMTRK ( ARRAY , M ) 

arguments: array NAME OF THE M-8Y-M MATRIX 

M DIMENSION OF THE SQUARE ARRAY 

subroutines REQUIRED: NONE 

COMMON USED: HALF 


subroutine SYMTHZ 

PURPOSE: fill IN THE UPPER-RIGHT TRIANGLE OF A SYMMETRIC SQUARE 

MATRIX whose lower-left TRIANGLE WAS INPUT 

CALLING sequence: CALL SYMTRZ ( P »K » N) 

arguments: P square MATRIX TO BE COMPLETED (N-HY-N) 

K NUMBER OF ROWS IN USE 
N actual DIMENSION OF MATRIX 

subroutines required: NONE 


SUBROUTINE TfME 


PURPOSE: convert A TIME IN SECONDS INTO DAYS* HOURS. MINUTES* SECONDS 

CALLING sequence: CALL T IME ( TSEC » JDA Y . JHR » JM I N , XSEC ) 

ARGUMENTS: TSEC TIME IN SECONDS 


jday 

NUMBER 

OF 

DAYS 

JHR 

NUMBER 

OF 

HOURS 

JMIN 

NUMBER 

OF 

MINUTES 

XSEC 

NUMBER 

OF 

seconds 
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TINE-A 


SUBROUTINE TINE 


PURPOSES TO COMPUTE THE JULIAN DATE, EPOCH 1900, FROM THE 

CALENDAR DATE OR TO COMPUTE THE CALENDAR DATE FROM THE 
JULIAN DATE. 


CALLING SEQUENCE* 

CALL 

TINEtOAY, IYR,MO, lOAY ,IHR,M1N,SEC, ICOOE) 

ARGUMENTS DAY 

I/O 

JULIAN DATE, EPOCH 1900 

lYR . 

O/I 

CALENDAR YEAR 

MO 

O/I 

CALENDAR HONTH 

lOAY 

O/I 

CALENDAR DAY 

IHR 

O/I 

HOUR OF THE DAY 

MIN 

O/I 

MINUTE OF HOUR 

SEC 

O/I 

FRACTIONAL SECONDS 

ICOOE 

I 

operational MODE 


r-lvWI. 

= t, INDICATES THE JULIAN DATE IS INPUT, 
CALENDAR DATE IS OUTPUT 


= 0, INDICATES the CALENDAR DATE IS INPUT, 
JULIAN DATE IS OUTPUT 


SUBROUTINES SUPPORTED 1 DATA 


SUBROUTINES REQU?RE08 
LOCAL SYMBOLS* I A 

IB 

IP 

ID 

IR 

IS 

XT 

lU 


NONE 

NUMBER OF CENTURIES 
YEARS IN PRESENT CENTURY 

NUMBER OF MONTH I BASED ON MARCH AS NUMBER 
2EROJ 

NUMBER OF YEARS 

NUMBER OF CENTURIES DIVIDED BY 4 

NUMBER OF YEARS SINCE LAST 400 YEAR 
SECTION began 

NUMBER OF LEAP YEARS IN PRESENT CENTURY 
NUMEER OF YEARS SINCE LAST LEAP YEAR 


number of days in LAST YEAR 
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IX INTERMEDIATE INTEGER 

j Intermediate integer 

JD number of days in JULIAN DATE 

P JULIAN DATE 

R FRACTIONAL PORTION OF DAY IN JULIAN DATE 
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TRAKM-A 


SUBROUTINF TraKM 

PURPOSE; compute the observation matrix and its augmentations 


CALLING sequence; 


Call trakm jrv*itrk,nr) 


arguments; 


RV state vector at TIME OF MEASUREMENT 
ITRK TYPE OF ^MEASUREMENT 

NP number of rows in observ.ation matrix 

(I.E.. THE DIMENSION OF THE MEASUREMENT) 


LOCAL SYMBOLS; 


AL 

ALAT 

ALON 

CE 

COAL 

CP 

DEL 

62 

63 

G5 

GELS 

lA 

1C 

ICO 

lENO 

IK 

IL 

IR 

NA 

NC 

PA 

PAR 

PATl 

PAT2 

RAONEP 

RF 

HFMAG 

RRATE 

RS 

SE 

SIAL 

SP 

ZTAR 


RADIUS OF STATION lA 

latitude of station IA 
longitude of station IA 

COSINE OF EARTH OBLIQUITY 
COSINE OF STAR-PLANET ANGLE 
COSIME OF PAT? 

time IN DAYS FROM INITIAL TIME ON FILE 

Y-COMPONENT» EQUATORIAL LOCATION* STATION lA 

Z-COMPONENT* equatorial LOCATION* STATION lA 

YOOT-COMPONENT* EQUATORIAL LOCATION* STATION lA 

GEOCENTRIC EQUATORIAL LOCATION OF STATION IA 
INDEX OF STATION MAKING MEASUREMENT 
COLUMN NUMBER IN MATRIX AL 

INDEX TO LOCATE AUGMENTATION PARAMETERS FOR IA 
LAST AUGMENTATION PARAMETER TO BE CHECKED 
COLUMN NUMBER IN MATRIX AM 
COLUMN NUMBER IN MATRIX AN 
reference index* STAR-PLANET ANGLES 
STAR-PLANET ANGLE CURRENTLY BEING COMPUTED 
LAST STAR-PLANET ANGLE TO BE COMPUTED 
ARRAY OF PARTIALS D ( RF ) /O ( GELS ) 

ARRAY OF PARTIALS D <RF ) /D ( R-PLANET ) 

INTERMEDIATE VARIABLE 

HOUR angle of STATION lA 

RADIUS OF EPhEMERIS PLANET (EARTH) 

state VECTOR WRT EARTH OR WRT STATION lA 

MAGNITUDE OF RF 

RANGE-PATE . - 

SPIN RADIUS OF STATION IA 

SINE OF .EARTH OBLIQUITY 

SINE OF STAR-PLANET ANGLE 

SINE OF PAT2 

UNIT VECTOR* OIRECTION OF STAR NA 
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TRAKM-B 


COMMON computed; al am an G H 

UL IN common block STM IS CALLED AAL IN TRAKM> 


COMMON USED; aLNGTH DELTM EPS 

SAL Slat slon 

UST VST WST 


lAUGIN OMEGA RAD RADIUS 

TM TRTMB TRTMl UNIVTM 
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TRAKM Analysis 


Subroutine TRAKM computes all observation matrix partitions for the 
measurement type indicated by ITRK. The number of rows, NR, in the 
measurement and the observation matrix partitions is also computed. 

The linearized observation equation can be written as 


y = Hx + Mx + Gu + Lv + Nw 
s 

where y is the observable, x is the spacecraft state, and x^, u, 

V, and w are solve-for, dynamic consider, measurement consider, and 
ignore parameter vectors, respectively. The function of subroutine 
TRAKM is to compute the observation matrix partitions H, M, G, L, 
and N, which indicate the sensitivity of the observable v to changes 
in X, x^, u, V, and w, respectively, in the error analysis/generalized 

covariance analysis program. The matrix N is computed only for a 
generalized covariance analysis. 

In the remainder of this section the measurement equation and all 
partial derivatives required to construct the H, M, G, and L observa- 
tion matrix partitions will be summarized for each measurement type. 


A. Range Measurement . 

A. range measurement has form 

(J = ^ ( X, R, ©, 0, t) 

where R, 6, and 0 are the radius, latitude, and longitude of the 
relevant tracking station. 
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More explicitly, 

p. [(x.;<^ 


X )^ + (Y - Y 
S E 


Y + (Z - Z 
S E 



where X, Y, Z = inertial position components of spacecraft 

X , Y , Z = inertial position components of Earth 
E’ E E 

X , Y . Z = station position components relative to Earth. 
S’ S S 


X , Y , and Z arc related to R, 9, and as follows: 
S S S 


X = R cos 9 cos G 
S 

Yg = R cos 9 cos € sin G + R sin 9 sin€ 

Z = -R cos 9 sin € sin G + R sin 9 cos € 
S 

where ^ is the obliquity of the Earth, and 

G = 9 + GWi 


where GHA 

is 

the Greenwich hour angle at 

time 

t . 


Partials 

of 

p wi th 

respect 

to spacecraft 

state are 

given 1 

dp 
a X 

= 


\ 



_ap 

ax 

= 0 


dP 

d Y 

= 

■ 

Y - 

E 

V 


dp 

dY 

= 0 


dp 
d Z 

= 


^E- 



dP 

di. 

= 0 


Partials 

of 

P with 

respect 

to R, 9, and 

9 

are given by 

d P 
d R 

= 

± 9 - . 

^R 

+ 

dP 

5Yg 5R 

+ ■ 

dP 

5R 

dp 


dp 



dp 5Yg 


dP 


59 


-s' 

50 


5Y 59 

S 

T 

dZ 

S 

50 

dp 


dP 

axg 


ap ' 


ap 


d9 



d9 

•T 

a Yg a 9 

T 

’ 

<90 
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where 


d P 


dp 

d X 


dp 

d'i„ 


dp 

d'i 


3P 

dz_ 


dp 

dZ 


and the ne^ntives of the partials of , and with respect to 

K, 9, and are suiniiiari /.ed in tlie subroutine STAPRL unalysi.s. 

h, Ran^', e-rate nieasiireiiicnt p . 

A range-rate rneasurement has form 

P ^ P iX, R, G, 0, t) 

where all arguments have been defined previously. More explicitly, 



• 

D ■ 

= P 

v-Px^PyPy^ 

Pz Pz 





H 


P 





where 

P 

X 

= 

X - 

E S 

K 

= ■ 

X - 

X - X 
E S 




Y - 

E S 

P 

y 

= 

Y - 

• « 

Y - Y 
E S 


P. 

= 

Z 

E S 

i 

= 

Z - 

Z - Z 
E S 

Xs, Ys, and 


are related to 

R, 0, and 

9 

as 

fo 1 lows 


X 

s 

• 

= 

- R cos 0 s in G 







U) R cos 0 cos 

€ cos G 





z 

s 

= 

-WR cos 0 sin f cos G 




where (p 

is the rotational rate 

of the Earth. 



Partials 

9 

of p with respect to 

spacecraft 

state are given 

dP 


9 

Px 

Px P 

dP 


Px 


ax 


p 

P2 

dx 


P 


dP 


9 

Pjl 

. PyA 





d'i 


P 

p2 

dY 


P 


dP 


Pz 

. p-iA 

dk 

B 



d z 


p 

p2 

d z 


P 
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The partial 

of p with 

respect to R is given by 

dp 

dP dh ^ 

di> „ ap i!s 

d R 

6>Xg d 

R 

aYg aR azg aR 


dP 

ax's 

1 ap ^ dP d\ 


ax 

s 

■ a^ 

a^s a^s ’ a^ 

where 




dP 

1 



dx^ 

dx 

6 CC • 


s 






etc . 

f 


^x 




# • « 

The negatives of the partials of X^, Yg, Zg, Xg, Yg, and 2g with respect 

to R, G, and 0 are summarized in the subroutine STAPRL analysis. Partials 
of p with respect to 9 and 0 are treated similarly. 

C. Star-planet angle measurement Ct , 

A star-planet angle measurement has form Ct 

ce^cs ( X, a, e, i, . M) 

where a, e, i, (j , and M are the standard set of target planet 
orbital elements. 


If we differ p = (p P >P position of the target planet 

relative to the spacecraft and (u, v, w) to be the direction cosines of the 
relevant star, then 


a = 


cos 


-1 


-i-(up + vp + VI P ) 


where 


P = X - X, p = Y - Y, p = Z - 2, 
X p y p z p 


and (Xp, Yp,, 2^) represent the position coordinates of the target planet. 
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Partlals of Cl with respect to spacecraft state are given by 


da ^ ^ '' / _u_ 

d X sinOJ \ ^ 

1 f ^ 

Sy sin cs \ ^ 

do 

6 Z 

where 

sin Of 

D. Apparent 



da 

p'^ / 

'dk 

Py COS Of ^ 

del 

p2 ) 

6y 

P^ cos Of \ 

da 

p2 / 

dz 


sin Of 


planet diameter measurement 0 


= 0 


= 0 


= 0 


An apparent planet diameter measurement has form 
0 = 0{X, a, e, i, Q , (J , 

where all arguments have been defined previously. 


Defining P = (P , P , P ) to be the position of the target planet relative 
X y z 

to the spacecraft and Rp to be the radius of the target planet, die apparent 
planet diameter can then be written as 


0 = 2 sin'l 

Partlals of 0 with respect to spacecraft state are given by 




2 R P 
P ^ 


o 

u 

dx 

P^ 

■ V ] 


dk 

1 



2 R P* 
p y 


ii . 0 

d Y 

P^ 




d0 


2 R P 
P z 


M = 0 

dz 


■ V : 


dz 


1 
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TRJTRY-a 


SUBROUTINE’ TRJTRY (PRELIM ENTRY POINT) 

PURPOSE; dummy link WITH NON HALO ORBIT OPTIONS 
CALLING sequence; CALL TRJTRY 


subroutine trnsps 

purpose; to form the transpose of a matrix 

CALLING sequence; CALL TRNSPS (A»B»M.N) 

arguments: 

A I matrix to be transposed 

B 0 RECEIVING MATRIX 

M I NUMBER OF ROWS IN A AND COLUMNS IN B 

N I NUMBER OF COLUMNS IN A AND ROWS IN B 

SUBROUTINES REQUIRED: 

NONE 


SUBROUTINE ZERMAT 

PURPOSE: TO ZERO OUT A MATRIX 

CALLING SEUUeNCE: CALL ZERMaT ( ARR»NR»NC) 

ARGUMENTS: ARR ARRAY TO BE ZEROED OUT 
NR NUMBER OF ROWS 
NC number OF COLUMNS 
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